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The Distribution of Deformation in Lead Fatigued in vacuot 


By K. U. Snowpen}t 
The Baillieu Laboratory, University of Melbourne, Victoria, Australia 


[Received J uly 27, 1960 ; and in revised form October 20, 1960] 


ABSTRACT 


Observations have been made on the nature of the deformation, hardening, 
and recrystallization in polycrystalline specimens of high-purity lead fatigued 
in reverse-plane bending in vacuo (approx. 5x 10~° or 2x 10-§ mm Hg). 

Fatigue straining induced slip, and grain boundaries to migrate to more 
stable positions at + 45 degrees to the principal stress axis. 

The early distribution of slip-trace angles was in good agreement with the 
angular distribution which Hedgepeth (1952) predicted on a basis of the 
maximum shear stress criterion and by neglecting the slip and grain-boundary 
interaction. Evidence for this latter interaction was found as hard zones up to 
109 microns wide near grain boundaries. These zones were the sites for 
recrystallization and grain growth. 


§ 1. LyrropucTION 


THE fatigue strengths of a number of metals have long been known to 
be affected by the presence of the atmosphere during cyclic straining 
(Haigh and Jones 1930, Gough and Sopwith 1932, 1935). More recently, 
in postulating and evaluating theories of fatigue in terms of fundamental 
concepts, such as dislocations and lattice defects, it has become important 
to separate the part played by the alternating strain from that by the 
attack of atmospheric gases. Moreover, Wood (1959), and Segal and 
Partridge (1959) have shown that the difference between unidirectional 
and fatigue deformation becomes more distinct at small alternating 
strains, and Wood considers that it is the deformation at small alternating 
strains which is typical of fatigue damage. It therefore seems likely that 
under these conditions of strain with the associated long life, the attack 
by atmospheric gases may have sufficient time to contribute significantly 
to the fatigue damage. 

Recent investigators have examined the effect of atmosphere on fatigue 
life and deformation (Thompson et al. 1956, Wadsworth and Hutchings 
1958, Snowden and Greenwood 1958b, and Hempel 1959). These have 
shown that the lives of copper, aluminium and lead are affected by the 
presence of the atmosphere whereas the lives of gold and steel are not. 
The observations on copper and aluminium suggest that the initial crack 
formation is much the same in air as in vacuum but that the rate of crack 
growth is accelerated by the presence of the atmosphere. On the other 
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hand, the early observations of Gough and Sopwith (1935) indicate that 
the manner of fatigue failure of lead in air is different from that in vacuum. 
These observations are supported by more recent work (Snowden and 
Greenwood 1958b) which shows that lead fatigued in vacuum (about 
5x 10-mm Hg) is more severely distorted than lead fatigued in air for 
the same number of cycles. 

High-purity lead specimens have been fatigued at an air pressure of 
either ~5 x 10-8 or ~2 x 10-* mm Hg, and the nature of the deformation 
and hardening is examined in more detail below. It has been shown 
(Snowden 1959, 1961) that at about 5 x 10-? mm Hg, the fatigue deformation 
and life of lead at 500 cycles per min is no longer affected by further reduction 
of air pressure. The observed changes are therefore thought to be due 
almost entirely to the alternating strain. It is hoped to publish observa- 
tions on the influence of air pressure on fatigue life at a later date. 


§ 2. EXPERIMENTAL METHOD 


Specimens of the shape shown in fig. 1 were tested in a reverse-plane 
bending machine at the rate of 500 cycles per min. The specimens were 
prepared from lead strip (99-999% pure), 0-160 in. thick, rolled from 
rectangular ingots | x 1 x 9 in. cast at 400°c. Prior to fatigue, specimens 
were cleaned with benzene, placed on flat steel plates, annealed in air 
at 100°c for 24 hrs, and chemically polished in a 70:30 mixture of glacial 
acetic acid hydrogen peroxide (Worner and Worner 1940). The specimen 
shape was designed to produce a uniform macroscopic strain over the 
test length, which is indicated by the shaded area in fig. 1. 


Fig. 1 


Specimen shape. 


For the present observations, specimens were fatigued either at ‘high’ 
or ‘low’ alternating strains which were arbitrarily chosen as: <0-140 oS 
to 20:075% for the high strains and 0-034% for the low strain. The 
endurance corresponding to these strains was from 1:50 108 to 
25 x 10® cycles, and > 2x 107 cycles respectively. 


§ 3. EXPERIMENTAL RESULTS 
3.1. Grain-boundary Migration and Slip Formation 


As described in a previous paper (Snowden 1958), the initial effect of 
the alternating strain was shown by the appearance of slip traces and 
grain-boundary migration. The rate of migration increased with 
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increasing strain and was highest in the early part of the life, ie. in the 
first 104 to 10° cycles. At the low strain, figs. 4 and 5, the grain structure 
tended to stabilize, except for the sporadic migration mentioned below. 
Boundary displacements of one or two hundred microns, such as those 
shown in fig. 4, were found to be common. However, some boundaries 
including coherent twin boundaries, did not move. 

A striking feature of the migration at both high and low strains was 
that boundaries tended to migrate to positions at +45 degrees to the 
specimen axis (figs. 4 and 5) as in lead fatigued in air (Snowden and 
Greenwood 1958a). However, in specimens fatigued in vacuum there 
was no sign of the intercrystalline cracks that formed in the specimens 
fatigued in air. This type of boundary migration resulted in the formation 
of an orthogonal grain structure which was better developed at the higher 
strains, fig. 27. 


Fig. 3 
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Theoretical and experimental slip-angle distribution for a specimen fatigued. 
to 10° cycles at a strain of +0-092%. Theoretical distribution, curve 2, 
after Hedgepeth (1952). 


In view of the above tendency for boundaries to align themselves at. 
+45 degrees to the specimen axis, the angular distribution of slip traces 
was examined to see if a similar relation existed for slip. The angular 
distribution of 289 traces was determined on the surface of a specimen 
fatigued to 10° cycles at.a strain of +0:092%. The result is shown in 
fig. 3, curve 1, in the form of a cumulative frequency diagram plotted 
at 5 degree intervals. The ordinate N(6) is the sum of frequencies of 
angles up to and including a given angle 0, where 0 is the angle between the 
slip trace and the transverse specimen axis. The curve shows that 
approximately 65% of traces made angles between 30 and 60 degrees with 
the transverse axis of the specimen. This suggests that there is a tendency 
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for traces to occur at angles near the +45 degree positions although the 
trend is not as clear-cut as that shown by the boundaries at the higher 
strains. 

A further feature of curve | is that no slip traces were found which made 
angles greater than 70 degrees with the transverse axis of the specimen. 
This feature and the shape of curve | are in accord with the theoretical 
slip-angle distribution, curve 2, which Hedgepeth (1952) predicted for 
unidirectional straining on the basis of the maximum shear-stress criterion 
and by neglecting slip and grain-boundary interaction. The agreement 
suggests that the crystal aggregate deforms initially as if each crystal were 
free to deform independently of the other crystals. However, Taylor 
(1938) and others have pointed out that, for unidirectional deformation 
(at high strain rates and normal temperature), grains remain continuous 
at the boundary with their neighbours during deformation, i.e. there is 
an absence of cavities. Evidence for this latter restriction in fatigue 
at the low strain is provided by the series of micrographs shown in figs. 4 
to 6, in which the initial slip distribution is modified near the grain 
boundaries, although the slip systems (indicated by double-headed 
arrows) operating in the central portions of grains remain unchanged. 
In order to bring out the structural changes more clearly, the specimen 
was repolished after various amounts of fatigue, thereby removing previous 
fatigue traces and observing the distribution of deformation due to 
subsequent fatigue. Figure 4 shows the area after 5 x 10* cycles. Previous 
observations on this area showed that the boundary of the central grain 
had migrated from the trace A to the final position B. The specimen was 
repolished at 7-5 x 10° cycles, given 4:55 x 10® cycles further, i.e. a total 
of 5-3 x 108 cycles, and photographed again, fig. 5. Figure 5 shows that 
after repolishing, the deformed area which was produced by subsequent 
fatigue, tended to stop short of the grain boundaries. Approximately 
80°% of boundaries showed zones of low slip density of up to 100 or so 
microns wide. The slip density in the boundary zones tended to increase 
slightly with increasing number of cycles; however, after an intermediate 
polish at 5-3 10® cycles and further fatigue the zones reformed and 
tended to be wider than before repolishing. This can be seen by comparing 
figs. 5 and 6. 

Finally, examination of the area at 8 x 10° cycles showed that the zones 
of low-slip density underwent deformation on a new slip system not 
previously evident in the adjoining grains, fig. 7. This behaviour was the 
result of local recrystallization which produced a new grain along the 
original grain-boundary zones. The strip of new material readily deformed 
and built up its own zones of low-slip density as was shown by repolishing 
and further fatigue, fig. 8. The relative positions of boundaries before 
and after recrystallization are illustrated in fig. 9 which shows that the 
new grain in fig. 8 straddled the original orthogonal grain boundary. 

Although the zones of low-slip density showed a steady development 
at low strains, the onset of localized recrystallization occurred spon- 
taneously and irregularly after lengthy periods of cyclic straining. At the 
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higher strains, however, similar structural instability near grain boundaries 
was detected earlier and at a strain of +0-092°% was even found during 
the initial period of rapid grain-boundary migration, i.e. in the first 104 
to 10° cycles. An example of this is shown in figs. 10 and 11 in which the 
migrating boundary moved from positions A to B between 10% and 
2x 10* cycles. In this case the specimen was not given intermediate 
polishes. 


Fig. 9 


NEW ORIGINAL POSITION 
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Diagram shows position of recrystallized grain (figs. 7 and 8) with respect to 
original grain boundary. 


Recent observations by de Fouquet (1958) have suggested that an 
analogous process of grain growth along the intergranular regions may 
occur during the fatigue of iron in argon at 700°c. This kind of effect 
has also been found under other conditions of stressing by Brinson and 
Hargreaves (1958). 

§ 4. Discussion 


The observations described above throw some light on the nature of the 
slip deformation and of the process of work -hardening near grain boundaries 
in fatigue. The early stages of the deformation suggest that dislocations 
are set in motion on octahedral planes of high resolved shear stress, i.e. on 
the so-called ‘natural’ systems (Hargreaves 1956), as in unidirectional 
straining. During the first half cycle, some of these dislocations would 
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interact with grain boundaries and form pile-ups there or be transmitted 
to adjoining grains according to the orientation difference. However, if 
interactions were limited to simple pile-ups, then during stress reversal 
these dislocations would tend to run back on their glide plane so that the 
hardening would be dispersed after each half cycle. Evidence for this 
type of reversed movement during a partly-reversed stress cycle has 
been found by Takamura and Mura (1958) in the case of an a-brass 
bicrystal. Nevertheless, the present results indicate that when the 
strain is cycled for a large number of times, some of the boundary hardness, 
as indicated by the decreased tendency of the boundary zone to undergo 
slip after repolishing and fatigue, persists and in fact tends to extend away 
from the boundary (figs. 5 and 6). 

In this case, the pile-ups at a boundary are thought to activate, near 
the boundary in the adjoining grain, extra slip systems which interact 
and form barriers to slip on the natural systems. Examples of slip on 
multiple systems near grain boundaries in fatigued lead have been 
described previously (Snowden 1958). Further, Berghezan and Fourdeux 
(1959) have observed the dislocation interaction near grain boundaries 
in aluminium during unidirectional straining. They found that multiple 
slip operated and resulted in an entangled network of dislocations which 
contained a high proportion of jogged screw dislocations. A similar process 
may occur during fatigue as the latter dislocations would be difficult to 
move (Read 1953) and could therefore harden the boundary zone. Gough 
et al. (1923, 1934, 1936) have observed similar boundary zones in fatigued 
aluminium and steel. 

The initial grain-boundary migration may be explained on the basis 
that, in general, different amounts of deformation are built up on each 
side of boundaries so that, as Forsyth (1953-54) and Liicke and Detert (1957) 
have proposed, the migration is produced by unequal numbers of dis- 
locations on opposite sides of boundaries. Such a mechanism would 
indicate that the diminished rate of migration after 104 to 10° cycles 
results from an almost complete balancing of forces across migrating 
boundaries. A further possible reason for the diminished rate is suggested 
by the tendency for boundaries to align themselves in the directions of 
maximum shear stress. This suggests that the diminished rate is associated 
with enhanced grain-boundary sliding at orientations in the shear 
directions, i.e. at +45 degrees to the specimen axis. 

The observed features of the boundary zones and the intergranular 
grain growth are in accord with the model for recrystallization of a cold- 
worked metal proposed by Gay et al. (1954). In this model the onset of 
recrystallization is associated with the tilting in excess of a limiting angle 
of a ‘particle’ in a region of high dislocation density. The present results 
suggest that the tilting process, which evidently occurs more rapidly at 
high strains (figs. 10 and 11), takes place in the heavily strained boundary 
regions where the driving force for growth would be provided by 
absorption of the disoriented material. 
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Some New Habit Features in Crystals of Long Chain Compounds+ 
Part I. Paraffins 


By A. KELLER 
H. H. Wills Physics Laboratory, University of Bristol, Bristol, 8 


[Received July 20, 1960] 


ABSTRACT 


Habits of paraffin crystals were examined with the optical and electron 
microscope. In particular striations along different crystallographic 
directions were noticed which could be produced systematically by heat 
treatment. It was found that such striations often corresponded to geo- 
metrically regular ridges. It is proposed that these ridges represent different 
phases within the crystals corresponding to structures of differing obliquities 
all based on the same subcell. Accordingly a variety of structures can be 
realized within the same crystal by different amounts of simple shear on 
different planes along the ¢ direction where the shear, while regular, need not 
be uniform. A number of other less regular morphological features are also. 
described. 


$1. InTRODUCTION 


CrysTaLs of paraffin’s have been studied extensively in the past by 
various methods (e.g. Miller 1927, Dawson and Vand 1951, Teare 1959). 
Our own interest in these crystals originates from the recently discovered 
similarity between the crystal habits of paraffins and polymers, poly- 
ethylene in particular (Keller 1959). Accordingly, crystals of polyethylene 
and those of paraffins can have identical habits. Further research on 
polyethylene crystals has revealed a number of additional morphological 
details. In view of the hitherto unsuspected close connection between 
paraffins and polymers it has become important to separate those features 
which can be explained on the basis of the crystallography of paraffins, 
from those where the long chain nature of the polymers has to be 
specifically invoked. Consequently, the question arises how far the 
more detailed morphological features in polymer crystals have corre- 
sponding counterparts in paraffins. In particular, we were interested 
in the origin of certain striations which frequently appeared in polymer 
single crystals. These were very elusive and there appeared to be little 
hope to clarify their nature by the continued study of polymers only. 
Therefore we turned to paraffins where evidence of similar morphological 
effects was found which appeared more readily interpretable. The 
present paper will deal with these investigations on paraffins, at this 
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stage with the morphological aspects, although it is hoped that these 
investigations will be supplemented later with diffraction results. In the 
following paper of this series the results obtained with paraffins will be 
utilized in the elucidation of morphological features in polymers. 


§ 2, EXPERIMENTAL 


2.1. Materials and Techniques 


The majority of the investigations were carried out on nC,,H,,. All 
the figures and illustrations in the present publication refer to this latter 
paraffin except one photograph (fig. 9+) which is of nCz>Hgs. The 
phenomena under discussion, however, are quite general and were found 
to occur also in all the other normal paraffins, namely, Cy,H,,, C5 9H¢o, 
Ca Hog, CazH og, Cz4H 9, handled by us in the course of this work. The 
paraffins were dissolved in low boiling petroleum ether (b.p. 40—60°c) 
or occasionally in benzene. On prolonged standing, or on cooling the 
originally warmed up solution, crystals precipitated forming a suspension. 

The crystals were deposited on a microscope slide by evaporation of the 
solvent and examined first by the optical microscope. Most of the crystals 
of interest for this work were thin enough to profit from phase contrast 
illumination. All optical microscopy was restricted to the powers of 
dry objectives as the crystals could not be immersed in a mounting 
medium. 

The crystals, selected under the optical microscope, could also be 
examined with the electron microscope either by precipitating them 
originally on a slide coated with carbon, or by carbon coating the crystals 
subsequently. If diffraction patterns were required or effects due to 
diffraction (Bragg fringes, etc.) were to be examined in the image, very 
weak beam intensities had to be used in order to minimize the effect of 
the beam. We found that the beam destroyed the crystallinity without 
affecting the crystal shapes, an effect which has already been known in 
the case of polythene (Agar et al. 1959). This restricts the working to 
lowest powers, unless the double condenser technique can be applied 
(see reference above), which was not possible on the instrument available 
for this work. 

Transmission electron microscopy was only possible on the thinnest 
crystals. In the case of thicker crystals replication techniques had to 
be resorted to. For this purpose the crystals were first metal shadowed, 
coated with a moderately thick layer of carbon and the films were floated 
off in water and picked up on grids. By bringing the films in contact 
with xylene the paraffins were dissolved and the replicas were ready for 
examination. 

In most cases the crystals were heat treated for these studies. This 
was done by placing the microscope slide with the crystals on them on 
a temperature gradient bar. The seemingly irreversible changes happening 
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at different temperatures could be examined subsequently after the 
removal of the slide from the bar. The changes on heating could algo 
be followed in situ with the help of a hot stage, both under the optical 
and electron microscope. Hot stage work in the electron microscope 
was necessarily restricted to transmission, hence to the thinnest crystals. 
Even so, it was made very difficult by the beam damage which was 
more severe at the elevated temperatures. 


2.2. Observations 

The crystals were the usual lozenge-shaped tablets ranging from 
macroscopic dimensions (1 mm or more) to a few microns). As already 
stated, some crystals grew while in suspension on standing or cooling. 
Further crystals formed when the solvent evaporated on the slide. The 
crystals grown in suspension were usually larger than the ones formed 
on evaporation. The acute angles of the lozenges could be either in the 
range of 75° or 67°. The former corresponds to the monoclinic, the 
latter to the orthorhombic variety (Dawson and Vand 1951, Schearer 
and Vand 1956, Teare 1959). Crystals of C,,H,, grown from benzene 
were always monoclinic while both orthorhombic and monoclinic crystals 
were obtained from petroleum ether which was the most frequently used 
mode of preparation for this work. It was found in the course of a 
systematic study of habits carried out by Baltd-Calleja in this laboratory 
(details to be reported later) that the crystals which grew while in 
suspension were always monoclinic while those which formed on 
evaporation were orthorhombic, consequently the two could be separated 
by filtration. Besides regular lozenges many other habit modifications 
like twins, dendrites, lath-shaped crystals, etc., could be observed, many 
of which have been reported in the literature earlier and will not be 
itemized any further here. 

Electron microscope examinations revealed that the crystals thickened 
through spiral terraces indicative of the screw dislocation mechanism 
in line with earlier findings (Dawson and Vand 1951). It may be 
noteworthy to point out that all transmission electron microscopy 
so far has been carried out on specimens, the crystallinity of which was 
destroyed by the electron beam. Whenever the damage could be 
avoided, by means already referred to, a variety of extinction effects 
were obtained in the image. Figure 1 shows one example. The crystals 
in this photograph being unshadowed, would have appeared practically 
structureless—as in fact they did become after electron bombardment— 
without the Bragg fringes. The Bragg fringes correspond to those parts 
of the crystals which are in a suitable orientation to diffract electrons. 
These electrons will be screened off by the microscope aperture and thus 
will be absent from the corresponding parts of the image. Consequently 
the existence of Bragg fringes implies that the structure within the 
crystal is not uniform. The fringes can be indexed, i.e. the reflections 
producing them can be determined by taking dark ground photographs 
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through chosen reflections or by the ghost image method (Agar et al. 
1959), and with this information it should be possible to map the lattice 
orientation within the whole crystal. It was found, however, that the 
fringe pattern was changing in the electron beam—in addition to the 
trend of disappearing gradually—consequently it could not be readily 
assessed how far the fringe patterns which could be photographed reflected 
the original characteristics of the lattice and how far they were conditioned 
by beam effects. For this reason this particular line was not pursued any 
further. 

When scanning extended areas of the preparations under the optical 
microscope occasionally straight lines could be seen in the crystals along 
one or the other lozenge diagonal. These were rare in C,,H,, but in some 
of the other paraffins, particularly in C,,H,, (studied by Balta-Calleja) 
lines along the short lozenge diagonals could be very frequent. It is the 
study of these lines which forms the main part of the present investigation. 
In the following some of the most typical observations will be enumerated. 

It was found that the lines in question could be produced systematically 
by heat treatment. If C,,H,, crystals as grown from petroleum ether 
were heated to temperatures between 45-—60°c lines which were readily 
visible under the optical microscope (fig. 2) formed along the long lozenge 
diagonals (to be referred to as a lines). In the region of 60°c and above, 
these lines receded and the crystals could become quite structureless 
optically. Above 65°c striations appear again but this time parallel to 
the short diagonal (to be referred to as b lines). There was a considerable 
variability as regards the temperatures. It could happen that the b 
striations appeared while some a lines were still present. As the melting 
point 76°c was approached the b striation pervaded the whole crystal. 
Near the melting point striations could also develop along other 
directions, e.g. parallel or occasionally approximately perpendicular to 
the lozenge faces ({110) and (130) directions by the orthorhombic 
subcell structure; see later) the former being illustrated by fig. 3. As 
seen, the lozenge in this photograph also contains b lines. 

The present account will be restricted to the striations of the kind just 
described. It will only be mentioned that in the immediate vicinity of 
the melting point even more drastic changes happened to the crystals 
which were still unfused. Thus the b striations, while becoming more 
pronounced, also became less regular and the crystals could buckle, 
fold or even crack, usually along their short diagonal. When the melting 
point was exceeded the crystals fused into blobs which, on cooling, 
recrystallized in a strongly striated form. It is hoped to illustrate these. 
effects in a later publication. ; 

The regular a or b striations appeared in the orthorhombic form only. 
Whenever they were found in preparations which originally contained 
monoclinic crystals it could be assessed from the lozenge angles that the 
crystals have become orthorhombic by the time the striations appeared. 
Whenever striations of the kind as in fig. 3 were seen the crystal (or part, 
of the crystal) had become hexagonal (interfacial angle ~ 60°) 
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All these effects were in first approximation irreversible, shown by the 
fact that they could be observed after the specimens had cooled to room 
temperature. However, the a lines which formed at the lowest 
temperatures were sometimes observed to recede slowly over a period 
of days. 

Optical observations of the effect of heating in situ with the hot stage 
revealed that the ridges appeared suddenly as short lines which then 
grew lengthwise across the crystal. In the course of the same experiments 
changes of the lozenge shapes could also be directly observed. Those 
corresponding to the monoclinic habit (interfacial angle 75°) were 
sometimes seem to change into the orthorhombic (angle 67°) and 
later into the hexagonal macroscopic shape. A transformation front 
was seen to sweep over the whole crystal. All these changes were 
irreversible, no reversion to the original form being seen on cooling. 

A solid substrate was not essential for the striations to form ; they 
also appeared when the crystals were floating on the surface of water 
which was heated to the desired temperature. 

The crystals which showed the striations best, were usually too thick 
for direct transmission electron microscopy, therefore they were examined 
by replication. Figures 4-9 show some examples of these investigations. 
As seen, the striations detected first by the light microscope are really 
regular ridges or roofs. Figures 4 and 5 show examples of a lines which 
were usually the largest and most regular. They are mostly roof-shaped. 
Figures 6 and 7 contain b lines in addition to some a lines which have 
still not fully receded. The b lines formed at these lower temperatures, 
where they may exist together with a lines, are flat topped. When formed 
at higher temperatures also the b ridges are roof-shaped but on a smaller 
scale than the a roofs. The b roofs are usually close together giving a 
corrugated appearance to the whole crystal. This state has already set in 
on the top right of the crystal in fig. 7. This erystal was evidently in the 
process of transformation. In crystals, as in figs. 6 and 7, the b ridges 
continue along the same straight line even when intersected by a lines. 
This is rather remarkable in view of the fact that they developed when 
the a ridges were already present. Figure 8 shows the centre of a crystal 
containing both b ridges and ridges parallel to the lozenge faces ({110)). 
Lines along (110) are also visible in fig. 9, where otherwise regular roof- 
shaped ridges perpendicular to the crystal face ({130) in the case of a 
hexagonal subcell in terms of orthorhombic indices) represent the 
dominating feature. Here, as seen from the shadowing, some of the 
ridges are reversed, i.e. they correspond to grooves and not to protrusions. 
This makes it plausible that the ridges even when seen as protrusions are 
in fact hollow in their interior. 

Even if crystals, as in figs. 4-9, were too thick for good transmission 
electron micrographs it was possible to obtain some kind of image in 
transmission before the crystallinity was destroyed by the beam. The 
images were necessarily very weak on the negatives and were conse- 
quently of poor pictorial quality. Figure 10 shows a lines in such a 
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transmission picture. It can be seen that some lines reveal diffraction 
contrast as one side of such ridges could be darker than the other in 
unshadowed specimens. In the case of really thin crystals suited for 
transmission work, the regularity of the striations was not so obvious. 
Such crystals could crack along the ridges, as shown by fig. 11, which 
was taken with a heating stage within the electron microscope. In spite 
of the difficulties encountered in such hot stage work, referred to earlier, it 
was possible to observe the formation of these cracks in situ. In other 
cases of such thin crystals, the striation showed up as pleats, particularly 
those along the b direction (fig. 12), which in the earlier stages became 
apparent through Bragg extinction fringes (fig. 13). In some cases it 
could be established from the corresponding diffraction patterns that 
in the course of this pleating the lattice is tilted around the short lozenge 
diagonal (b axis). 


§ 3. Discussion 


The main part of this discussion will be confined to the regular ridges. 
It will be attempted to relate these ridges to the arrangements of the 
paraffin chains in the crystals. In view of the ever-increasing complexity 
of the literature on paraffin crystals it may not be superfluous to lay down 
a few basic principles first. 

In the crystals of paraffins the molecular chains are in a parallel packing. 
The question of the details of this packing can be resolved into three 
parts: (1) the lateral packing of the CH, groups, that is the main body 
of the chains, (2) the packing of CH, groups, that is the relative position 
of the chain ends, in adjacent chains, (3) the end-to-end arrangement of 
consecutive chains. As the normal habit is of the layer type, with the 
layers consisting of parallel stacks of paraffin molecules, the above three 
points correspond to the structure within the interior of an individual 
layer, to that at the surface of an individual layer and to the mode of 
packing of subsequent layers respectively. 

(1) The paraffin chain consists of a repetition of CH, groups which 
defines a small repeating unit, the subcell with edges as, bs, cs within 
the true unit cell a, b, c, the latter comprising the full length of the chain 
or possibly several chains (Vand 1951, v. Sydow 1956). One of the 
subcell edges corresponds to the repeating CH, groups (two CH, per 
repeat) along the chain, the other two edges being determined by the 
relative positions of adjacent chains. Accordingly, the subcell structure 
defines the lateral packing of the chains. Two kinds of subcell are known 
to exist in paraffins. One, the more general orthorhombic, contains two 
molecular repeats, the other of triclinic symmetry contains only one 
(see e.g. v. Sydow 1956). Only the orthorhombic subcell is relevant to 
the crystals to be discussed here. For the case of an infinite paraffin 
chain (polyethylene, where the subcell is usually regarded as the true 
unit cell) it is shown by fig. 14. At high temperatures the orthorhombic 
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subcell can acquire higher symmetry and become hexagonal in which 
case a= 3/4/36 in fig. 14. 


Fig. 14 


(b) 


Orthorhombic subcell structure of a paraffin of infinite chain length (poly- 
ethylene, after Bunn), (a) perspective view; (b) projection along the c 
axis (molecular chain). 


(2) If the adjacent chains are displaced lengthwise with respect to 
each other by an amount corresponding to the subcell edge along the 
chain direction the subcell structure is maintained unaltered except that 
the terminal CH, groups will be at altered levels. If they are all in one 
level the planes defined by them will be perpendicular to the chains 
(fig. 15a), otherwise they will be oblique (fig. 156). In the first case 
the structure may be orthorhombic, in the latter always monoclinic 
or triclinic (see below). Clearly a large number of obliquities can be 
realized by shearing the structures along the c¢ direction on different 
shear planes through different amounts. The basic requirement is that 
the terminal planes must be indexable in terms of the subceli. Thus 
in fig. 15 a, the terminal plane is (001)s, while in fig. 15 6 it is (h, ko Ws 
where at least one of # and k is not zero. The suffix s is used here and 
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subsequently to denote subcell indices. While this scheme for generating 
paraffin structures was implied in some of the early works on paraffinoid 
substances, it was Miiller (1927) who first stated explicitly that the 
cross-sectional area of the unit cell projected on a plane perpendicular 
to ¢ was identical for different monoclinic angles. The relation between 
different oblique structures having the same chain packing was 
systematized by Schoon (1938). The subcell concept was first introduced 
by Miiller (1928), without actually naming it as such, in order to facilitate 


Fig. 15 


(b) 


Simplified two-dimensional representation of the arrangement of the paraffin 
ace (4) a a vertical stucture; (6) in an oblique structure. Solid 
circles are , open circles are CH, groups. Th i 
indicate the subcell. a ona: Ree es 
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structure determinations, and was utilized later under this name by 
Vand (1951), followed by v. Sydow (1956) to simplify the classification 
of paraffinoid structures in the way presented in this paper. 

(3) The fitting of consecutive chains end to end is of lesser importance 
for the present discussion. It will be added here for completeness. Here 
we have the possibility of distinguishing between the structure of odd 
and even paraffins. Figures 16 (a and b) shows the same kind of end-to-end 
relation for odd and even paraffins. In both cases the planes containing 
the terminal groups are perpendicular to the chains, nevertheless (a) is 
orthorhombic while (6) is monoclinic. 


cena 


Fig. 16 


eae < “ 


(b) 


Simplified two-dimensional representation of the unit cell in a vertical structure 
for the case of (a) odd; (b) even paraffins. Solid circles are CH,, open 
circles are CH, groups. 


When determining the symmetry different criteria are used by different 
investigators. It is customary to distinguish between the orthorhombic 
and monoclinic forms on the basis of the ¢ repeat directly derivable from 
the 001 reflections. While this would give the correct answer for cases 
like fig. 15 (a and 5), fig. 16 b would be classed as orthorhombic on this 
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basis because with sufficiently long chains the shortening of the ¢ period 
may not be normally noticeable. 

As structures in fig. 16 (a and 6) are clearly distinct from that in fig. 15 b 
we propose the distinction of vertical structures (which would also 
include fig. 166) and oblique structures (fig. 15) in addition to the 
strict classification by lattice symmetry. In the vertical structures, 
the planes containing the CH, groups should be indexable as 001; (in 
terms of the subcell), while in the oblique structures hs and ks would 
not be both identically zero. In the literature hexagonal phases are also 
referred to. Evidence for this is usually obtained from the hk0 reflections 
which show that a=3/,/3b. However, this only gives the subcell 
structure. The scheme outlined above (fig. 15) could still be operative 
also with a hexagonal subcell. Thus many oblique structures could 
correspond to a hexagonal subcell. Whether the structure is truly 
hexagonal or not can only be determined from /&/ reflections. 


Fig. 17 


UM 


(b) 


Schematic representation of a crystal layer having (a) a vertical structure; (b) 
an oblique structure. 


The lozenge-shaped paraffin crystal layers consist of parallel stacks of 
molecules. In the case of the vertical structure the molecules are 
perpendicular to the crystal planes (fig. 17a), in the oblique case 
(fig. 176) they are not. Since molecular packing requires shears by 
discrete amounts, and intermediate amounts of shear involve higher 
free energies, transition should not occur simultaneously over the whole 
crystal but spread out from wherever it was first initiated. It will now 
be proposed that such a progressive transformation could occur by the 
molecules slipping past each other while maintaining their original 
orientation with respect to the geometrical outlines of the crystal. The 
result will be a sloping basal plane which could lead to a ridge as in 
fig. 18a. It is propsed further that the ridges observed in the crystal 
are of this origin. Accordingly, a striations would arise through a simple 
shear on (010) planes, b striations by simple shear on (100) planes and 
other striation directions would have their corresponding shear planes. 
Consequently in single crystals like the ones in figs. 6—9 several oblique 
structures differing by the amount of shear and by the plane of shear 
are represented. 
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It will be assumed that the orthorhombic or, if at high temperatures, 
hexagonal subcell is preserved during this transformation. Subcell 
theory requires that the planes containing the end groups, hence the 
planes of the ridges, should be defineable in terms of rational subcell 
indices. We shall now proceed to index some of these planes from the 
electron micrographs. The a ridges were the most suitable for this 
purpose, being the largest and most regular, being clearly separated by 
flat areas of the unaltered crystal. These ridges are terminated by 
oblique faces or by a sequence of such faces (fig. 5). These terminating 
faces must have originated through a shear on (100) planes, hence must 
have the same origin as the b ridges, consequently they should also be 
indexable with reference to the subcell coordinates. 


Fig. 18 


(a) 


(b) 


Schematic representation of a crystal layer with a localized transformation. 
(a) from a vertical to an oblique structure; (>) from an oblique structure 
to one with a larger obliquity. 


While without further information the indices themselves cannot be 
assigned, the ratio of the indices of planes meeting along a line can. 
The subcell indices of the sides of the ridges are (Okl)s, those of the 
terminal planes are (h0l)s. The line along which they meet is the axis 
of the zone to which both belong. By measuring the angle enclosed 
by this edge in projection and a given crystal axis the component of this 
zone axis in the (001) plane, hence the U and V zone coordinates and 
thus h/k can be determined. 

The angles between these interplane edges and the a axis (the direction 
of the ridges were measured in several cases. Most of the figures obtained 
fell in two main groups. Those which were the most regular, hence most 
accurately measurable, were between 50 and 56°. A number of measure- 
ments gave figures 17-23°. However, other values were also obtained. 
Such were 24-28°, and much smaller angles like 10-11°, and <5° 
possibly 3°. The angle 53° corresponds to h/k=2, which clearly defines 
the first group. The angle of 20° gives h/k=1/2 which defines the second 
group. 24° gives h/k=2/3, 28° 3/4, 10-11° 1/3 or 1/4 and an angle less 
than 5° would require at least h/k=1/10 or possibly a still smaller ratio. 
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For obtaining the indices themselves the heights of the ridges are 
required. These were obtained by shadowing the specimens at an angle 
of 10° and by measuring length of the shadows cast by suitably oriented 
ridges. From this and from the ridge width the slope of the planes could 
be directly assessed. 

It was found that the inclinations of the side faces of the a ridges 
were mainly of two kinds. One corresponded to the range of 25-29° 
while the other to that of 15-18°. Intermediate values such as 19°, 
22° were also measured ; in fact the slope could vary along the same 
ridge and could also be somewhat different on both sides of the same 
ridge, although in general, the ridges were symmetrical. 

The slope of 27° corresponds to a shear on (010) by 2-5A per unit cell, 
i.e. by cs/as. Accordingly the basal planes would be (011)s. The slope 
of 15° requires a shear of $c, per unit cell, the corresponding subcell 
index being (012),. We conclude that these indices define the planes 
of the corresponding ridges in the crystals. Intermediate values would 
correspond to indices like, e.g. (023). Knowing the ratio h/k from the 
interplane edge directions the terminal faces could also be similarly 
indexed once the side planes of the ridges were known. In suitable 
photographs, indices of the terminal planes could also be determined 
directly from the ridge height. The indices of the terminal faces found 
by the above methods were (201); ; (101)s; (102),; (104); or (103), ; 
(203), or (304)s, ~(1, 0,10). In the case where alternatives are given 
both types of planes are permissible within the accuracy of our 
measurements. 

One may now consider the structures giving basal planes as determined 
from the electron micrographs. With (010) or (100) as shear planes, 
and with a homogeneous shear, preserving the same subcell structure, 
the displacement of successive planes of molecules must either be a 
multiple of cs, or an odd multiple of $cs. accompanied by a rotation 
through 180° about cs of the molecules in every second plane. The 
largest shear observed in fact corresponds to the latter case, with 
displacements of }cg between successive planes. 

We must infer either that the rotations through 180° of half the 
molecules occur, or that the shears are not homogeneous (rejecting as 
improbable the change of subcell structure from Pnam to Pbnm known 
as a paraffinoid subcell structure only in 2D, methyloctadecanoic acid, 
Abrahamson 1959). The smaller observed shears cannot be homogeneous, 
but appear to be regular, with displacements between every second, 
or every third plane, etc., since the inhomogeneity is not resolved in 
the electron microscope, and a principle of rational indices does seem 
to apply. 

In general, let (hkl); be the basal plane and hg be expressed as nH 
and ks as nA, where n is the largest common factor of hg and ks. Then 
the molecules within every /sth hkO plane will need displacing by n . cs/2 
along the c-axis. Clearly a very large number of oblique structures can 
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be realized in this way for any particular shear plane. Such possibilities 
were not included in earlier schemes for generating paraffin structures. 

High index planes implying long-range order would not be expected. 
Nevertheless, such are occasionally seen to exist amongst the terminal 
planes of the ridges where slopes so low that they could not be accurately 
measured any longer on our photographs could be observed corresponding 
to /s indices at least as high as 10. The terminal parts are really 
accommodation structures forced into being by the main body of the 
ridges, consequently they can be far removed from equilibrium structures. 
The existence of very high index planes is therefore not surprising, 
neither is the occasionally observed curved boundaries of the terminal 
faces indicative of a continuous variation of the slopes. Nevertheless, 
there is a tendency towards equilibrium structures even within the 
terminal parts, as indicated by the frequent occurrence of the low index 
faces like (101); and (102). 

The foregoing implies that the molecular direction remains unaffected 
in all these ridges, consequently the same hk0 reflections should be 
either in or out of reflecting position within the different structures in 
the same crystal. This follows from the fact that in first approximation 
the hk0 reflections are determined by the subcell which should be identical 
in all the structures, and should also retain the same orientation in these 
transformations. As with paraffins, the diffraction contrast in electron 
. micrographs is predominantly due to hAk0 reflections one would not 
expect such a contrast in the image of the ridges. However, near 
diffracting position only a very small tilt, let us say 1°, towards or away 
from diffracting position would already be sufficient to produce a strong 
contrast effect. Such small deformations could occur if only for reasons 
of mechanical rigidity in the ridges. Also it has to be remembered that 
though the molecular displacements in a transformation of crystal 
structure are well described, within a layer, as simple shear displacements 
along the axes of the chains, the resulting changes of structure at the 
surfaces of the layers are likely to cause small readjustments in the 
stacking of one layer on another, the consequence of which will be small 
rotations of the subcell structure, which can well suffice to account for 
the observed diffraction contrast between the two sides of a roof. This 
might be the case in some of the ridges in fig. 10 where one side of some 
of the roofs is darker than the other. That the difference in orientation 
can only be slight, and consequently does not invalidate our interpre- 
tation of the roof structures is shown by the fact that in some ridges the 
dark zone is in the centre of the roof and in others it is absent altogether. 
Also the dark sides can be opposite ones in the different ridges. This 
would not happen if the different parts of the ridges were associated with 
distinct molecular orientations which differ from that in the untransformed 
crystal and also differ for the two sides of the ridges. 

In none of these cases did we consider the change in the subcell structure 
due to the change of temperature. It was mentioned that the ridges 
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were unambiguously observed in the orthorhombic (vertical) crystals. 
In principle, it is not excluded that crystals with an oblique structure 
as in fig. 176 should form ridges. Such ridges would have to be 
asymmetrical (fig. 18 5). 

The regular ridges were most readily observed in crystals about 
twenty—forty layers thick. In thinner crystals, more suited for trans- 
mission electron microscopy, the striations lacked the geometrical 
regularity. It is to be expected that thin layers will not be able to 
form large regular ridges, if merely for lack of rigidity. As the shear 
spreads the crystals could eventually split along the shear planes. 
Apparently this has occurred in the crystal of fig. 11. The thinness of 
the crystals would also lead more readily to bending or pleating of the 
layers (fig. 12). which would also cause noticeable disorientation to show 
up in the diffraction patterns, also revealed by the corresponding Bragg 
fringes in fig. 13. Nevertheless, there appears something systematic 
in the pleating along the b direction which is noticeable also in thicker 
crystals. This can correspond to a tilt of the unit cell around b, as 
indicated by some of our experiments. The large expansion along a 
might be partly responsible for this effect but there is evidence that this 
could not account for it fully. Transitions from vertical to oblique 
structures, such as in fig. 176, involving a tilt of the molecule cannot 
be excluded. This would increase the surface area and if it occurred in 
parts or progressively it could push up the material in the form of pleats. - 
While such less regular morphological features have been mentioned here, 
a more detailed discussion or description would be beyond the scope of 
the present paper. 

Besides the purely structural implications, the present observations 
raise a number of further questions. A few of them will just be 
mentioned without attempting an explanation at this stage. The 
formation of ridges implies a change from the vertical to an oblique 
structure, that is from a system of higher to lower symmetry on raising 
the temperature. This is contrary to the usual trend in crystal trans- 
formations. Similarly puzzling is the behaviour that such structures 
should become vertical again and change into ones of a different kind of 
obliquity. Further, the structural changes in question are essentially 
irreversible. Crystallographic changes in paraffins below the melting 
point are usually referred to as reversible in the literature, both as regards 
the changes confined to the true unit cell (monoclinic-orthorhombic, or 
rather oblique and vertical in the present nomenclature) and as regards 
the changes also occurring in the subcell (orthorhombic—hexagonal). 
Clearly the thermodynamics of these systems is still largely obscure. 
Finally, the conditions under which the present effects appear are very 
sensitive to a number of variables. Besides the generally recognized 
criterion of purity (e.g. Shearer and Vand 1956, Teare 1959), the size 
and thickness of the crystals, the rate of heating, substrate, etc., are 
apparently important factors. 
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It will be seen from the next paper of the series that the effects described 
and discussed here have their counterparts in polymer crystals, which 
without this background with paraffins could hardly have been 
interpreted or even followed up effectively experimentally. 
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ABSTRACT 


The morphology underlying the characteristic striations observed in poly- 
mer crystals has been examined in more detail than hitherto. The striations 
were found to correspond to corrugations in the crystal layers. Both pleat 
and roof-ridge corrugations were observed in polyethylene crystals. These 
latter are pictorially similar to those in paraffin crystals discussed in Part I. 
The analogy between the two cases is pursued with the folds replacing the end 
groups in the polyethylene case. A possible common origin for both types 
of corrugation is proposed to be the collapse of non-planar crystals. A 
particular non-flat-based pyramidal model is proposed. Evidence for the 
existence of such crystals is presented. A variety of observations, 
including the non-planar crystals, is shown to be explainable in terms 
of a particular type of packing of the folded molecular chains. Further, 
it is suggested that if folding also occurs in the bulk polymer, then roof ridges 
of the type discussed might account for its well-known four-point low-angle 
X-ray pattern. 


§ 1. [INTRODUCTION 


THE first observations on polymer single crystals showed that they 
consisted of flat plates growing perpendicular to themselves by screw 
dislocation mechanisms (e.g. Keller 1957). Polyethylene single crystals, 
in particular, have similar habits to those of the orthorhombic n-paraffins. 
Both the electron microscopic images and the diffraction patterns were at 
first sight indistinguishable for crystals of the two materials. Thus in 
polyethylene, also, the molecular chains were perpendicular (or nearly 
so) to the crystal layers. As these polyethylene layers are only 80 to 
160A thick and the molecules microns long, it was necessary to postulate 
that a regularly and sharply folded molecular configuration existed 
within them (Keller 1957). More detailed examination of polymer 
single crystals has revealed a number of new features, in particular 
characteristic striations (e.g. Agar et al. 1959). As described in Part I, 
striations were, in fact, also found in paraffins, of which at least the most 
regular could be given a satisfactory interpretation. The present paper 
considers some examples of striations in polymer crystals and draws 


parallels with paraffins as far as is possible. 
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§ 2, EXPERIMENTAL 
2.1. Materials and Specimens 


The experimental material was mainly the linear ‘polyethylene, 
Marlex 50, but reference will also be made to polyoxymethylene (CH,0),. 
In all cases, the polymers were dissolved in hot solvents and crystallized 
by cooling these solutions. In this way a suspension of single crystals 
formed which could be examined after drying down on a microscope 
slide. By using carbon-coated slides, and subsequently floating the film . 
off on water the crystals could be examined electron microscopically. 
The illustrations shown here were selected from a body of material which 
has accumulated whilst following several different lines of research. 
Because of this, to detail the various preparative methods used would 
take too long here and we shall merely quote the types of specimens 
used. For polyethylene, these were four. 


(1) Crystals which were truly lozenge-shaped, obtained by various 
methods of crystallization at temperatures below 80°c. 

(2) Truncated lozenge crystals, formed at slightly higher temperatures 
(the relative size of the truncating faces {100} increasing, for a given 
solvent, with the temperature of crystallization). Those mostly referred 
to were formed at 90°c and were filtered off at that temperature before 
cooling. 

(3) Crystals heated above their original crystallization temperature 


while still in suspension and subsequently subjected to various conditions 
of cooling. 


(4) Crystals heated dry on a microscope slide and examined after 
re-cooling. 


The polyoxymethylene crystals referred to were formed as a solution 
cooled to room temperature. 


2.2. Some New Observations 


As first observed, the striations in polyethylene crystals of type (1) 
were irregular, wavy lines, approximately perpendicular to the faces of 
adjacent crystal quadrants (Keiler and O’Connor 1958, Agar et al. 1959). 
These were observed as Bragg fringes in suitably taken electron 
micrographs, but were later found to correspond to an actual morphology 
(see, e.g. Keller and Bassett 1960). These physical corrugations were 
for the most part, transient and largely disappeared soon after deyine 
the crystals down. When the regular roof-shaped ridges were recognized 
in paraffins (Part I), the possibility of drawing a parallel between these 
and the polymer striations presented itself. For this, a clear pictorial 
evidence was required. The experimentation which followed aimed at 
obtaining such evidence by improving the methods of specimen 
preparation and observation. In the course of this work, the regularity 
of the striations became increasingly apparent. In the first instance, 
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the regularity was observed in the Bragg fringes described above, which 
could take the forms of straight lines or bands (figs. 1, 2+). These were 
seen best in specimens examined immediately after drying -down. 
Moreover, the regularity and definition of these fringes deteriorated 
continuously on exposure to the electron beam. The underlying 
morphology, however, could be preserved at least partially by shadowing 
the specimens more heavily than is usual, immediately on drying down. 
This makes the structures more permanent and also more resistant 
to electron damage. Using this method photographs at higher magnifica- 
tion could be obtained. Figure 3 shows a central part of a crystal obtained 
in this way. As seen there is a line structure corresponding to the 
direction of the extinction fringes in figs. 1 and 2. The most pronounced 
lines appear to be pleats, but others parallel to them, regular ridges 
as best recognizable from fig. 4 which is an enlarged detail of such a 
crystal. As seen from figs. 3 and 4 such ridges are crossed obliquely 
by another striation. This second striation direction is best visible 
where it passes through the first-mentioned ridges but can also appear 
alone in other parts of the crystal. In some other photographs this 
latter feature could be seen more prominently. Figure 5 is an example 
of a multilayer crystal, with a few striations on the top layers only. 
It has been a general experience that the striations on the higher layers 
were progressively more permanent. The decorating effect of small 
latex-type particles of unknown origin along the crystal edges may be 
noteworthy. 

The regular roof-shaped nature of the striations is also evident in the 
truncated lozenges of type (2). Here they are usually most pronounced 
in the sector bounded by the truncating face where they are parallel to 
the short diagonal of the crystal (fig. 6). The detail in fig. 7 reveals the 
geometrical regularity of these ridges. (The overgrowth corresponds to 
crystals formed at a lower temperature (see Keller and Bassett 1960, 
fig. 31); the study of such effects is part of another programme of work 
and will not be commented on further here.) There are usually signs of 
the ridges continuing beyond the sector boundary changing their direction 
to that of the kind found in the true lozenges described above. This 
effect became more pronounced when the suspension containing such 
erystals was heated above 90°c in which case these ridges could be 
observed to extend further into the adjacent sectors until eventually 
they pervaded them entirely (Keller and Bassett 1960, figs. 28, 29). In 
these latter cases the direction of the lines in the adjacent sectors could 
change somewhat so as to become more nearly parallel to the short crystal 
diagonal. 

The transient nature of the striations has already been mentioned. 
These previous comments referred to the disappearance of the lines after 
drying down. The suspensions themselves did not alter. On drying 
down the same kind of crystals were obtained even after the suspension 
had been kept for months. However, those striations which appeared in 
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the truncated lozenges on heating the suspension could only be observed 
if the crystals were deposited within a very few hours of the end of the 
heat treatment. 

Apart from such rather special effects, in general the effect of heating 
was to enhance the striated appearance of the crystals. Here again 
the striations could be either pleats or regular roof-shaped ridges or 
stages between these two, even within the same crystal. 

As reported elsewhere (Bassett and Keller 1959, Keller and Bassett 
1960) striations could also appear on heating the dry crystals on a micro- 
scope slide. On the first layer, in contact with the substrate, the lines 
were either perpendicular to the lozenge faces or parallel to the short 
diagonal; on the higher layers they were always of the latter type. While 
these usually manifested themselves as pleats, in the early stages they 
showed a definite geometrical regularity. In all the heating experiments, 
in suspension or dry, the striations formed first on the top layer of the 
crystals. 

Observations of pleats are not new. These have been reported in 
polymer crystals under a variety of circumstances. In some truly lozenge- 
shaped crystals they outline a small lozenge at the centre similar in shape 
to the main crystal, but rotated through 90° with respect to it. Im most 
other cases these pleats lie along the short crystal diagonal, though 
occasionally they are found also along other directions. These latter 
pleats are usually triangular in shape and amongst our preparations could 
be most clearly seen in the truncated lozenges (fig. 8). 

Striations and a corresponding morphology do not appear only in 
polyethylene, but may well be a general phenomenon in long chain 
crystals. Polyoxymethylene, for example, forms large, rather irregular, 
crumpled hexagonal crystals often containing radial pleats or thickenings. 
Electronmicrographs revealed regular striations which were shown up 
best by diffraction contrast (fig. 10) but there is evidence for an under- 
lying morphology (fig. 11). 

For the sake of continuity the description of a number of further 
observations will be included in the appropriate part of the Discussion. 


§ 3. Discussion 


At first the basic structural data will be briefly recapitulated. The 
subcell (often called the unit cell) of polyethylene is orthorhombic (Bunn 
1939, see fig. 1, Part I). In the crystals the long lozenge diagonal is along 
the a, the short along the b axis. The prism faces always present are 
{110}. The truncating prism faces present in crystals grown at higher 
temperatures are {100}. The ¢ direction, that of the chains, is approxi- 
mately perpendicular to the lamellae. The molecules are expected to 
fold along the boundary prism faces of the crystals. This implies that the 
crystal consists of as many distinct sectors as there are prism faces, the 


sectors being defined by the direction of the folds within them (Keller 
and O’Connor 1958, Bassett et al. 1959). 
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The striations provide a clear example of sectorization. Thus, they 
change in direction at the sector boundaries (figs. 1, 2, 6) and the different 
sectors are distinguished by the different directions of the striations 
within them. The existence of distinct sectors is confirmed by diffraction 
effects since different diffraction conditions may be satisfied in the different 
sectors (Agar et al. 1959, Bassett et al. 1959). This‘ selective diffraction 
allows different sectors to be brought into Bragg extinction positions in 
different orientations to the electron beam. The crystal in fig. 8 has two 
{110} sectors in Bragg contrast; we have other micrographs with the {100} 
sectors alone diffracting—see e.g. fig. 2 in Bassett et al. (1959). Figure 1, 
with opposite pairs of quadrants in different contrast, illustrates the same 
effect. | 

Until recently striations were only mentioned as far as they served to 
confirm the hypothesis of sectorization. The nature and origin of these 
lines was left unexplained. The transient and rather irregular line struc- 
ture did not lend itself to further interpretation. The present results 
(figs. 4, 7) have revealed an unmistakable similarity between these lines 
and the paraffin ridges reported earlier. In paraffins these ridges were 
interpreted as resulting from new phases, the different parts of the ridges 
containing structures with different obliquities, the sloping planes of the 
ridges being the planes containing the end groups. In the case of polymers 
consisting of folded chains there are no end groups, but the folds should 
play a similar part. They define a basal plane which should be indexable 
in terms of the subcell (Part I), and accordingly could give rise to vertical 
and oblique structures (figs. 2, 3 in Part I). By geometrical arguments 
alone all the possibilities discussed in connection with the paraffins should 
be realizable. 

Having obtained photographs with more regular striations than before 
the directions could be more accurately measured. They were found, in 
contradiction to earlier reports, to be always crystallographic. In the 
true lozenges they were either (130)s, (250), or (120), where the suffix 
s refers to the subcell. The common feature of all these is that they are 
nearly perpendicular to the {110} faces bounding the adjacent sectors. 
The second kind of striation sometimes seen to cross the ones mentioned 
above (figs. 3 and 4) is (230)s. This is distinguished by the fact that 
[230], and [230], are almost at right angles to each other, the corresponding 
lines in the different sectors defining a square. In the {100} sectors of the 
truncated lozenges the lines were along b, [010]. While the above directions 
were the principal ones they are by no means unique as less frequently 
others were also observed. Such were the ones parallel to a, [100], and in 
one instance along the lozenge faces, i.e. (110). 

The lines along (130)s, (250), (120) define the sectors by changing 
directions at the boundaries. The b lines in the truncated lozenges 
define the {100} sectors by being confined to them. (The sectors are 
often also delineated by elevations along the boundaries (fig. 6)). Some- 
times, however, the b just as the rare a lines may pass through the whole 
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crystal. A possible reason for the particular choice of shear planes in the 
different sectors will be given later. In the {100} sectors of truncated 
lozenges they are parallel to the assumed fold plane; accordingly layers 
of folded molecules would only need to stack differently without any need 
to change the arrangement of the folds within the layers. In other cases, 
however, the shear plane would have to pass through folded layers which 
would affect the folding of the molecules themselves. In the case of 
paraffins the ridges formed mainly on heating. There is a parallel in 
polymers where heating enhanced the striated effect. 

Regular striations are not restricted to polyethylene, as shown by our 
observations on polyoxymethylene (figs. 10, 11). It is possible, therefore, 
that the same explanation as above would also be valid for other polymers, 
polyoxymethylene in particular. 

One may enquire about the origin of the ridges formed in untreated 
crystals. One possibility is that in analogy with paraffins they represent 
phases stabler at the higher temperatures where the crystal originally 
formed. When the crystals cool to room temperature, while still in 
suspension, parts of them might change into the modification stabler at 
the lower temperature. This issue could only be decided if the crystals 
could be examined at the temperature at which they were grown. 

The ridges and pleats would also originate from a dished pyramidal 
habit of the crystal layers. That polymer crystals may indeed be non- 
planar was inferred from the triangular pleats seen in crystals as in fig. 8 
(Bassett et al. 1959, Niegisch 1959). ‘These pleats, it was suggested, were 
formed by the regular flattening of a hollow pyramidal crystal as it dried 
down on a substrate, the pleat accommodating the excess material. 
Alternatively, the crystal could crack, which has also been observed (both 
pleats and cracks are seen simultaneously in fig. 5). Other kinds of pleats, 
e.g. the square, round, or fourfold star-shaped ones in poly-4-methyl- 
pentene-1 (Frank et al. 1959) and the six-fold stars in polyoxymethylene 
reported here are equally compatible with the idea of a collapsed pyramid. 
There were, in fact, some theoretical reasons which gave support to this 
suggestion. Accordingly the dished pyramidal habit could be a con- 
ceivable consequence of the sectorization of the crystals (Bassett et al. 
1959). Consider a crystal consisting of four {110} sectors (true lozenges). 
The existence of folds along {110} planes would lower the symmetry of 
the lattice in a sector because the two 110) directions would not be 
equivalent any longer as in the truly orthorhombic cell. In case the 
diagonal containing the fold became shorter this could make the central 
angle by which the quadrants meet less than 90°. Consequently the 
four quadrants would not join while staying in a plane and a ‘ dished ’ 
pyramid would result. The hollow pyramidal habit could also be arrived 
at by considering the possible modes of packing of the layers containing 
the folded molecules. The folds need not define a basal plane normal to 
the chains but could form a staggered configuration. This would lead to 
oblique basal planes with different directions of obliquity in the different 
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quadrants and thus to hollow pyramids. This explanation suggested to 
us by Reneker and Geil (private communication) is, in fact, an extension 
of the ideas underlying the formation of ridges. 

However, such theoretical considerations, while very attractive cannot 
be considered as compelling unless backed by unambiguous experimental 
evidence. The triangular central fold was certainly a strong indication, 
but was not in itself decisive evidence as it was conceivable to account for 
it otherwise. Clearly direct observation of the pyramid in its uncollapsed 
state was required. This was achieved by the following experiment. To 
eliminate the effect of a possible flattening a drop of suspension containing 
the same kind of crystals as in fig. 8 was placed on the surface of a viscous 
solution of polyvinyl alcohol. After evaporation of the xylene, crystals 
could be seen floating on the polyvinyl alcohol solution with the optical 
microscope. It was noticed that the different sectors could appear bright 
or dark even as seen optically (see fig. 20, Keller and Bassett 1960). This 
effect was preserved after the polyvinyl alcohol solution had solidified. 
At this stage the surface was metal shadowed, coated with a thick layer of 
carbon, and the polyvinyl alcohol dissolved in water. It was assumed 
that the heavily carbon-coated crystals examined electronmicroscopically 
would retain the configuration they had when deposited on the liquid 
surface originally. In fig. 9 one of the first photographs taken of this 
subject is seen. The central fold is absent and the different sectors are in 
different contrast near the centre as required for a pyramid. The obser- 
vation that the different sectors could be brighter or darker as seen with 
the light microscope would follow from the fact that such sectors together 
with the underlying substrate material would be prisms which deviate 
light so as to increase or decrease the corresponding image intensity in the 
microscope. 

Hollow pyramids could flatten in two distinct ways. In one case, if all 
the molecules are parallel to the pyramidal axis, they could slide past each 
other as visualized in the ridge formation. If the sliding past of mole- 
cules is not uniform over the whole crystal and the flattening is not 
complete, elevated regions, ridges in particular, could form in the course 
of collapsing. Accordingly, a ridged structure could be an intermediate 
one between a pyramidal and a flat crystal. In fact it has been observed 
that ridges disappear rapidly on standing after the drying down of the 
crystals which accordingly could represent the final stage in the collapse 
of the pyramid. No definite change of area on the substrate was observed 
on the disappearance of the ridges in agreement with these ideas. In 
contrast to the regular striations, the pleats do not disappear in time. 

Alternatively, one could consider the molecules in the crystal layers as 
fixed. In this case the layer would tilt as a whole when the pyramid 
collapses. This, in fact, would lead to excess material within the basal 
area of the collapsed pyramid, which would result in folds as, e.g., in fig. 8 
being pushed up. These folds would be associated with a tilt of the 
molecules which should be noticeable from the diffraction patterns. This, 
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in fact, was confirmed both by x-ray and electron diffraction. X-ray 
patterns were taken of aggregates obtained by sedimentation where the 
crystals sedimented horizontally (Keller and O’Connor 1957). Such 
ageregates formed oriented thin films. As reported previously, the 
molecules were approximately perpendicular to the film planes. The 
x-ray patterns obtained with the beam parallel to the film plane were 
essentially fibre photographs with all {k0} reflections on the equator. 
The present experiments on aggregates of crystals like those in fig. 8 
revealed that the {200} reflection was split azimuthally by some 20°. The 
{020} was not split while in the {110} the presence of a smaller split was 
poRsible but could not be ABNHEaINEA with certainty owing to the large 
inherent disorientations of the arcs. For similar reasons, the inter- 
pretation of the patterns is necessarily imprecise. It is clear, however, 
that the b,* axis lies essentially in the specimen plane while a,* is rotated 
out of it. (The asterisk denotes a quantity in reciprocal space.) This 
implies that the molecules are now inclined to the fibre axis. In terms of 
the single crystals forming our specimens this would involve crystal 
layers tilted to the specimen plane for a vertical structure (fig. 15 a, 
Part I). For an oblique structure the molecules would be tilted within 
the layers which could themselves be lying flat. Physically, we expect 
the hollow pyramids to collapse during the formation of the x-ray 
specimen. This and other reasons which will shortly become obvious, 
lead us to prefer the latter alternative. 

The main features of the electron diffraction patterns from crystals as 
in fig. 8 were consistent with the x-ray evidence. It was reported pre- 
viously (Bassett et al. 1959) that diffraction was selective, different sectors 
adding different reflections to the pattern given by the crystal as a whole. 
The details are as follows: with the beam normal to the crystal, the {110} 
reflections were by far the most prominent. They originated primarily 
from the {110} sectors, each sector producing only that pair of {110} 
reflections, from planes parallel to the faces of adjacent sectors. The 
(020) reflections were mostly present in strengths, which varied from 
crystal to crystal. The {200} reflections, though, were either absent or 
very weak, in the latter case split azimuthally into two or sometimes three 
spots. The narrow {100} sectors were also diffracting selectively as seen 
from the Bragg fringes (Bassett ef al. 1959) but the reflecting conditions 
could not be determined with greater accuracy. 

The observed selection is what one would expect from the collapse of a 
hollow pyramidal crystal containing triclinic {110} sectors with chains 
parallel to the pyramidal axis; (the existence of {100} sectors will be dis- 
regarded for the moment). We shall assume that the pyramid is not flat- 
based (fig. 12). Such a pyramid can be characterized by the direction of 
the lines along which a plane normal to the pyramid axis intersects the 
pyramid. We shall refer to the contiguous pyramid based on these lines as 
the equivalent flat-based pyramid. We shall assume further, that in our 
case these lines are (130); (or some nearby direction such as (250), or 
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(120); (see Keller and Bassett 1960)). It will be demonstrated that a 
number of observations can be accounted for on this basis. 

If a pyramid is to be flattened without plastic deformation, its faces 
must be tilted as unaltered structural entities. For a non-flat-based 
pyramid in which the surfaces containing the folds are {hkl\s, then a 
condition of flattening is that the new orientation of a sector can be 


Fig. 12 


pyramid 
axis 


Diagram of the non-flat-based hollow pyramid. The four basal apices A, A’, 
B, B’ are not coplanar. The lines connecting A and A’, and B and B’, 
are perpendicular to the pyramid axis but are on different levels. Lines 
BD, DB’, B’D’ and D’B define the directions along which planes 
normal to the pyramid axis intersect the pyramid. It is postulated 
that the crystals are pyramids of this kind where these directions are 
<130>,. The base BDB’D’ and apex C define the equivalent flat based 


pyramid. 


related to the old by a single rotation about an axis near to one of the 
<kh0)s directions. In our case, these directions would be of (130); type 
which would account for the observed selection of the {110} reflections 
and with a suitable choice of planes {31/};, for the observed tipping of 
approximately 20° of the {200} reflection out of the specimen plane in the 
x-ray patterns mentioned previously. Planes of the type {312}, and 
{313}, give the best fit with the experimental results here and else- 
where. The presence of the {020} reflection is not explained on this 

model, however, as it would still be tipped at least 10° out of the specimen 
plane, i.e. out of reflecting position (cf. Agar et al. 1959) but not so far as 
is (200). Its presence could be due to misalignments within the crystal, 
formed possibly as irregularities in the collapse of the pyramid or by the 
electron beam. This is consistent with the variable intensity which may 
well be indicative of its accidental origin. In any case, the orientation of 
the {020} planes in a {110} sector is intermediate between that set of { 110} 
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planes which remains essentially untipped and the {200} plus the other 
set of {110} planes. Thus the indecisive behaviour of the {020} reflection 
compared with the clear cut extremes of {110} and {200} is at least broadly 
compatible with the picture proposed. 

The above mode of collapse would push up excess material, i.e. would 
produce pleats, along the corresponding (130); directions. Such lines 
are seen sometimes in the truncated lozenges (e.g. fig. 6) and are systematic 
features in the large non-truncated lozenges (figs. 3, 5). If the postulated 
pyramid collapses through the sliding past of molecules and the formation 
of roof ridges as discussed previously, these ridges would have to be along 
(130); directions, and along different ones in different sectors, which is 
again the observed fact (figs. 3, 5). (For a flat-based pyramid, these 
directions would be {110)s.) Thus the hitherto puzzling direction of the 
striations appears to be accounted for. The choice of pyramid assumed 
can, moreover, be based on considerations of molecular packing. 

We have made three assumptions in our consideration of polyethylene 
lozenges so far. (1) That the molecules fold along the growing faces, i.e. 
{110} planes in the {110} sectors. (2) That the lozenges are pyramids 
with the chain direction parallel to the pyramid axis. (3) That the 
pyramid is not flat-based but has an equivalent flat-based pyramid 
bounded by {310}; prism faces. These three requirements are incom- 
patible with a vertical structure in the folded molecular ribbons lying in 
the {110} growing surfaces. The structure in the ribbon must, in fact, be 
oblique and the adjacent ribbons must pack obliquely, too. For a 
pyramid whose faces are {311}s, {312}s, {313}s ete. adjacent chain segments 
within a ribbon ({110} prism faces) must be displaced vertically by cs, des, 
3¢g and so on respectively. The staggered stacking of adjacent ribbons 
requires the same segment displacements respectively between successive 
segments along (b). The displacements of contiguous chain folds within 
each ribbon are twice these amounts. The first two of these pyramids can 
be achieved uniformly with each fold displaced in the same way relative to 
its neighbours, the third requires non-uniform displacements (see the dis- 
cussion for paraffins in Part I). As already mentioned, the planes {312} 
and {313}; agree best with the experimental evidence. Of the two, {312} 
would give the simpler packing. As, however, the detailed crystallography 
of the folds is not known, {313}s cannot be ruled out @ priori on this basis. 
In these ways it seems we can account for the characteristic striation 
directions in the {1 10} sectors in terms of a particular packing of the folds 
5 those sectors. Similar considerations apply to the b striations in the 
{100} sectors, when considering a corresponding six-sided pyramid. In 
addition, the experimental observation of the striations developed on 
flattening in crystals heated in suspension being nearer to b than are the 
(1380s directions would be accounted for if the pyramids became steeper 
on heating (by a suitable change in fold packing). 

When ridges and pleats are both present, the two modes of collapse 
could be operative simultaneously. The regular ridges would not 
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produce appreciable changes in the diffraction conditions, consequently 
erystals in which the ridge structure is preponderant are not expected to 
reveal appreciable selection in the diffraction pattern. This is largely in 
agreement with electron diffraction results on crystals of the kind in 
figs. 1-5. One would, however, expect selection around the localities 
of the pleats and the fact that selection of the expected type exists in 
zones parallel to the appropriate directions (i.e. ~ (130)s) in the different 
sectors has already been reported (Agar et al. 1959). 

The most definite observations on selective diffraction were made on 
the {110} sectors of truncated lozenges; these were the results referred to 
previously. In these crystals, however, a triangular pleat along b is 
the dominant morphological feature (fig. 8). This is a more complicated 
type of collapse than that previously dealt with and requires a rotation 
around [001]; as well as about a particular (130); direction to relate- 
the sector orientations before and after collapse. The reason for adopting 
this method of collapse rather than that referred to earlier is not obvious. 
In our experience, however, the central fold is most regular in the truncated 
lozenges, though it also appears in some true lozenges, particularly those 
grown at higher temperatures (cf. figs. 2 and 8). Crystals which collapse 
in this way to give a well-developed central fold show few, if any, other 
pleats in the {110} sectors. This is to be expected as the excess material 
in a sector is already taken up in the central fold. The angle of the 
triangular fold (about 5° in practice) can be calculated subject to the 
assumptions (2) and (3) above. The best fit is given by assuming planes 
{312}, and {313}, for the pyramidal surfaces, the former giving 63°, the 
latter 3°. These are the two planes which agreed best with the x-ray 
evidence mentioned previously, which gives confidence in the argument. 
developed here. 

In the foregoing, the chains were assumed to lie along the axis of the 
uncollapsed pyramid. If this condition is relaxed, a variety of other 
mechanisms could be constructed. However, as seen, the simplest 
assumption suffices to explain the experimental observations. 

It needs pointing out that diffraction effects in the crystals undergo 
pronounced changes during electron irradiation apart from their tendency 
to disappear altogether. The sequence will not be described here, it 
will only be mentioned that the selectivity, where it is originally present, 
gradually decreases and the crystal becomes uniformly diffracting. One 
characteristic stage, frequently encountered unless great precautions are 
taken to minimize the beam damage, is when all {hk0} reflections are present, 
with (200), still very weak, split azimuthally, while other reflections are 
split both azimuthally and radially (corresponding to a change both in 
orientation and spacing). This split decreases towards 6;* in the dif- 
fraction pattern, (020) being single. Such a pattern was first brought to 
our notice by Mitsuhashi (private communication) and was quoted as 
evidence for a sheared, non-orthorhombic subcell due to the presence of 
the fold along one of the cell diagonals (Bassett e¢ al. 1959). ‘The existence 
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of the pyramid alone with the chains parallel to the pyramid axis (origi- 
nating in a way such as that proposed by Reneker and Geil, see above) 
requires merely that the true unit cell should be of lower symmetry than 
orthorhombic, i.e. that it should possess an oblique structure. The 
departure from exact orthorhombic symmetry of the subcell itself must 
be inferred from splittings of diffraction spots like those just quoted. 
This splitting is particularly apparent in such patterns. As the selectivity 
is much reduced, each sector contributes with every main {hk0} reflection, 
consequently, the non-equivalence of (hk0) and (hkO) is apparent by 
inspection. More complex splittings have also been observed during 
electron bombardment. In particular, triple {400} reflections are charac- 
teristic of truncated lozenges. The third spot is probably due to the 
{100} sectors. The detailed effect of electron bombardment, though, will 
neither be described nor discussed further here. 

According to the above explanation of the roof ridges, the molecular 
orientation should be preserved in them and no essential change in the 
diffraction conditions for the strongest {hk0} reflections and hence no 
Bragg fringes would be expected. As pointed out in Part I, however, a 
small disorientation, of order 1°, would be sufficient to produce such 
fringes if the crystal as a whole is near diffracting position. It was shown 
that such small disorientations should occur particularly if the crystals 
are thin which is always the case with the single layer polymer crystals. 
Also, slight changes, in the chain orientations due to possible. molecular 
readjustments within the ridges (see Part I) may occur. Thus the 
appearance of extinction fringes in the ridged crystals is not incompatible 
with the basic mechanism proposed in this paper. Large disorientations, 
however, would certainly not be expected. Whenever these accompany 
the striations as, for example, the rotations around the b axis in the 
diffraction patterns of heat treated crystals (Bassett and Keller 1959) the 
present explanations cannot suffice any longer. It will be recalled that 
similarly to the analogous cases in paraffins (Part I) in such cases definite 
pleats are seen in the crystals, which are different from the regular ridges 
although they might be initiated by them. This folding of the crystals 
and the accompanying diffraction phenomena promise to be of significance 
in connection with spherulite formation and the heat relaxation of 
oriented systems (Bassett and Keller 1959) but is beyond the scope of the 
present discussion. 

The interpretation of the ridges was based on the idea that different 
crystal phases with the same subcell can exist in polymers. This situation 
has been well known in paraffins even if the morphological consequences 
in question have not been pointed out previously. In polymers, however, 
such a concept was of no meaning while the chemical repeat unit, which 
is defined as the subcell in the present discussion, was considered as 
the true unit cell. It becomes meaningful only with the introduction 
of chain folding. This makes the crystallography of polymers, poly- 
ethylene in particular, analogous to that of paraffins, the folds in the 
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former corresponding to the end groups in the latter. There is a distinc- 
tion though; the length of the paraffin chain is invariable, while the fold 
length in polymers can be changed by physical factors How far this 
difference may affect the phase changes in question still remains to be 
seen. 

In paraffins the existence of different phases is normally assessed by 
diffraction methods. In principle, this should also be possible with 
polymers. The fold length can, in fact, be assessed from x-ray diagrams, 
the corresponding periodicity showing in discrete low-angle reflections 
when aggregates of single crystals are examined (Keller and O’Connor 
1957). Such reflections at low angles have been known to exist in bulk 
polymers, particularly in drawn fibres, for a long time (e.g. Hess and 
Kiessig 1944) where they correspond to a periodicity along the fibre axis, 
i.e. along the chain direction. These are not normally considered as 
belonging to the crystal structure and are usually believed to be due to an 
alternation of amorphous and crystalline regions. Such explanations 
are by no means completely satisfactory and we are under the impression 
that the whole phenomenon is even qualitatively still unexplained. Once, 
however, such large periods have been associated with chain folding in 
the particular case of single crystals it seemed to be natural for us to look 
for a similar explanation also in other cases where such low angle 
reflections were found. This view, however, is not generally shared, and 
owing to the difficulty of obtaining morphological evidence from bulk 
samples it cannot readily be put to test. Keller and O’Connor (1957) 
found that a film of sedimented single crystals of polyethylene gave a 
‘ fibre-diagram ° corresponding to a chain orientation normal to the film, 
and small angle reflections on the meridian when the x-ray beam was 
parallel to the plane of the film, as would be expected for a vertical 
structure. If the structure became oblique through the formation of 
ridges each single low-angle diffraction maximum would be expected to 
split into two spots along a layer line, the separation of the spots becoming 
larger with the increasing obliquity of the structure. It would be difficult 
to test this directly as individual crystal layers cannot be used and in the 
case of aggregates the nature of the individual crystals is not under control. 
However, analogous effects are well known in the case of oriented fibres. 
Here, under the effect of various physical treatments, the meridional low- 
angle reflections can split along the layer lines while the chains remain 
parallel to the fibre axis. The nature of this so-called ‘four point 
diagram ° (Meibohm and Smith 1951) is still quite unexplained apart from 
the recognition that an oblique arrangement of large periodicities is 
required (Statton 1959). It is tempting to suggest therefore that vertical 
and oblique phases in the sense discussed in these papers might exist in 
the fibres and would gives rise to these low angle diagrams. This would 
naturally imply that the folded configuration would be present in the 
bulk specimen, a conclusion which is suggested by other morphological 
evidence (Fischer 1957). 


358 On Habit Features in Crystals of Long Chain Compounds 


To summarize, it has been shown how very few assumptions, about the 
packing of the folded chains, namely the simplest modes of staggerings 
within and between the ribbons within the crystal sectors, are sufficient 
to relate a great diversity of experimental observations on polymer single 
crystals. 
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ABSTRACT 


The electrical charge flow associated with an inhomogeneous plastic 
deformation, as first described by Fishbach and Nowick, was further 
investigated. Microindentations were made in order to localize the 
deformation to a single small volume. In these conditions the observed effect 
corresponds qualitatively to a transport of negative charge by moving disloca- 
tions. Small additions of cadmium ions enhance the effect slightly. Flexion 
gives strong but non-reproducible signals. 


§ 1. INTRODUCTION 

Fishbach and Nowick observed an electrical signal when rocksalt crystals 
were indented, i.e. inhomogeneously deformed (Fishbach and Nowick 
1955, 1958). We examined this effect further in relation to experiments 
carried out in this laboratory about charged dislocations (Amelinckx et al. 
1959, Remaut et al. 1960, Remaut and Vennik 1961. Pure and 
doped rocksalt crystals, grown in the laboratory, have been used and some 
tests with crystals submitted to flexion were also performed. As far as 
pure crystals are concerned, our observations are different from those of 
Fishbach and Nowick; the origin of this discrepancy could be found. 
In this paper, we describe some of the results which have been obtained 
so far. 


§ 2. EXPERIMENTAL 


Two Leitz microindentation hardness testers were used. The first 
was a standard type instrument, the second had some special features. It 
had a stainless steel point indenter which was electrically insulated from 
the rest of the instrument and could be connected to an outside circuit 
without hampering its movement. The shape of this steel point was 
similar to that of the diamond identer of the current type apparatus. 
Both instruments allowed controlled loading of the indenter and micro- 
scopic observation of the produced pattern. In most of the experiments 
the applied weight was 100g, the indentation process slow and of 
the order of 10,/min. By indenting etched crystals, which were 
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etched again afterwards, we made sure that only the typical rosette 
patterns were produced. ‘The set-up and the electrical circuit is shown 
in fig. 1. Two alternative arrangements were used. In one a 
gold film, evaporated on one of the faces of the crystal, was connected to 
the grid, and through Rg to the earth. In the other, the probe was the 
brass holder of the diamond point, or the steel point. 
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(a) Electrical circuit; (6) recorded signal; (c) equivalent circuit. 
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The crystal was fixed with araldite to a rigidly mounted insulator. All 
necessary precautions to avoid effects of humidity and other possibilities 
of disturbance were taken. The experiments were performed in a small 
dry box, which also acted as a Faraday cage. 


§ 3. REsuLts 
3.1. Pure Crystals 
A considerable number of experiments were made to ascertain that 
the signal was not due to other causes than the plastic deformation. 
The gold film was evaporated on different faces of the crystal and the 
distance of the indentation point varied with respect to it. In this way 
a variety of signals was recorded; a typical one and the arrangement 
used are shown in fig. 1. This type of signal is reproducible in the stated 
conditions, provided that the distance of the indenter to the edge of the 
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crystal is greater than 0-5 mm. In principle, the indentations were always 
made near the centre of sub-grains. The M shape is due to a slowing 
down of the indenter. This only happens when the overall speed of 
penetration of the point is slow, otherwise no bump is observed. That 
. the M shape is due to this could be checked by performing a similar 
indentation test on an Instron machine. The recorded * stress—strain ° 
type curve showed two sharp changes of slope corresponding to the 


Fig, 3 
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(c) 
(a) Sharp indenter, crystal deposited on a flat base; (b) flat indenter, crystal 
deposited on a flat base: (c) flat indenter, crystal deposited on a rough base. 
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decrease and further increase of the electrical signal. A similar M-shaped 
signal is recorded when the indenter is the probe (fig. 2). However, the 
sign of the induced charge is now reversed and its magnitude somewhat 
greater. As seen in fig. 1, the signal is always negative when the indenter 
is not the electrode. This is opposite to what Fishbach and Nowick 
observed. In order to trace back the origin of this, we tried to reproduce 
their arrangement. The crystal (thickness 2mm) was deposited on a 
flat steel base acting as probe (fig. 3(a)). A flat indenter (6=2 mm), 
connected to the earth was used, the applied load was 200 g (fig. 3 (b)). 

As we presumed that the deformation was not localized around the 
indenter in this case, the crystal was placed on a rougher base (a grid as 
used in electron microscopy) in order to enhance this effect. Figure 
3 (c) shows that the signal is now completely reversed. 


3.2. Doped Crystals 

A number of crystals, grown from melts to which quantities ranging 
from 0-09% to 0-6% of CdCl, were added have also been examined. The 
charge flow, obtained in similar conditions and normalized for indentations 
of equal sizes, rises with the cadmium content; the increase is however 
small. As an example we quote the values, obtained for pure rocksalt 
and the crystals containing the highest percentage of cadmium. They 
are respectively 1-39 + 0-13 and 1-664 0-05 ¢ 10-4. 


3.3. Bending 


As a complement to these experiments, several bars of NaCl with a 
thin film of gold deposited on one face acting as active electrode have been 
bent. Very strong signals have been observed of positive and negative 
sign. The magnitude of the signal changes from experiment to experi- 
ment. When the bar is deformed slowly and the deformation stopped 
when the electrical effect starts, new slip bands are always observed. When 
deformation was stopped before an electrical signal could be observed, 
no movement of dislocations or creation of new slip bands could be 
found on subsequent etching. 


§ 4. CONCLUSIONS 


From the described experiments we can deduce that : 

(i) A negative charge flow is observed when dislocations, formed on 
one side of the crystal, move towards the probe. Expanding loops 
thus disturb the electric field in the swept area. Whether the dislo- 
cations become charged during their motion or not, needs further 
investigation. 

(ii) The results of Fishbach and Nowick are due to a deformation 
involving creation of dislocations on opposite faces of the crystal. In 
this way, a mixed effect is observed, the dislocations running in opposite 
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direction to those formed around the indenter play a more prominent 
role, because they are nearer the probe (cf. figs. 1 and 3). 

(iii) Addition of very small quantities of divalent ions influences the 
effects. 


Experiments, now in progress, concerning the effect of temperature, 
and simultaneous measurements of ionic conductivity may allow the 
interpretation of the described results. 
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ABSTRACT 


Emission band forms of several elements obtained with a photon-counting 
spectrometer are compared with previously published spectra. The sodium 
L,;,J(#) vs \ curve departs from the free electron parabola at energies less than 
about 1 ev from the Fermi limit in accord with the theoretical form. The 
sodium extrapolated band width is 2-6+0-3 ev. In beryllium K there is 
no contribution of p states from the bottom of the overlapping second zone. 
Zone overlap is estimated at 1-2 ev. The extrapolated band width for boron 
K is estimated at 14:1+0-8ev. The L,, band of silicon differs from those 
previously published but is apparently not affected by the temperature or 
the purity of the silicon. The extrapolated width is 13-4+0-8ev. A strong 
lne at 2534 and a blue fluorescence which accompanies the lithium K 
spectrum are shown to be connected with the reaction with lithium of 
hydrocarbon molecules not condensed by the cold trap. 


§ 1. INTRODUCTION 


In addition to the material already published on the soft x-ray emission 
spectra of Li, Al, Mg and K (Crisp and Williams 1960 a, b, Crisp 1960) the 
photon-counting spectrometer (Fisher et al. 1958) has provided new 
information about the band forms of the K or L,3 spectra of several other 
elements of low atomic number. The direct registration of intensities, the 
use of low-beam currents to reduce the power dissipated in the target to a 
few watts, and the ability of the photon-counting instrument to follow 
changes of intensity with time have made possible confirmation of features 
not always detectable photographically, as well as corrections to some of the 
band forms obtained previously by photographic photometry. The 
instrumental window, or resolution, in all cases was 0-9 A on the previously 
adopted criterion (Fisher 1954). All the spectra reproduced here, except 
those for Li, were obtained with a glass grating. The pressure in the target 
chamber was 1-2 x 10-* mm Hg and the target was almost surrounded by a 
copper baffle cooled to about —110°c by external conduction to a liquid 
air reservoir. 
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§ 2. SopIuM 


The L,, spectrum of Na at 110°x, published by Skinner (1940), shows a 
departure from the free electron parabola just before the high-energy limit 
isreached. Skinner attributed this feature to zone overlap, but could not 
observe it at 300°x. The Na spectrum published by Tomboulian (1957) 
does not show this feature, and it has been tacitly assumed that there is no 
departure from the free electron shape. Figure 1+ shows a sodium Ly, 
spectrum with a distinct departure from the free electron shape near the 
high-energy limit. This feature can always be obtained from a fresh 
target cooled to about 300°k, but is less apparent when the target surface is 
contaminated. 

A curve of I(£) cc £4? has been fitted on a wavelength scale to the middle 
portion of the spectrum to indicate the extrapolated band width, 2-6 + 0-3 ev. 
To the high-energy side of about 420A the spectrum diverges from the 
E12 curve, following the trend of the theoretical V(Z) vs # curve proposed 
by Jones and Mott (1937). Their curve is included in the diagram, with a 
broken vertical line indicating the approximate position of the Fermi limit. 
This interpretation accords with other evidence that the Fermi surface in Na 
is nearly spherical and that as a consequence the occupied states do not 
come closer than about 1 ev to the nearest zone boundary. 


§ 3. BERYLLIUM 


Measurements on fourth-order K spectra of Be show that the metal 
edge extends to between 25 and 30% of the intensity at the peak. The 
width of this edge estimated from measurements on second, third and fourth 
order spectra is about 0-2 ev which agrees with Skinner’s (1940) result since 
the temperature of the surface of a solid Be target is probably appreciably 
higher than that of the cooling water. According to Skinner the contribu- 
tion of states in the second Brillouin zone “ can just be detected as a slight 
upward curvature at the edge”. None of our Be K spectra show any such 
upward curvature, though the upper limit of the metal edge is quite clearly 
marked. ‘This is taken as evidence that there are no low lying p states 
in the second zone. The prominent spike in the MgL,, spectrum shows 
that s states are present in the bottom of the second zone, and this can be 
assumed to hold for Be. By extrapolating the J(Z) vs A curve of Be past 
the metal edge, the area corresponding to the unoccupied part of the zone 
can be estimated approximately. This part of the zone is not filled 
because s states in the overlapping second zone are occupied. If a 
hypothetical J(H) vs A curve is drawn (with I(H) oc H3!) for these s states, 
having the same area up to the metal edge, the amount of zone overlap can 
be estimated as between 1 and 2eyv. 


ee eee 
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§ 4. Boron 

A fourth-order boron K spectrum is reproduced in fig. 2. The points 
shown have ordinates proportional to J(#)2/%, so if for the bottom of the 
band /(£) oc H®?, the intersection of the straight line fit with the axis should 
define the lower band limit. The ‘band end’ recorded by Skinner (1940) 
on the high-energy side involves intensities less than 1° of those at the 
peak and the existence and extent of this ‘band end’ are both uncertain in 
our spectra. Assuming that the band end extends to 65-5 4, the extra- 
polated width of the boron K band is then 14:1+0-8ev. Thisis comparable 
with the 17ev obtained by Coster and Hof (1940), but is much less than 
Skinner’s estimate of 30+2ev which included a long tail as well as the 
‘band end’. The short-wavelength fall of the band shown in fig. 2 is 
distrinctly concave towards high energies, not straight as other published 
spectra show. Its width is 2-6ev (instrumental width 0-6ev) which is 
consistent with emission involving the valence band of a semi-conductor. 


§ 5. SILICON 

Silicon L,, spectra from evaporated Si have been published by Skinner 
(1940) and from solid crystalline Si by Coster and Hof (1940). The Si 
spectrum in fig. 3. was obtained from a target of powdered AR grade Si 
which glowed red-hot under the action of the electron beam. An identical 
spectrum was obtained from a single crystal of 30-ohmecm p-type Si. 
The temperature of the crystalline Si varied from about 100°c for the 
lowest target powers (2Ww) to possibly several hundred degrees for the 
spectrum of fig. 3, but no changes with temperature have been observed. 
The spectrum of fig. 3 resembles that of Coster and Hof (1940) in that the 
maxima are the same in number and position, but the indication is that in 
the photograph taken by these authors the emulsion was seriously saturated. 
Skinner’s spectrum, which closely resembles the SiO, spectrum of Coster 
and Hof, shows a single broad maximum which cannot be explained by 
saturation of the emulsion. 

The points plotted are proportional to J(H)?, assuming I(H#)ocH1? 
for the bottom of the band so the intersection of the line with the axis gives 
the lower band limit. The extrapolated width of the band is then 
13-4 + 0-8 ev in agreement with Skinner’s l4ev. 


§ 6. LitHIuM 

A strong line at 253 A accompanying the Li K spectrum was attributed in 
an earlier paper (Fisher et al. 1958) to contamination by kerosene in which 
the solid Li had been stored. Subsequently it has been found that the line 
was produced from an evaporated Li target, when the Li had previously 
been vacuum melted and the target was nearly surrounded by a copper 
baffle cooled to about —110°c. A blue fluorescence is closely associated 
with the appearance of the line and both, needless to say, develop only on a 
lithium surface. A strong evaporation of Li which coats the baffle as well as 
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the target is found to prevent the appearance of either fluorescence or 
line (Crisp and Williams 1960 b). 

Figure 4 (a) shows the decrease with time of the peak intensity of Lik, 
and also, since the band shape does not change, of the intensity in the tail 
of the Li band at 2534. The decrease is roughly exponential, the intensity 
falling to half in about 1}min. The variation with time of the intensity 
at 253 Ais shown in fig. 4(). When the exponential decrease is subtracted, 
there is left a component whose intensity increases roughly as [1 — exp (— 2) | 
which represents the growth of the Li line. A continuous chart record 
showing successive traversals of the LiK band and Liline for 60min 
following the recording of fig. 4 (c) showed that after reaching an equilibrium 
value the line intensity remained constant, while the Li K peak intensity 
fell by 90%. 

The fall in intensity of the Li K band is mainly caused by the deposit of a 
layer of carbon on the target following decomposition of condensed hydro- 
carbon molecules by the electron beam. The exponential fall indicates 
that the carbon film thickness grows at a constant rate. Carbon con- 
tamination of Mg (Crisp 1958) was almost eliminated when the baffle was 
cooled below a critical temperature of — 90°C, indicating that a substantial 
proportion of the hydrocarbon molecules were removed by condensation on 
the cooled surfaces. However, fluorescence, carbon contamination, and 
growth of the Li line were still observed for Li targets. While it is possible 
that Li,O and Li,N could be formed on the target, both would be decom- 
posed by electron bombardment, and Catterall and Trotter (1959) have 
observed this for Li,0. It seems more probable that hydrocarbon mole- 
cules which do not condense on the cooled surfaces react with Li, but not 
with Mg. 

The discovery that a heavy evaporation of Li inhibited development of 
the Liline supports thisexplanation. When the baffle surfaces surrounding 
the target are coated with Li this will help remove hydrocarbons from the 
residual vacuum and less of these will be available with the Li on the target 
to produce fluorescence and the Li line. There is no marked decrease in 
pressure which would be expected if nitrogen were the reactant, but there is 
a marked slowing of the rate of decrease in intensity of the LiK band, in- 
dicating a much slower rate of deposition of carbon. Figure 4 (d) shows a 
lithium K spectrum in which no line has appeared a few minutes after a 
heavy evaporation of Li. A continuous record of the LiK band and the 
region of the Li line for 100 min following this evaporation showed that no 
line developed and that the LiK peak intensity decreased by only 33%. 
This is one-quarter of the contamination rate for the conditions of fig. 4 (c). 

It is concluded therefore, that removal of hydrocarbons by the widely 
spread Li film, decreases the thickness of the hydrocarbon layer reacting 
with lithium on the target surface to an extent that makes the line unde- 
tectable and the associated fluorescence invisible, and consequently lowers 
the rate of decomposition of hydrocarbon molecules on the target by the 
electron beam so reducing the rate at which the Li K band intensity drops. 
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ABSTRACT 


A method is proposed for calculating correlation energies and pair functions 
associated with an electron gas at any density. Upper bounds to the exact 
energy, lying however at least as low as the corresponding plane wave results 
are guaranteed and exact solutions are obtained in the high- and low-density 
limits. The technique used is that of allowing a periodic electron density 
and a proof emerges of the recent conjecture of Overhauser that there exist 
one-particle states energetically more favourable than the usual plane wave 
orbitals. It is concluded that in the ground state a single-particle approxi- 
mation implies static spin density waves of wavelength proportional to the 
usual separation electron parameter 1s, at least for rs sufficiently large. 


§ 1. LyrropuctTIon 


THE modern theories of the electron gas, for example those of Wigner 
(1938), Bohm and Pines (1953), Gell-Mann and Brueckner (1957) and 
Cowan and Kirkwood (1958), have proved to be illuminating and 
valuable within their own restricted domains of validity. However, 
each has its own peculiar difficulties and it seems fair to say that there 
does not appear to be a simple, reliable method applicable over a wide 
density range. It is therefore of considerable interest to note the recent 
speculation of Overhauser (1960 a, b), and of Gross (1960) that the Hartree 
equations appropriate to an electron gas have solutions of a periodic 
nature corresponding to a total energy below that of the usual plane 
wave solutions. In the following we will prove the latter assertion, show 
that the periodic orbitals are related to Wigner’s lattice and demonstrate 
that the Hartree-Fock method promises to give accurate information for 
an electron gas below a certain density. In particular we show that exact 
results are forthcoming in the limits of very high and very low densities. 

The technique used is the variational principle employing a Slater 
determinant of orbitals corresponding to the von Weizsacker approxima- 
tion (March and Young 1958). The latter orbitals are capable of 
achieving the high and low density results mentioned in the previous 
paragraph. Should it prove necessary, a discussion using more general 
orbitals is entirely possible but it is much easier to follow the behaviour 
of the various dynamical functions as the density is widely varied in 
the present approximation. The use of the von Weizsacker orbitals 
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also reveals an interesting relation with the results of Young and March 
(1960) (subsequently referred to as YM) from which indeed, the following 
work stemmed. YM used a method suggested by Mayer (1951) and specified 
a second-order density matrix characterizing the system. The advantage 
of the method was its wide applicability, its disadvantage was that the 
density matrix was not obviously derivable from a many-body wave 
function. In this respect the following constitutes an improvement over 
the YM investigation. 

The work is divided into two parts, § 2 dealing with a one-spin and 
§ 3 dealing with a two-spin electron gas. It has been thought worth while 
to present the ferromagnetic theory since the effect of allowing a periodic 
electron density alone may thereby be studied. In § 3, the additional 
complication of the juxtaposition of the two kinds of spin density waves 
is resolved. 


§ 2. THe ONE-SPIN GAS 
It has been shown (Mayer 1955) that the investigation of a degenerate 
electron gas for all densities 3/4775* is equivalent to considering the ground 
state of the Hamiltonian 
[he 1 1 Nc. ade, dr,dr 
Pa a2 yy a+) pn | patie ta os + { [Soe ‘| 1 
n> aax lear = [evr] J J Iron i 
of a plasma of N particles in a volume JN for all values of the coupling 
parameter A= (47/3)'rs. The appropriate energy (in atomic units) is 
then 


EF ES Re ar tater : hen! fl—2Pi(r'r)],., ‘ 
aS =a [V,2y(r’ | r) Je ord 5 | [Rae dr, (2) 


where y and P are respectively the usual first-order spinless density matrix 
and pair function (for precise definitions see YM). For further details of 
the density matrix formalism reference may be made to Léwdin (1955), 
Mayer (1955) and Koppe (1957). Formula (2) has been used previously by 
Mayer and YM under the usual Sommerfeld condition of constant particle 
density but it is readily seen to retain its validity under the more general 
periodic conditions. 

We now recognize the desired density periodicity in the following way. 
Suppose we consider the high coupling limit A oo. Then by an argument 
of Wigner (1938) the potential energy becomes all-important and the 
particles take up positions on a body-centred cubic lattice. On relaxing 
the coupling we obtain a periodic particle density and we are now free to 
employ the conventional methods available for such problems. Strictly 
the particular choice of lattice should be chosen on energy grounds at any 
given coupling, but as this seems impracticable we will restrict our con- 
siderations to the body-centred cubic case. It should be remembered 
however, that an upper bound to the energy below the plane wave value is 
obtained whatever lattice is used. . 
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We now choose as trial wave function the single Slater determinant of 
Bloch orbitals 


1/2 
bu(r) =" : sous RS eta 


a NA 

= iB See) aris 
which have been discussed by March (to be published). Here the k are 
the NV vectors of least moduli having integer coordinates and J is the 
Jacobian of R(r), the latter having the properties 


Ries BS Rir) RS) = Rr), hk (4) 


y being any translation vector between equivalent lattice points. R is 
now available for energy minimization using the variational principle. 
Computing matrices in the usual way we obtain 


eit (r)) (rn) Re Rey Riv) 3 2 «\ cin 20 (5) 
and 
2P(r’r)=J(r’) J(r) (1—I{R(r’)—R(r)}] 2... (8) 
where 
U(r) = 32-3 (sin€—Ccosl), f= (672). PRY een Te) 


Hence, using the periodicity condition (4) we obtain from (2), 


€ 1 6n*)? 7S op (VI |g WGI ve 
s=5{] 10 1.210%) nek aif |er | oe ee E 5 (8) 


where 
F(é)= | 2P(r+Erjdr= | Te+B) T(r) LPR +E)—R(P)}Idr, (8) 
ie 


C denoting a unit volume Wigner-Seitz cell defined by the body-centred 
cubic lattice. From (5) we see J plays the réle of the particle density. 

It is now possible to see that exact descriptions are possible in the two 
density limits using (8). When A 0 only kinetic energy is important and 
R=r gives the usual plane wave result. When A-—- co only potential energy 
is important and the Wigner lattice of point charges is given by 


J => d(r—X,), EaLAbe al shamcueratte ome LO) 
the X; denoting all the lattice points. (To obtain (10), one defines R(r) 
within C by the coordinates of the point where the straight line from the 


centre of the cell to r meets the cell boundary, (4) being used to extend the 
definition to all space.) The average pair function then becomes 


HG Nae Orie on lt A) 


the exchange term cancelling the charge at the origin since R(0) =0.. Use 
of (11) in (8) then gives the exact energy per particle in this limit of » 


-3({ By \e-— eae 5 (15) 
2A\ J Alspace. & Xx, +0X 4 Ps 


the latter following from a Madelung calculation (see Wigner 1938). 
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Since, for large r,, (12) is lower than the plane wave result, we have proved 
Overhauser’s conjecture for the ferromagnetic gas, at least for sufficiently 
low densities. 


§ 3. Tur Two-sPIN CASE 


In setting up this case we must consider the way in which the two kinds of 
particles are to be associated with the Wigner lattice. In order to obtain 
the desired results, it seems necessary to consider the body-centred cubic 
lattice as consisting of two interpenetrating simple cubic sub-lattices, 
each associated with a given spin. In coming to this conclusion the 
results of §8 of YM were particularly useful and reference to fig. 4 of that 
article should prove helpful to the reader. 

Accordingly we describe the situation by a determinantal wave function 
of N/2 orbitals J¥?/N1? exp(27i)(k.R/N1%) associated with one spin and 
N/2 orbitals K4?/N1? exp(271)(k.$/N 3) associated with the other, the k 
consisting of the V//2 vectors of smallest moduli having integer coordinates 
and J and K being the Jacobians of R and S respectively. We also require 
the periodicity conditions 

R(r py.) = 9S (6) 4 ea ee ad 
and 
R(r+p.)=R(r)+y.,, R(—r)=—R(r) . . . . (14) 
where p, and yw, are translation vectors between equivalent points in the 
body-centred and either simple cubic lattice respectively. The spin 
densities are $J and 3K, the total density being }(J +). We obtainina 
sear eewacd manner 


ye’ |r) = 3 LF (1) J (r)} 12 Re’) — R(r)} + Bf K(r’)K(r) 12 US(r’) — S(r) } 


and 
2P(r’r)=f{J(r') + K(r')} {J (r) + K(r }} — gd (r') J (r)P{R(r’) — R(r)} 
—ZK(r')K(r)P{S(r’)—S(r)} (16) 


where here 
Ur)=3f-* (sing = cost), C= (8°) ek ea) 
Equation (2) is again applicable to give 
@imal a) aie 1.(VJ)? 1 l—-F 
Pee Z dr—1 erst), 
N al, (ee 2 2(V a6 J } : al m space = ve re 
where 


(=| Pie+end rab [ (er +8)I( r)[1—2{R(r+€) — R(r)}] 


+J(r)K(r+ &)\ dr, BS ec ig ee) 
D denoting two adjacent Wigner—Seitz cells defined by the body- -centred 


cubic lattice. The apparent lack of symmetry with respect to R and $ in 
(19) is due to the use of (13). 
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At high densities R=r gives the plane wave results, while at low densities 
a suitably defined step function for R leads to a distribution of point charges 
on one of the simple cubic structures for $J. An analysis somewhat 
similar to the one-spin case then reproduces the Wigner result (12). 


§ 4. Discussion 


We have seen above how one may profitably modify the usual one-particle 
model of an electron gas to obtain a scheme which is rigorously based on 
the variational principle in such a way that at all densities an upper bound 
to the exact energy at least as low as the plane wave value is assured. This 
upper bound property is by no means a trivial consequence of the use of 
the variational principle, as it is well known that a trial wave function 
chosen without due regard to the physical conditions of the problem can 
lead to energies which are too low. Since furthermore it has been shown 
that exact results are obtained in the high- and low-density limits, one 
might hope for accuracy at intermediate densitiest. While detailed 
numerical calculations in the region corresponding to real metallic densities 
are not as yet available, work is currently being undertaken in this Depart- 
ment in which the aim is to obtain the accurate predictions of this 
approach and that of YM over a wide density range. 

A point of interest implied by the results of § 3 is the form of the spin 
density waves in the ground state. It is felt that by relating them to the 
Wigner lattice, as described above, we obtain a more accurate description 
within a one-body framework than that of Overhauser (1960 a, b), useful 
though the latter was in demonstrating the existence of new low-energy 
states. The existence of such spin density waves in the ground-state 
electron gas would appear to lend support to the theory of 
antiferromagnetism proposed by Overhauser (1959), although clearly the 
quantitative predictions of this theory will have to be reworked{. It 
is hoped to report further on this at a later date. 

The connection of the present approach with the conventional philosophy 
of the Sommerfeld model also merits attention. The Sommerfeld method 
may be summarized as saying that since we have a translationally invariant 
Hamiltonian we may suppose all points in the gas to be equivalent in the 
sense that the wave functions differ by at most a trivial phase factor indepen- 
dent of the electron coordinates. It follows that the Nth order and all 
reduced matrices are translationally invariant within the Sommerfeld 
framework. If we consider for notational convenience the ferromagnetic 
case, it is readily seen that eqn. (8) results fom the use of the translationally 
invariant Nth-order trial density matrix 


WP meey UKE LEE oN EBM E rat Bot B)GE, (20) 
Volume V 
+ See, however, the recent work of Koun, W., and Netret, S.J., 1960, Phys. 


Rev. Letters, 5, 8. 
{ See, however, MarsHa1, W., 1960, Phys. Rev., 118, 1519. 
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where 7 is the Slater determinant of the orbitals (3). Thus while many 
methods (for example Jastrow (1955) and Clarke and Feenberg (1959)) use 
variationally trial wave functions which are translationally invariant, we 
have seen it is permissible to use the Nth-order matrix (20). The mathe- 
matical simplicity thereby introduced is in marked contrast with the method 
based on wave functions. The interpretation of the matrix elements will 
however be different in the two cases. For example, the Sommerfeld second 
order matrix based on the use of (20) is 


i v(t; +6 ir +6) y (ry +8| r.+6) . 
Tl eaaey Hole LE |r+2) Sieben 


y being defined by (5), the pair function thus being 


1 
mal 2P(E+nn)dn= | 2P(§ +n) dn. 
Volume NV Cc 


More generally, diagonal elements of matrices on the Sommerfeld model are 
obtained by averaging the corresponding periodic elements. This is in 
accord with the essential distinction between the methods, the one placing 
all points on an equal footing, the other inducing a periodicity by reformula- 
ting the problem. 

Finally, it is interesting to note the degree of similarity between the 
results presented here and those of YM despite the use of fundamentally 
different approaches. This is perhaps best illustrated by an examination 
of the pair functions. From (6), for example, we may say 


J (r+ §)J(r)[1 —P{R(r + §) — R(r)} ] 
is the probability that there are simultaneously particles at r and r+ &. 
Since J(r) is the particle density, J(r+&)[1—/{R(r+ &)—R(r)}] is the 
probability that there is a particle at r+ § when there is a particle at r. 
Since J is greatest at the lattice points, the most likely value of 


J(e+ B)[1—-P{R(r + 8)—R(r)}] is J(E)[1—-P{R(E)}] 
and this is just the YM Sommerfeld pair function. Clearly a complete 
elucidation of the connection between the methods would be of interest. 
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ABSTRACT 


The principal features of the experimental results of Blewitt et al. (1957) 
and of Suzuki and Barrett (1958) on deformation twinning in face-centred 
cubic metals are explained in terms of an extension of the ‘ prismatic ’ source 
mechanism of Cottrell and Bilby (1951). Suitable prismatic dislocation 
sources are thought to be long jogs in the dislocations of the two conjugate 
slip systems in a work-hardened metal and prismatic dislocation loops 
produced by vacancy condensation in a neutron-irradiated metal. The 
dissociation of slip dislocations is also considered and is found to be an unlikely 
mechanism for the production of deformation twins in copper, silver or gold. 


§ 1. INTRODUCTION 


Brewirt ef. al. (1957) have observed twinning as a mode of deformation 
of copper single crystals in low-temperature tensile tests. Suzuki and 
Barrett (1958) have made similar observations on a series of silver—gold 
alloys. 

Figure 1 indicates the principal features of the results on copper of 
Blewitt et al. (1957). Twinning, when observed on the principal slip 
plane (111), took place by shear in the [112] direction and not by shear 
in the [211] direction. Twinning was not observed at 4:2°K when the 
tensile axis was between [311] and [100], and was confined to a small 
range of orientations near [111] at 77°K. In specimens which had been 
neutron-irradiated to an integrated flux of 1-5 x 10°°nvt at 80°o, twinning 
at 4-2°K occurred on the secondary slip plane (111), although the stress on 
the primary twin system (111)[112] was higher than on the conjugate 
system (111)[121]}. Specimens which had been irradiated to lower 
fluxes at 30°c twinned at 77°K in a similar manner to work-hardened 
crystals at the same stress. Similar results were obtained for work- 
hardened silver—gold alloys by Suzuki and Barrett though twinning often 
occurred on both the conjugate slip planes (111) and (111). 

The resolved shear stress for twinning is high, typically ~ 15 kg/mm? for 
Cu and ~6-10kg/mm? for the Ag—Au alloys. The onset of twinning is 


+ Communicated by the Author. 

+ Blewitt et al. (1957) state that for crystal 416, whose tensile axis had not 
rotated to the dodecahedral plane when twinning began, the stress is higher on 
the secondary slip plane. This is not so for the two twin systems quoted above. 
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accompanied by a load drop of about 10% of the apphed load, but as the 
area of the twinned region is typically about 10% less than that of the 
non-twinned region of the specimen, the shear stress for the propagation 
of twinning throughout the crystal is about the same as that for twin 
initiation. 


Fig. 1 


© CATEGORY A 
®@ CATEGORY 8 
S TWIN AT 77.3°K 


© DID NOT TWIN AT 
77.3°K 


(Reproduced by courtesy of the American Institute of Physics.) 
Stereographic projection showing the tensile axes of the copper crystals tested 
by Blewitt e¢ al. (1957). Crystals of class B twinned in tensile tests 
at 4:2°K whereas those of class A did not. The shaded area of the 
stereogram contains those tensile axes which had rotated to between 

[311] and [111] after 50% extension of the crystals. 


Any model of deformation twinning must account for the above results, 
in particular for the twin systems observed and the orientations of the 
tensile axis necessary for twinning to occur. In this paper a dislocation 
model of the production of deformation twins is proposed and the results 
of Blewitt et al. (1957) and of Suzuki and Barrett (1958) are interpreted in 
terms of the proposed model. 


§ 2. PosstsLz Moves or Twinnine in F.C.C. Merars 


If a twin is generated by the spiralling of a twinning dislocation around 
a suitable pole, then sources of the twinning dislocation, and a pole dis- 
location, whose component Burgers vector perpendicular to the {111} 


twin plane is equal to the spacing of these {111} planes, must be present 
in the metal. 
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It is convenient to use the reference tetrahedron of Thompson (1953) 
to designate the slip planes and Burgers vectors in the f.c.c. lattice. The 
orientation of the tetrahedron is shown in fig. 2. In work-hardened 
copper, with the tensile axes of fig. 1, the principal slip and twin plane is 
(a) and the observed twin direction aC. The twinning dislocation aC can 
be produced in the dissociations of the whole dislocations BC, DC and 
AC. BC and DC may be called ‘glide’ sources, as they are glissile in (a), 
and AC a ‘prismatic’ source as it cannot glide in the twin plane. Suitable 
pole dislocations for twinning on (a) are AB, AC and AD. 


OlO 


Thompson’s notation for slip planes and Burgers vectors in a face-centred cubic 
crystal. Planes (a), (5), (c) and (d) (opposite corners A, B, C and D) of the 
reference tetrahedron correspond to (111), (111), (111) and (111). The 
mid-points of the faces (a), (6), (c) and (d) are denoted by a, f, y and 6. 


Cottrell and Bilby (1951) have considered the dissociation of the 
prismatic dislocation AC (fig. 3 (a)) on plane (a) into the Frank sessile Aa 
and the Shockley partial aC. Frank (1951), and Frank and Nicholas (1953) 
have designated the two possible types of stacking fault in face-centred 
cubic metals ‘intrinsic’ and ‘extrinsic’ fault. It will be assumed through- 
out, with Cottrell and Bilby, that intrinsic faults occur and that they 
have a much lower energy than extrinsic faults. There is now considerable 
experimental evidence that this is the case (Howie 1960). The stacking of 
(111) planes across the intrinsic fault behind a€ has the sequence 231/312 
(fig. 3(b)). After one revolution around the pole dislocation, the aC 
dislocation meets the sessile Aa again (fig. 3(c)), and, according to Cottrell 
and Bilby, cannot proceed further as the stacking sequence 231/123 
would be produced, in which nearest-neighbour distances are not preserved. 
However, when aC meets Aa, we have the original dislocation AC and a 
sheet of stacking fault. If this dislocation AC glides prismatically along 
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the glide cylinder (G-G’, fig. 3(c)) in either direction through a distance 
equal to the Burgers vector than the dissociation 


AC >Aa+aC 


can take place again. This dissociation produces the double stacking 
fault sequence 231/3/23 or 23/2/123, depending on whether AC glides up 
(as in fig. 3(d)) or down respectively. Repeated operation of this source 
then gives the required twin. The production and operation of such a 
source is considered in § 4 for work-hardened metals and in § 6 for neutron- 
irradiated metals. 


Fig. 3 
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The extension of the mechanism of Cottrell and Bilby (1951). 
(a) The unit dislocation AC lying in the (111 )planes (a). The positive 
direction of the dislocation is into the paper, along CB. 
(b) The dissociation AC+Aa+aC, producing intrinsic stacking fault on 
plane (a). 
(c) Recombination of the dislocations Ae and aC after one revolution 
of the Shockley partial. 


(d) Prismatic glide of AC along G-G’, followed by the dissociation (b) 
on the next atomic plane. 


It is also conceivable that twinning could originate from the slip dis- 


locations BC and DC by the normal splitting on plane (a) (Heidenrich 
and Shockley 1948), 


BC->Ba+aC and DC-+Da+ac. 
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Under certain conditions the dislocations Ba and Da could be immobilized 
and these reactions are considered in § 3. 

It will be shown below that the separation of glide sources into partials 
is unlikely to provide a suitable twinning mechanism, and that the facts 
recorded in the introduction can be understood in terms of the extension 
of the mechanism of Cottrell and Bilby given above. 


§ 3. OPERATION OF ‘ GLIDE’ SouRcES 
3.1. General Considerations 


If a slip dislocation, Burgers vector b, which splits up into Shockley 
partial dislocations b, and b,, is to be a twin source, a resolved shear 
stress on the slip plane, 6, must be applied which is sufficient to separate 
these partials. The configuration of the dislocation under stress may 
approximate to that of fig. 4(a). The front partial is b,, which is bent to 


Fig. 4 


D B 


Glide sources. 
(a) Configuration of a glide dislocation under stress. 
(b) Burgers vectors of the dislocations which are glissile in plane (a). 


a radius 7, under the stress o, and the back partial b, has a radius of 
curvature r,. The dislocations are pinned at the points P and Q, which 
are a distance / apart. Unit lengths of the two partial dislocations (in 
the centre of the arcs PQ) are in equilibrium if the following two equations 
are satisfied : 


2 
o.b,+F(d)-y- Te ae CERT SIT} 


- 


2 
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F(d) is the repulsive force per unit length between the partials, which 
depends on their spacing d, and wis of order unity for large r. The stacking 
fault energy is y and p is the shear modulus. If the back partial by also 
operates as a Frank—Read source then slip occurs rather than twinning. 
The back partial cannot follow the front partial beyond the unstable 
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semicircular configuration, if eqn. (3.2) can still be satisfied when the 
front partial is at a large distance from the back partial. When the 
separation of the partials is large, #(d)=0, and eqn. (3.2) reads, since 
Sri, 


o.b,+y— 22 <0. 35), ORs) 


Thus twinning may originate from a glide dislocation if eqns. (3.1) and 
(3.3) are satisfied and the dislocation can only be a source of slip if they are 
not satisfied. 

The whole and partial dislocations which are glissile in (a) are shown in 
fig. 4(b), which is the face (a) of Thompson’s tetrahedron. The dis- 
sociations are as follows: 


(1) BC-—> Ba+aC, 
(2) DC-+Da+aC, 
(3) DB ->Da+aB. 


The tensile axis is between [311] and [111] when twinning occurs, and 
hence the direction of greatest resolved shear stress on (a) is between BC 
and aC. By applying the conventions of Frank (1951) and Thompson 
(1953) it may be shown that aC is the front partial in both dissociations 
(1) and (2), given intrinsic stacking faults. The stress is such as to separate 
the partials, independent of the sign of the Burgers vector. The stress in 
case (3), although it tends to move the partials Da and «B in opposite 
directions, tends to decrease the width of intrinsic stacking fault, and so 
the dislocation DB is not expected to be an important twin source. The 
operation of the sources BC and DC is considered in the next section. 


3.2. The Stresses Required to Operate ‘ Glide’ Sources 


The eqns. (3.1) and (3.3) give the restrictions on the stress 6 and source 
length / necessary for twin production from ‘glide’ sources. In order to 
obtain estimates for o and /, the following simplifying assumptions will 
be made. It may be seen by more detailed calculations that these esti- 
mates are not in fact restrictive enough. 

It will be assumed that 6 can be represented by an effective local tensile 
stress 7. Thus o.b, will be equated with 7b, cos y cosA,, and 7b, cos y cos Ay 
with o.b,, where x and A are the angles between the tensile axis and the 
glide plane normal, and the tensile axis and the glide direction respectively. 
Since the stress on the dislocation BC or DC under consideration tends to 
separate the partials, the separation of the partials, d, will be greater than 
the equilibrium value. ‘Therefore F'(d) will be smaller than its equilibrium 
value, y. If the value of F(d) only varies slowly with r,, then the front 
partial b, will be unstable when 2r,=/, in a similar manner to an ordinary 
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Frank—-Read source. Thus eqn. (3.1) now reads, putting F(d) <y 


7b, COS eos, > WA, Ber ee CS) 
and eqn. (3.3) reads 
Son caeea eam ie h = 
1 COS y CosAg, 5 y, where |b,|= |d,|. 4 as) (3:0) 


These equations give the following condition on J when cos), and cosa, 
are both positive: 


2 
pore (1-= Rarer nett veo) ree (56) 


Table 1 gives the minimum values of 7 and the maximum values of | 
necessary to operate the twin sources BC and DC, as a function of tensile 
axis orientation. 


Table 1. Values of the source length, /, and of tensile stress, 7, necessary 
for glide dislocations BC and DC to act as twin sources in copper. 
The observed tensile stress, tops, when twinning occurs is 
~ 35 kg/mm? (Blewitt et al. 1957) 


[y=40erg em, »=4-2 x 1011dynecm-, Db, =1-:5x 10% cm] 


Minimum tensile 


aw Spb i 
ensile Restriction on J, stress, 7, eqn. (3.4) 


axis eqn. (3.6) j (A) (kg/mm?) 


DC 


+ A compressive stress is required to separate the partials. 
t Value estimated from eqns. (3.1) and (3.2) with #(d)=0. 


The main feature of table 1 is that, whereas the DC dislocation source 
requires a lower shear stress than the BC source for operation, separation 
of the partials Da and aC will be most easily effected with a tensile axis 
near [100], when twinning is not observed. The dislocations of the 
principal slip system have Burgers vector BC, and these dislocations may 
react with the few dislocations CD which are present, to form BD dis- 
locations. In this case only the BC dislocations will be available as glide 
sources, and these sources, while giving approximately the observed 
dependence of twinning frequency on tensile axis orientation, require 
very large stresses and prohibitively short source lengths. 


PvE. 2B 
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Cottrell and Bilby (1951) pointed out that the pole dislocation, around 
which the twinning partial rotates, must not move under the stress 
necessary to move the twinning dislocation. It seems very unlikely 
that the dislocation network would stand the stresses given in table 1, 
and unless the back partial can be immobilized in some other way than 
by making the source length small, then it is doubtful whether glide sources 
can operate, except in metals of very low stacking-fault energy. The 
further difficulty of the partials overcoming their mutual stress field after 
one revolution of the twinning dislocation has been considered by Seeger 
(1956), for the case of generating a h.c.p. from a f.c.c. lattice using a pole 
mechanism, and will not be considered here. 


3.3. ‘Glide’ and ‘ Prismatic’ Twin Sources 

The back partial (Ba or Da) may most effectively be immobilized by 
reacting with a dislocation which does not lie in (a), and which cannot glide 
in this plane. This limits the choice of the interacting dislocations to the 
possible pole dislocations AB, AC and AD. If the interactions with the 
sources DC and BC are to have sufficient stability, the Burgers vectors 
of the interacting dislocations must be at 120° to each other (Whelan 
1958), and hence the total Burgers vectors of these new sources must be 
BA, DA or AC. ‘These sources, however, are no longer glide sources but 
prismatic sources, as the total Burgers vectors do not lie in the twin plane 
(a). 

Hence it appears that, on several counts, the slip dislocations are un- 
likely to be effective twin sources. Prismatic sources will be considered 
in the following sections. 


§ 4, Prismatic Twin Sources IN WORK-HARDENED F.C.C. Mprvats 
4.1. Possible Prismatic Twin Sources 


Prismatic dislocations which lie in plane (2) may have Burgers vectors 
AB, AC or AD. These dislocations may become sources of twins on each 
of the three twin systems on plane (a) by dissociating in the manner 
considered by Cottrell and Bilby (1951): 


(4) AB > Aa+oB, 
(5) AC->Aa+aC, 
(6) AD->Aa+aD. 


These dissociations may take place when the stress on the twinning dis- 
locations is high. For tensile axes between [100] and [111] the resolved 
shear stress on (a) is greatest between the Ba and aC directions, and is 
never high in the aD direction, so that the dissociation (6) is unimportant. 
If it is assumed, as before, that intrinsic faults have substantially 
lower energy than extrinsic faults, then the twinning dislocations 
aB and aC separate on one particular side of the Frank sessile 
Aa. (see fig. 5). By using the convention of Thompson (1953) again, it 
can be seen that the stress is such as to expand the intrinsic fault behind 
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aC and contract the intrinsic fault behind «B. The expansion or con- 
traction of the fault is independent of the sign of the Burgers vectors. 
(The difference between reactions (4) and (5) may be seen simply. ‘The 
dislocations #B and aC are on the same side of Aa, but the stresses are 
high for the movement of Ba and aC dislocations in the same direction, 
and hence #B and aC must move in opposite directions (fig. 5).) 


Fig. 5 
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The dissociation of prismatic dislocations. 


(a) AC+Aa+ac. 
(6) AB-Aa+aB. 


The positive direction of the dislocations is along CB, into the paper. 
Under the tensile stress #B will move to the left and aC to the right. 
The Burgers vector aB points into and ac points out of the plane of the 


paper. 


Thus the AC dislocation is the only important source for twinning on 
plane (a) and twinning can only take place by shear in aC direction on this 
plane. The result of Blewitt et al. (1957), that twinning on the (111) 
plane takes place in the [112] direction only, can therefore be explained. 
in terms of the geometry of the twin source dislocations, if only intrinsic 
stacking faults are produced. 


4.2. Production of Prismatic Dislocations in the Work-hardened Metal 

The dislocations of the secondary glide system have the Burgers vector 
AC and lie in plane (b). Lengths of these dislocations can be transferred 
to plane (a) by glide processes in three ways: 

(1) The dislocations can lie along CD. 

(2) They can cross-slip onto plane (d) and then lie along BC, the 
direction which is also in plane (a). 

(3) They can be intersected by dislocations of the primary slip system. 
If Frank—Read sources of ‘plane’ type operate (Bilby 1955), 
repeated intersection of the dislocation AC will produce a long jog 
which lies along the primary slip direction BC in planes (a) and (qd). 

These three processes will be considered separately. 
2B2 
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If the dislocations AC lie along CD they can react with the primary 
slip dislocations, Burgers vector CB, to form Lomer—Cottrell locks AB, 
lying along CD. It seems unlikely that considerable lengths of the 
dislocations will lie along CD unless they are part of Lomer—Cottrell locks 
(see Whelan 1958). 

If process (2) is important then twinning cannot occur below the cross- 
slip stress. The beginning of stage III of the work-hardening curve has 
been associated with the onset of cross-slip by several authors, and the 
initial stage III stress, extrapolated to 0°K, Ty,;,, is above the twinning 
stress at 4:2°K (Seeger et al. 1959). The copper crystal which was typical 
of the crystals which twinned (Blewitt ef al. 1957, fig. (3)) was work- 
hardening linearly at the twin stress, indicating that cross-slip is not 
important. Also as for process (1), there appears to be no special reason 
why the dislocations AC, on planes (6) or (d), should lie along that 
particular direction which is also in plane (@). 

If process (3) occurs, then the long jog with Burgers vector AC is con- 
strained to lie along BC and the dislocation is pinned where it leaves the 
plane (0). This source should withstand the high stresses needed for 
twinning as envisaged by Cottrell and Bilby (1951). 


4.3. Operation of the Prismatic Twin Source 


The operation of the long jog in AC as a twin source will now be considered 
in more detail. If the jog is formed by the motion of the slip dislocations, 
Burgers vector BC or CB, on the same atomic plane, the jogged dislocation 
AC (or CA) has the configuration of fig. 6(a@). The dissociation of the 
jog AC is shown in fig. 6(6). The twinning dislocation «€ rotates in 
opposite directions around the poles X and Y (fig. 6(b)) until it meets the 
sessile Aa again. The partial aC approaches Aa between either the upper 
(A) or the lower (B) pair of atom planes which bound the core of the sessile 
dislocation, as shown in fig. 3(b). The actual position of the partial dis- 
location aC (or Ca) is between that pair of planes which is further away 
from the pole dislocation, and hence the configuration after one revolution 
may be that of fig. 6(c), if the twinning dislocation spirals round both 
poles X and Y. The partial aC which rotates round X and Y will meet 
one atomic plane apart along RS, as the pole dislocation has a component 
Burgers vector perpendicular to and equal to the spacing of the (111) 
planes (a). 

Seeger (1956) has suggested that a dislocation such as the twinning 
partial aC may acquire sufficient kinetic energy as it rotates around the 
poles X and Y to overcome the stress field of the configuration RS (fig. 6 (c)) 
dynamically. If this very large stress field (resolved shear stress 
o~180kg/mm?) is overcome after one revolution of the twinning dis- 
location, then the partials comprising. RS can rotate about the poles Y 
and X, generating a lenticular twin without further hindrance (see fig. 
6(f)). It is not, however, necessary that the twinning dislocation should 
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Fig. 6 
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The operation of a long jog in the dislocation AC as a twin source. 


(a) The dislocation AC which is jogged along CB. 

(6) Dissociation of the jog, with the Pantie dislocation aC rotating 
about the poles X and Y. 

(c) After one revolution of the twinning dislocation. A unit jog has been 
created at N, and is moving along the long jog, reforming the twin 
source Ty. 

(d) Operation of the twin source T,, producing a second layer of fault 
above the first. A unit jog has been created at N, and is moving 
towards R, producing the source T,, which will operate (see (f)) to 
continue the twin downwards. 

(e) After a complete revolution of the twinning dislocation from T,. 
Annihilation of twinning dislocations of opposite sign has taken place 
along RS, leaving a continuous helical layer of fault. 

(f) The form of the twin after several operations of the sources T, and 
T,. The twin is a continuous helical fault and has a biconvex 
lenticular shape. It is viewed along aD in this diagram. 


390 J. A. Venables on 


travel fast enough to overcome this obstacle dynamically for twinning to 
occur from such a source. If a helical twin a few (~10) atomic planes 
thick can be formed by repeated operation of the twin source, as suggested 
in $2, then the partials p, and p, (fig. 6 (e)) which bound the twin above 
and below are sufficiently far apart to overcome their mutual stress field , 
statically under the applied stress alone. Once a twin this size has 
been formed it can grow by the rotation of the partials p; and p, around 
Y and X respectively. The final width of the twin will be limited by the 
back stress on the rotating partials and by the length of the pole dislocation. 
The prismatic glide which is necessary for the first few repeated operations 

of the twin source may take place in one of several ways. The two 
mechanisms considered here are 


(1) The passage along the long jog of a unit jog created at the node N, 
and/or at the node N,. This mechanism will only be possible if 
the extended pole dislocation on plane (b) is constricted at the nodes. 


(2) Reactions taking place along the length of the long jog. These 
may take place if the nodes N, and N, are stable pinning points for 
the twinning dislocation aC. It is not at present certain that the 
nodes N, and N, are stable in the configuration of fig. 6 (c). 


If a pole dislocation (X or Y, fig. 6(c)) is constricted at N where it is 
jogged, a unit jog can be formed which is glissile in the same surface (the 
plane (d)) as the long jog, thus moving the twin source through the 
required distance AC. Unit jogs (of length AC) will tend to be formed 
for two reasons. They are expected to have a lower energy than multiple 
jogs when the pole dislocation is not exactly in screw orientation; and 
the twin source which is reformed by the movement of a jog can split on 
plane (a) with lowest energy when the new stacking fault behind «C is on 
the next atomic plane to the old one (see fig. 3(d)). Because the sessile 
Ac is split into the stair rod dislocation a along BC and Aé& on plane (d) 
the dislocation configuration at the two nodes N, and N, is not sym- 
metrical (see fig. 7 (a) and (c)). A unit jog created at N, will have a higher 
energy than one from N,. Any shear stress concentration on plane (a) 
will force a jog from N, along the long jog and a jog from N, back to N,. 
Hence the source T, (fig. 6 (d)) may be the only one which operates, but 
T, is shown operating in the sequence of fig. 6 for the sake of clarity. In 
any case, whether T, and/or T, operate as twin sources, more helical 
layers of fault, and a more extended composite dislocation line RS will be 
produced (fig. 6(e)). The component partials aC of the dislocation RS 
which lie on the same atomic plane can annihilate in pairs, leaving only 
ik top and bottom partials p, and p, and a single helical fault between 
them. 

Reactions taking place along the length of the long jog are shown in 
fig. 7. The configurations (a) and (b) refer to lengths of the long jog near 
N, and (c) and (d) near N,. Because there is a shear stress concentration 
on the (111) planes (see §4.4), the §C dislocation produced by the 
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coalescence of aC and §« moves towards the AS dislocation on plane (d). 
It is energetically favourable for the reverse reaction 


§C + $a+aC 


to take place on the next atomic plane because of the presence of the first 
fault. If the configuration of fig. 7 (b) and/or (d) occurs a twin is produced 
which is equivalent to the one obtained by generating unit jogs at N, only. 
It appears most likely that this is the twin produced from this source. 

In all cases considered, this prismatic glide of the twin source can take 
place at a lower stress than is required to split AC dislocation on plane 
(a) initially. Therefore the critical stage for twin nucleation is the 
creation of the first ring of fault from this source. The stresses required 
for twin nucleation are considered in the next section. 
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Prismatic glide of the twin source AC. The lines represent intrinsic stacking 
faults, and the letters A and B give the positions of the faults with respect 
to fig. 3, which is in the same orientation. The reaction (a) > (6) occurs 
near the node N, and (c) — (d) near N,. See text for discussion. 


4.4. The Stress Required to Operate a Prismatic Twin Source 


The resolved shear stress on plane (a) required to operate this prismatic 
twin source AC may be given by the following equation: 


2 
o.by=y 4s aT Mee yA Are RAST) 
where J is the length of the source and b,=aC. This equation assumes 
that apb,2/l is greater than the short range force needed to separate the 
aC and Aa dislocations. For a copper. crystal with a [111] tensile axis, 
the stress required to operate this prismatic source will be less than that 
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for either slip source (table 1) if J is greater than about 50 A. Thus 20 
primary slip dislocations on one atomic plane intersecting a secondary 
slip dislocation produce a jog of sufficient length to make it the most 
favourable type of twin source in the work-hardened crystal. Since a 
stress equal to twice the observed twinning stress is needed to overcome 
the surface tension of the fault alone, some shear stress concentration on 
the twin plane is necessary for twinning to occur at all. The stress con- 
centration factors required for the operation of this source range from about 
five for a 504 source to below three for a source which is longer than 200A. 
These shear stress concentrations could easily be produced in the neigh- 
bourhood of the jog by a small pile-up of primary slip dislocations which 
have not penetrated the twin source. 


§ 5. CORRELATION OF THE PROPOSED MODEL WITH 
EXPERIMENTS ON WORK-HARDENED CRYSTALS 


The most important result of considering the various types of twin 
sources is that the observed twin systems and the tensile axis orientation 
relationship of Blewitt et al.(1957) and of Suzuki and Barrett (1958) may 
be explained, if prismatic sources operate. Only the [112] twin direction 
(aC) is observed experimentally when twins are produced on the primary 
slip plane (111). Twins with a [211] shear are not produced in tension 
tests and this is explained, by considering the stress acting on the aB twin 
source, and by assuming that extrinsic stacking faults have a high energy. 
Both ‘ glide’ and prismatic sources will account for this distinction between 
the [112] and [211] directions of twinning. Only the prismatic source 
AC will account satisfactorily, however, for the dependence of twinning 
on tensile axis orientation. 

To satisfy eqn. (4.1) the resolved shear stress must be high along [112] 
and there must be some stress concentration on the primary slip plane. 
A high stress concentration can only be maintained if cross-slip is un- 
favourable and this is the case when the tensile axis is around [211]. (For 
a [211] tensile axis there is no stress on the cross-slip plane.) Therefore 
the tensile axes of the crystals which twin must be concentrated near 
[111] and away from [100]. The range of tensile axes for which twinning 
occurs must also decrease as the temperature is raised; i.e. the tensile 
axes must cluster around [111] as cross-slip becomes more prevalent. 
This is precisely what was observed by Blewitt et al. (1957). Ifthe internal 
stress is considerably different from the external stress, the scatter in the 
values of the measured twin stress obtained by Blewitt et al. (1957) and 
by Suzuki and Barrett (1958) is easily accounted for. 

In a copper sample (169) whose initial tensile axis was near the (100)- 
(111) zone (fig. 1), Blewitt e¢ al. (1957) found twins on the two conjugate 
systems. Conjugate twinning was also a general feature of the twinning 
in the silver—gold alloys investigated by Suzuki and Barrett.(1958). This 
may be simply explained on the proposed model. Just as jogs produced 
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by primary slip dislocations BC in AC dislocations are sources of twins on 
the primary twin plane (a), so the jogs produced in BC by the secondary 
slip dislocations AC are sources of twins on the secondary twin plane (0). 
Thus in conditions when the internal stress is the same on both systems or 
when there is local overshoot (Suzuki and Barrett 1958), conjugate twins 
will be formed. 

The high shear stress for continued twinning may be explained if new 
sources are necessary to propagate the twin throughout the crystal. 
This is consistent with the width of each twin being limited to the length 
of the pole dislocation. If eqn. (4.1) is an adequate expression for the 
twin stress, the twin stress should not be strongly temperature dependent. 
This was observed over the temperature range 20-273°K by Suzuki and 
Barrett (1958) and between 4-2 and 77°K by Blewitt et al. (1957). 

The prismatic source model will also account for the twinning behaviour 
of neutron-irradiated copper (Blewitt et al. 1957). The production of 
mechanical twins in neutron-irradiated f.c.c. metals is considered in the 
next section. 


§6. THe Errect oF NEUTRON-IRRADIATION ON TWINNING 


Recently it has been discovered that neutron-irradiated copper contains 
prismatic dislocation loops (Silcox and Hirsch 1959) similar to those 
observed in quenched metals by the same authors. A summary of their 
results and the results of Blewitt et al. (1957) on the effect of neutron 
irradiation on the twinning properties of copper single crystals is given in 
table 2. 


Table 2. Twinning in neutron-irradiated copper 


From Silcox and Hirsch (1959) 


Trradiation temperature (°c) 30 30 80 
Neutron flux (nvt) 6-7 x 1017 5-6 x 1018 1-4 x 102° 
Average loop size (A) 75 150 300 


From Blewitt et al. (1957) 


Irradiation temperature (°C) 30 80 
Neutron flux (nvt) jp i AU) 125.107" 
Crystals investigated 396, 419 416 
Twinning characteristics at 
77°K Twinned ina similar | Did not twin as 
manner to o<13kg/mm? 
work-hardened 
crystals 
Twinning characteristics Not reported Twinned after small 
4:2°K extension on secon- 


dary slip plane 
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When the irradiated crystals are strained, dislocations on the primary 
slip system (Burgers vector BC on plane (@)) will interact with the 
prismatic loops. Interaction of the glide dislocations with large loops 
may cause the loops to act as twin sources. In particular the annihilation 
reaction 

BC —loop BC=0 
on (a) on (bd) 


loop BC-+BB+BC 
on (b) on (5) 


The twin source in neutron-irradiated metals. 


(a) A dislocation of the primary slip system BC on plane (a) held up by a 
prismatic loop BC on plane (0). 


() The slip dislocation intersects the loop by moving to the left and 
displaces the two halves of the loop by a Burgers vector BC. The 
slip dislocation is now continuous with the open half of the loop. 


(c) When the stress on (b) is high enough the dissociation 
BC-- BB +8 C occurs. 


The twinning dislocation BC spirals round the pole X and creates a 
jog at N,. The new twin source T, is formed by movements of this jog 
around the loop, and it can dissociate again producing a second layer 
of twin. The pole Y may be inoperative because of the presence of the 


sessile loop BB. The growth of the twin is analogous to that described 
in fig. 6. 
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produces a source of secondary twins, as pictures in the sequence of fig. 8. 
Secondary twins were observed by Blewitt e¢ al. (1957) in erystal 416 at 
4-2°K, 

If it is assumed, following Silcox and Hirsch (1959) that the prismatic 
loops are formed by the condensation of vacancies}, then the opening in 
the loop produced by the annihilation reaction enables the intrinsic 
stacking fault behind BC to expand freely (fig. 8(c)). After the loop has 
been cut, and the fault expanded on plane (b), the twin can grow by 
prismatic glide of the loop, in a similar manner to the twins in the work- 
hardened crystals. If unit jogs are created at the nodes where the source 
dislocation changes from plane (b) to plane (a), the mechanisms for work- 
hardened and neutron-irradiated crystals are anologous. 

If the stress concentration for twinning on the secondary system is small, 
the twin will only nucleate from large loops. In the case of crystal 416, 
the tensile axis had not rotated onto the (311)-(111) zone when twinning 
occurred and so the stress concentration on the secondary slip plane cannot 
have been large. Thus it is not to be expected that crystals containing 
small loops twin on the secondary slip plane. The crystals 396 and 419 
(table 2) were observed to have similar twinning properties to the work- 
hardened crystals at 77°K, but their behaviour at 4-2°K has not yet been 
reported. 


§ 7. CONCLUSIONS 
The main conclusions of this paper are as follows: 


1. Prismatic dislocations are the most likely source of mechanical 
twins in work-hardened and neutron-irradiated metals. These sources 
can operate by a modification of the mechanism of Cottrell and Bilby 
(1951). 

2. The twin systems observed and the tensile axes necessary for twinning 
to occur in the experiments of Blewitt et al. (1957) and Suzuki and Barrett 
(1958) can be explained in terms of the operation of the ‘prismatic’ 
sources. The stress concentration factor required for their operation can 
be as low as 2 or 3, given a source length of about 200A. 

3. Glide dislocations have been considered as twin sources. It is 
considered that the experimental results cannot be explained in terms of 
the operation of such sources. 

4. It is thought that suitable prismatic dislocations are produced in 
the work-hardened material by the jogging of the dislocations of the two 
conjugate slip systems. A long jog (250A) in a secondary slip dislocation 
is a source of twins on the primary slip plane. 


+ Note added in proof.—Since this paper was written, Barnes and Mazey 
(1960) have shown that the dislocation loops observed by Silcox and Hirsch 
(1959) are probably produced by the clustering of interstitial atoms. The 
model described here is also applicable to the case of interstitial loops, the 
fault (fig. 8 (c)) nucleating outside instead. of inside the loop. 
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5. In the neutron-irradiated copper crystals, prismatic dislocation 
loops produced by the collapse of discs of vacancies are observed to be 
present (Silcox and Hirsch 1959). Certain prismatic loops on the secondary 
slip plane can act as sources for twinning on this plane when they are inter- 
sected by the primary slip dislocations. Twins were observed on the 
secondary slip plane by Blewitt ef al. (1957) in some neutron-irradiated 
copper crystals strained at 4-2°K. 
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ABSTRACT 


Tron—0-5% carbon was fatigued in air and in vacuum. The knee in the 
S-N curve at the fatigue limit occurred at about the same number of cycles 
in vacuum as in air but at a higher strain. Above this strain specimens 
lasted about ten times as long in vacuum as in air. Slip occurred in a 
similar manner above and below the fatigue limit and cracks formed in the 
slip bands in both cases. Below the fatigue limit the cracks grew a few 
microns long. 


§ 1. INTRODUCTION 


Previous work (Gough and Sopwith 1932, 1935 and 1946, Wadsworth 
and Hutchings 1958, Wadsworth 1959 a) has shown that the presence of 
air can have a considerable effect on the fatigue life of various metals. The 
air increases the rate of growth of the cracks, apparently by some chemical 
action. All the experiments so far have been on metals without a fatigue 
limit in which presumably failure would occur eventually at all stresses. 
Steel, however, is peculiar in that it has a fatigue limit—that is a stress 
above which specimens fail in a finite number of cycles (usually < 107) and 
below which they apparently last forever. There is inevitably some scatter 
in the behaviour of different, nominally identical, specimens but it is clear 
that there is some real physical discontinuity in behaviour at the fatigue 
limit (see for example Clayton-Cave et al. 1955). In view of the close con- 
nection between strain-ageing phenomena, yield points and fatigue limits 
one might expect that the fatigue limit would be the stress below which slip 
could not occur. Hempel (1956) has shown that thisis not so. The discon- 
tinuity must occur later in the series slipcrack formation—failure. 
From the previous work, quoted above, one would expect that the presence 
-of air would speed the spreading of cracks once they formed, but have little 
effect at stresses where cracks never form. Thus a comparison of the fatigue 
limit in air with that in vacuum would provide further evidence about the 
nature of the fatigue limit. 


+ Communicated by the Author. 
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This paper describes the results of some experiments allowing such a 
comparison, together with the results of metallographic examination of 
specimens stressed above and below the fatigue limit. 


§ 2. MATERIAL 


The tests were made on a high purity iron with about 0-5% carbon added. 
The carbon was distributed uniformly throughout the specimen. Figure 1 
shows a typical microstructure. 


§ 3. APPARATUS AND METHOD OF TEST 
The specimens were fatigued in reversed bending between constant 
strain limits in the apparatus previously described (Wadsworth 1959b). 
The frequency was 100 c/s and the amplitude was built up to its maximum 
value in a few thousand cycles. The tests were run to completion without 
interruption unless otherwise stated. 


§ 4. RESULTS 

The specimens were fatigued at various strains in laboratory air or in 
vacuum. ‘The results are shown in fig. 2. In the higher strain range the 
life in vacuum was about ten times that in air. At lower strains there 
appeared to be a fairly well-defined fatigue limit both in air and in vacuum, 
The two limits occurred at approximately the same number of cycles but 
at different strains. There was a range of strain in which specimens would 
break in air but not in vacuum. 

The specimens were examined microscopically at the end of the test. 
The amount of slip visible was smaller on specimens subjected to smaller 
strains but some slip was visible on allspecimens. There was no qualitative 
difference between the appearance of specimens fatigued just above and 
just below the fatigue limit. Etching specimens after fatiguing showed 
that the slip occurred in the ferrite grains and not in the pearlite, as shown 
in fig. 3. 

Some specimens (not shown on fig. 2) were examined at intervals 
during the test. As any one test progressed the slip became more 
pronounced and obvious cracks appeared towards the end of the test. How- 
ever, cracks were present in some slip bands earlier in the test and could be 
detected by extending the specimen thus causing the cracks to open. A 
number of specimens fatigued at a strain of 1-34 x 10-8 were examined. 
The life in vacuum was about 10’ cycles and in air about 10% cycles. Cracks 
could not be detected in specimens fatigued for 10° cycles either in air or 
in vacuum but there were many present after 10° cycles in vacuum.’ 
Typically the larger of these cracks were 50-100 long and 5p deep. 

Specimens fatigued for between 5 x 107 and 7x10? cycles below the 
fatigue limit in air or in vacuum were examined after extension. In all 


} Figures 1, 3, 4 and 5 are shown as plates. 
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cases small cracks were found in a few slip bands. (A control test on an 
unfatigued specimen showed no similar cracks.) Those in the specimens 
fatigued in vacuum just below the fatigue limit were in general larger than 
those in specimens fatigued in air at a ELD stress. Figures 4 
and 5 show an example of each. 


Fig. 2 


2:0x 10> 


STRAIN 


15x 10° 


6 107 0° 


GY GEES 


3 
10x 10 
10% 


S-N Diagram. ¢testin air; Otestin vacuum; Atest in vacuum after fatiguing 
in air; [) test in air after fatiguing in vacuum. 


One specimen was fatigued for 7-5 x 10’ cycles at a strain of 1-26 x 10-8 
in vacuum and then the test was continued in air. The specimen was still 
unbroken after a further 5-6 x 10’ cycles although the life of fresh speci- 
mens strained at this level in air is only about 10° cycles. Two other 
specimens were fatigued for about 10% and 15% of their expected life in 
air at a strain of 1-18 x 10-3 and the test continued in vacuum. ‘The first 
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was unbroken after 3°3x107 cycles and the second survived until 
3 x 107 cycles when it broke possibly due to an equipment failure. 


§ 5, Discussion 


The results shown in fig. 2 show that the lives of specimens fatigued at 
high strains are about ten times longer in vacuum than in air. This factor 
is similar to those found for copper and aluminium (Wadsworth and 
Hutchings 1958). The interesting point, however, is that the fatigue 
limit occurred at approximately the same number of cycles in air as in 
vacuum but at a different strain. Thus there is a range of strains in which 
the specimens break after a few million cycles in air but seemingly last for 
ever in vacuum. It appears that fatigue cracks form in the slip bands at 
all the strains used in these experiments but that they can become non- 
propagating if the conditions are such as to give sufficiently slow initial 
growth. These conditions may be either a low strain (1-05 x 10~*) in the 
presence of air or a higher strain (1-3 x 10~*) without air. Once the crack 
has stopped propagating it is difficult to start again. Thus specimens 
fatigued below their fatigue limit in vacuum have small cracks in them but 
these do not propagate if the test is continued in air although fresh speci- 
mens would break at this strain. Presumably this difficulty in restarting 
a crack (or starting a new one) is also responsible for ‘ coaxing’ (the increase 
of fatigue strength produced by prefatiguing below the fatigue limit) 
(Sinclair 1952). Coaxing and the fatigue limit appear to occur only 
together and with strain ageing. 

It is interesting to consider if it is necessary for there to be continuing 
reversed plastic flow at the tip of a non-propagating crack or whether it is 
possible for the deformation at the tip to be elastic. 

The theoretical elastic stress concentration factor at the tip of a crack 
of depth d (small compared with the specimen diameter) and tip radius p 
is K,=1+2+/(d/p) (Neuber 1946). 

The maximum strain occurs over a very small region at the tip of the 
crack which may be very work hardened. Even so an elastic strain there 
of more than | or 2% is unlikely so K, cannot be more than 10 or 20 on 
these assumptions. For cracks 3. deep these correspond to root radii of 
about 15004 and 300A respectively which are not impossible. It is very 
doubtful however if the strains can be wholly elastic in the case of larger 
non-propagating cracks such as those formed at the root of sharp notches 
in mild steel. Frost (1959) summarizes the evidence on these and shows 
photographs of cracks whose depth including the notches, is 13mm and 
whose tip radius is less than 0-4. This indicates K,>100. These 
cracks did not spread under a nominal stress of YM/2000 which would 
produce a strain of 5% at the crack tip. This is improbably high and it 
seems almost certain that reversed plastic flow must have occurred in each _ 
cycle without the crack spreading; it is quite possible that the same 
occurred in the tests described here. 
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If we assume that the crack was sufficiently sharp for plastic flow to be 
necessary it is of interest to calculate the order of magnitude of the flow 
involved. The depth of the plastic region at the bottom on the crack 
depends on the depth of the crack, the applied stress and the yield stress of 
the work hardened metal under the triaxial stress conditions at the tip of 
the crack. Lack of knowledge of the last factor precludes exact calculation 
but clearly the depth of the plastic region will be of the same order as the 
crack depth. The local plastic strain in the zone varies from zero at the 
inner edge to a large value in a small region at the tip of the crack. How- 
ever, the whole region is surrounded by elastic material and the average 
plastic strain in it must be comparable with the elastic strain. Thus the 
total flow in the plastic zone is of the order of the crack depth multiplied 
by the maximum elastic strain. This is of the order of 1004. In other 
words in these specimens stressed just below their fatigue limit there was 
probably a region at the tip of each crack where some tens of dislocations 
move distances of perhaps half a micron each half cycle without producing 
further spreading of the crack. 

It seems from this work that the fatigue limit is the strain below which 
cracks cannot propagate rather than the stress at which they cannot form, 
but more work is needed before the mechanism by which they are stopped 
is clear. 
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ABSTRACT 


A perturbation solution of the differential equation controlling the 
migration of impurity atoms to a dislocation is given and boundary conditions 
appropriate to precipitation on the dislocation are derived. Numerical 
results within the range of validity of the perturbation analysis are obtained 
and used to discuss the ageing of carbon and nitrogen in iron. 


§ 1. INTRODUCTION 


Ir is now well established that there is a strong interaction between 
impurity atoms and dislocations. We envisage two distinct situations 
at the dislocation cores ; either (a) some of the impurities there are rendered 
immobile and effectively removed from solid solution, that is, local pre- 
cipitation is occurring on the dislocations (Bullough, Newman and 
Wakefield 1959, Bullough, Newman, Wakefield and Willis 1959, 1960, Ham 
1959); or (b) the impurities all remain mobile and the equilibrium state is 
a Maxwellian concentration distribution (Bullough and Newman 1959): 
this will be referred to as an impurity atmosphere. Situation (a) could 
arise if either the local impurity concentration exceeds the equilibrium 
solid solubility or two impurities with a strong chemical affinity for each 
other are present in the core, since this provides an ideal site for nucleation 
of a second phase. (6) will arise if the impurities are chemically inert 
and if the final concentration in the core region is still sufficiently small 
‘for the probability of occupation of sites to be controlled by Boltzmann 
statistics; it is also necessary that this final concentration should be less 
than the equilibrium solid solubility of the impurity. 
In the present treatment we shall consider a crystal initially containing 
a uniform concentration of impurities and investigate the kinetics of their 
migration. This is controlled by two processes, namely, diffusion and 
drift flow (Ham 1959, Bullough and Newman 1959). The solutions 
obtained of the differential equation take full account of the diffusion flow 
while a perturbation technique has been used to obtain the contribution 
of the drift flow when this is small. Boundary conditions at the dis- 
location cores are taken appropriate to either (a) or (6), and allowance is 
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made for a finite dislocation density. Also in the treatment of (a) due 
allowance is made for the possibility of a rate-limiting process occurring 
at the precipitate-matrix interface. This treatment therefore is essentially 
an extension of the previous theoretical analyses (Ham 1959, Bullough 
and Newman 1959, Cottrell and Bilby 1949, Harper 1951, Bilby 1955), 
each of which ignored at least one of the above factors. The general 
solution of the governing partial differential equation and a detailed 
formulation of the relevant boundary conditions are given in §§2 and 3 
respectively. The required solutions satisfying the precipitation boundary 
conditions are derived in §4 where the atmosphere solutions appear as a 
special limiting case. The solutions have been computed for various 
values of the parameters involved and are illustrated graphically in §5. 
Finally, the theory is discussed with reference to some experiments of 
Thomas and Leak (1955) on the ageing of steel. 


§2. Tue GENERAL SOLUTION OF THE DIFFERENTIAL EQUATION 
GOVERNING THE MIGRATION OF IMPURITIES NEAR A DISLOCATION 


The governing differential equation for the flow of impurities in the 
neighbourhood of a dislocation, where the interaction potential has the 
form 

d¢(r)= —A/r Pe Peer a ee 
has been shown to be (Bullough and Newman 1959) 


1 Oe(r,t)  Oc(r,t) | Loc(r,t) _ L[1dc(r,t) 1 
Di on ee he eee +Z[5 or | (2) 


where L=A/kT. Attempts to obtain an exact closed solution of (2) 
have been unsuccessful and recourse has been made to the following 
approximate analytic technique. The terms in (2) not involving Z 
constitute the well-known radial diffusion equation for which the general 
solution is known. We therefore look for a solution consisting of the 
diffusion solution plus additional terms of order Z/r. Such a solution 
will, of course, only be valid if L/r<0-2. For carbon in iron the value of 
L is about 404 (Cottrell 1953); however, much smaller values may be 
expected when the ionic radius of the impurity is similar to that of the 
atoms forming the parent lattice and where these impurities occupy 
substitutional sites. It is reasonable to associate the minimum value of 
r with the radius of the region at the centre of the dislocation inside which 
linear elasticity is not valid; this radius is probably less than about 10 A. 
Thus, it is clear that our first order solution} in L/r will not be directly 
applicable to the case of carbon in iron. But nevertheless it should 
indicate the relative importance of the diffusion and drift flows at various 
stages in the ageing process in iron and may well be directly applicable to 
other systems. ; 

RTE een aa eee et 


+ The second-order perturbation has been carried through and it is found 
that the coefficient of (L/r)? cannot be expressed in closed fara ee 
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Transforming eqn. (2) by means of a Laplace transform, and writing 


Clr A) = | c(r, t)exp(—qt)dt (where g=Dh)?) 
0 
we obtain the ordinary Ca aie Seem for ¢(r, A): 


Pe(r,A ] a 


7 r2 


where Cy is the initial uniform concentration and D is the diffusion coefficient. 
The general solution of (3), correct to order L/r, must have the form: 


e(r, A) = A(A)K (Ar) [1 + Lg(r, A)]+ BA)L(Ar) [1 + Lh(r, A)] 


+ 5 + Lf(r, d) big Geach Wa tins 2. (4) 


where A and B are arbitrary functions of X to be determined from the 
boundary conditions; (4), with L=0, is the well-known solution of 
Bessel’s equation to which (3) reduces when L=0. Substituting (4) into 
(3) leads to three elementary second-order differential equations from 
which the functions g; h and f can easily be deduced: 


CF ete ee ” Ko(At) Ky (At) 
ar N= 5m {La0r) | i eu ds 


| + K,(Ar) | eee ae (5) 
y r “? Ty(At)L,(At) 
rr eae as Kole) i aa 
+ 1,00) |’ aa, mah pet ede) 


fr, N= 7 +21, Qn) ik K,(xt)dt — K,(ar) (he "Ti0) ae ee 


The arbitrary constants associated with the subsidiary differential equations 
for g and h have been compounded into the A and B of the complete 
solution (4) and it can be shown that 


1K, (Ar) 

A eae E07) (8) 
recy) 

hr, d)= — + T.0w) ; | (9) 

flr, N= phd + Fetal) 4001} ee a1) 


where L, is a Struve function of imaginary argument (see Watson 1952, 
p. 329). However, because of the complexity of the boundary conditions: 
it was convenient to leave g, h and f in the integral form (5)—(7) until after 
inversion of the Laplace transform. 
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§3. THz BounDAaRy CONDITIONS 


The model of the core is similar to that used in our previous treatment 
for the growth of an impurity atmosphere (Bullough and Newman 1959). 
We surround the dislocation line by a virtual cylindrical surface of radius 
7) inside which the impurity atoms have very high mobility (in this region 
we designate an infinite value for D). Observations indicate that the 
precipitates nucleate and grow at discrete points along the dislocation. 
The basic feature of the model is that the rate of growth of these particles 
is assumed to be proportional to the concentration of impurity in the core 
region. The constant of proportionality will be varied between infinity, 
the fastest possible rate of precipitation and zero, corresponding to the 
formation of an atmosphere with no precipitation. 

In the derivation of this boundary condition, we shall neglect any back 
flow from the precipitate particles; such a flow would lead to a non-zero 
equilibrium distribution of impurity in the crystal. If, however, the 
equilibrium surface concentration at the precipitate matrix interface is 
small compared with the initial concentration c, then omission of this 
term is not important. We write the core boundary condition : 


C(1'9, t) =c,(t) 5) tie Sok Sh, hake poe ee 
where c,(¢) is the concentration of mobile impurities within the region 
r<Y. C,(¢) may be determined since its rate of change is simply the 
difference between the net flow of impurities into the core region F'(7, f) 
and their simultaneous removal as a result of the precipitation; thus 


deg(t) _. 
“ew A = ath ryt} —arge (Oh so) eae oo Wk 
where P is the rate-limiting constant governing precipitation and 
ZL 0 
F(r, t)=D | alr t+ X(r,, | 4 ie bide Us 
ae or 


Integrating (12) with the initial condition 
C,(0) =¢y re Meas Tes een) 


we obtain, with the aid of (11), the general core boundary condition: 
2 t 
C(1'o, t) =cyexp [— Pt]+ ~ exp [— Pt] | exp [Pt]F (ro, t) dt. (15) 
0 0 ; 
If P=0, (15) becomes: 
2 rt 
C(ro, t) =Cy+ =| F (ro, t) dt. <M de So EL) 
Yolo ; 


which is the core boundary condition for the growth of an impurity atmo- 
sphere. Similarly, if P+ oo, (15) becomes: | 
G(rg, A= 0 | 6 go cal eg tara Tan 


which is the core boundary condition (Ham 1959) appropriate to pre- 
cipitation without rate limitation at the precipitate matrix interface. 
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In a real crystal there will be a finite dislocation density p and we may 
associate with each dislocation a cylindrical surface of radius R= (mp)-"?2. 
The second boundary condition is simply that there should be zero net 
flow across this surface; thus 


Ho(R, + Zep, t)=0. ee are tes (18) 
For an isolated dislocation, p>0, (18) leads to the boundary condition 
c(r, t) cy, when r> o. 

The quantity eventually required is the number of impurities W(t) 
precipitated in time ¢ per unit length of a dislocation, which is given by 


W(\=nreP | gt) aw Sd hd 3) 
0 


where c,(¢) is defined by (11) and (15). It is not convenient to obtain W 
directly, and we first evaluate the number of impurities V(t) which have 
entered unit length of the core region r<r,in time ¢. We have, from (13), 


Lal ee Fe dc 
V(t) =2nryD | F E C(%o, t) + >, (79; | Cie tetas vet 20) 


This quantity rather than (19) is required for the case of the growth of an 
impurity atmosphere (Bullough and Newman 1959). By differentiating 
eqn. (20) and making use of (13), the relationship between W and V 
follows: 


¢ he Bs dV (x) 

W (t)=7192¢)[1 — exp (—Pt)] ail exp (— Pt’) dt | exp (P2) a cas 
0 ) 

(21) 

The Laplace transforms of the boundary conditions and the quantity 


V(t) are: from (15) and (13): 


Fenn af a”) (rooNe at a) ie 


where the dimensionless parameter a” is given by a2 =17,)2P/D; (18) becomes: 


2 2 Lie dé 
Aro N= ee + at | alr N+ Flr) | 22) 


Stim (ea Oue etme ta 2 2 (28) 
and ; 
Pay= =] are n+ Sire »| Beno) 


§ 4, THe REQUIRED SOLUTIONS 


The two arbitrary functions A(A) and B(A) in the general solution (4), 
with g, h and f given by (5), (6) and (7) respectively are obtained by satis- 
fying the boundary conditions (22) and (23). The resulting expressionst 
ee he Ae a EE EEE 

+ In the analysis L/r, is treated as a small quantity and only terms of this and 
zero order are conserved. 
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for G(r, X) and dé(r,, X)/dr, are then substituted into eqn. (24) for V(A). 
After considerable algebraic manipulation it can be shown that: 


V(A) “ a4 (79) 
Zmrgcg DA [a Wyo (79) +A*F9'Wr0("0) | 
L (792A? + ax?) | sis 3 Ai 
To Ds [221 9(7) a 7920 9(7'9) ig is Wy(Y) [a (W49(7o) W49(Y)) 
+)2r,2wio(ro) |ayly> = ee) 


where w;;(x) =I,(RA)K,(xA) + (— 1) K,(RA)L, (wa). 
The quantity V(t) is now given by the inverse Laplace transformation : 
c+io _ 
Vit)= Vy exp giidqen - lee anton) 
a7 J c—ico 
where A=1/(q/D), and the Bromwich contour appropriate for a single 
valued function of g is used; an unusual feature of this inverse Laplace 
transformation is the evaluation of the residues at the double poles of 
V(A). The integration over y is then performed; we obtain after con- 
siderable reduction : 


V(t) Be Aa? S exp (—A,,?2M)v,,(,,) ny 4k OL 
ies) oem Aci <0 oS oe 
© exp(—A,2M 
x SPMD fatat — KO yl + 08D, — (0? — KY) AB] 


(27) 


where V(0)=7¢)(k?—1,?), M=Di/R?, K=r,/R and X,, are the positive 
roots of 


Y (A) [a?J (KA) + 22. 4(KA)]— J4(A) [a? ¥9( KA) + K2A2Y,(KA)]=0. (28) 
The quantities A,, B,, C, and D, are given by: 
A, = (a? — KA,,7)(V99 — K04,) + 2KA,0 91, 


2 


9 Kx | 
Ve - re + H,(A,) — KA, +0 - Mies. [y.H (KA,,) — 9H, (KA,)], 


8 
Ga sats 2A LEA (Vir? + 107) — P1110]; 


Sav? 
D,=- AS + oP K [80417 + 8049? — 2An (0911 + Yo0%10) ] 
v1); 


— 07 [01029 + Vy9V01] + K?A,,2 [49% — Voo11} 
—= BA, 2 [Via M49 245 Po1hay asta) ne ee (29) 
where vy = Viz(An) = Y(A,)J;(KA,) —JA,) Y ;(KA,). 1 St SO) 


W(t), the number of impurities precipitated in time ¢ per unit length 
of the dislocation, can now be obtained by differentiating (27) term by 
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term with respect to ¢ and then substituting the result in (21), and per- 
forming the elementary repeated (term by term) integration. If we write 
V(é)/V( co) in the symbolic form: 


(t)/V(coo)=1+4+ SS X,, exp (— OWES 1A Magee — SS (Y,M+Z,,) exp (—A,?2M) 


n=0 "QO n=0 
(31) 
where the values of X,, Y,, and Z,, follow from (27) we obtain 
W (t) : a? — AE 
aoe te 9 AS — 27 ot [Ke WARS EE cat ee 
Wa) exp (—a?m)+ red K ) 2 Ke, ae d?,M) 
—exp(—am)}+ = (1K?) 
x { s [exp ( ae x") — exp ( tn MILY See ee =n) 
n=0 tae] KA 7) 
(a?/K*—d,?) YM exp (—A,?2M) 
u (c2/K2—n,2)2 } sed 


where W(cc)=7R*c, and m= Dt/r,?. 

If «+0 in (27), that is precipitation is inhibited and an impurity 
atmosphere forms round the dislocation, then V(t) N(t) given by eqn. 
(49) of our previous atmosphere calculation (Bullough and Newman 1959). 
It should be noted that when taking the limit as « tends to zero in (27) the 
first root Ay of (28) also tends to zero; in fact Aj>a—>0. In the present 
notation we have for the number W(t) of impurities having entered the core 
region r<7, in time t: ' 


N(t) 
N( 2) 


=1— a, 2, Nel i eee (93) 
ee 0 
where 


eeu. ie ie ee DRED - +3 (22,00) — Kd Kq)) | 


§,, = KA,?[V29—-Kv,,] and A, are the positive roots of 
JAN Y (KAY— JS (KAN Y{(ANH00 0. oe (84) 
N(co) in eqn. (32) is given by 
N(0)=7Lrye,(1—K)?. Ne ee Nee (OD) 


If the dislocation density is sufficiently low, any interaction between 
adjacent dislocations will be negligible and the kinetics of precipitation 
and atmosphere formation on such dislocations may be obtained by 
considering an isolated dislocation in an infinite medium. The he aes 
expression for (A) is again obtained by substituting ¢(r9, A) and dé(rq, A)/dr 
into (24). In this case, however, the arbitrary function B(A) is set equal 
to zero; that is there is only one arbitrary function A(A) in what is now 


410 R. Bullough and R. C. Newman on the 


the bounded general solution. A(A) is determined from the core boundary 
condition (22). The final result is: 


VA) ao Ky (r A) P(A) _ LZ (792A? +. 0?) 


Ne eee ss Eee 
2a Co D* a yy One DAS ) 
x {aro {" sO, E an ay | 
(36) 
where 
(A) = [02K (1A) +17 97A2K o(7 9A) |. os us et al) 


V(t) is given by the inverse Laplace transform using a ‘cut’ Bromwich 
contour appropriate for a multi-valued function, the only residue being 
at AX=0. After performing the integration over y we obtain for the 
number of impurities V(t) which have entered unit length of the core in 
time ¢: 


Vit) _ 20? (> [1—exp (—mz?)] { a2 — x? d 
47.2%) oT fF a A,?(x) + A, a 3 
+ =|. [1 exp (=m?) (e?—08) 04] Hy) — Y,(x)— ek 

0 7x 


19 


x*[ A,?(x) + A,?(x) ] 
my abe 2>[A,2(@) + A2(@) si 


(38) 
where m=Dt/rq?, A,(x)=07Yo(x)+2?Y,(x) and A,(%) =a?J9(x) +22J,(2). 
The quantities R(x) and S(x) are given by: 

R(x) = 2a[ Jy (x) Yy(x) + So(%) Yo(x)] — [J o(x) V(x) + Yo(x)Jx(x)], 
S(x) =a T(x) + Fo%(e) — Yy%(x) — Y¥q2(x)] — Jo(ux)Jy(x) Ge 
+ Yo(x) Yy(2). 
W(t) can be obtained directly from (38) by following the procedure given 
above. 
a can be set equal to zero in (38) and we get 
Vi) 4 (eee ee Y(e 


TY Co L Cu 


)=(2/ray1+22/3}} 4, 
0 a Jy?(x) + YQ?(x) 

+ (40) 
for the growth of an impurity atmosphere round an isolated dislocation. 


{In our previous treatment for the growth of an impurit 

y atmosphere 
(Bullough and Newman 1959) the term J,(x) in the numerator of the see 
(eqn. (43)) should be Y,(x) as shown in (40) above. The correct term was, 


hk included for the numerical computation described in this previous 
work. 
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§ 5. Numericat REesuurs 


The expressions (27), (32) and (38) have been computed for various 
values of the parameters. In order to compute (27) and (32), we have 
first to determine the positive roots of (28). At this stage we must allocate 
numerical values to the parameters ao? and K=r,/R. We choose. an 
arbitrarily large value of 604 for r, which for a dislocation density of 
3 x 10" lines/sq.cm leads to the value of K of 0-1838; by a reduction of 
79, the same K would, of course beapplicable toa greater dislocation density. 
The values of «? taken were: 


0-03, 0-1, 0-3, 1-0, 3-0, 10-0 and 1000, 


the final value corresponding approximately to no rate limitation (a? = 0) 
at the precipitate matrix interface. The roots of (28), for the above 
values of K and a”, were obtained by an iteration procedure involving the 
use of a high speed digital computer. The first four roots were found for 
each «? value and are given in the table. These roots were used to sum 
the first four terms of (27) and (32) for the range of M =2" where n runs 


First four roots of expression (28) with K =0-1838 


Ay As 3 


0-172100 3°871807 7-196025 10-675518 
0-308740 3:914209 7:221338 10-688623 
0-511182 4-028760 7:293896 10-727002 
0-815088 4-353662 7546533 10870264 
1-082666 _ 4-861475 8-161279 11-331201 
1-258545 5308838 9:005615 12-525635 
1-359912 5:569385 9-491260 13-372607 


from —7 to +9 in unit steps. The constant factor between successive 
values of M was taken to provide a uniform spacing of points when JM is 
plotted on a logarithmic scale. Clearly, for M>1, only the first term in 
(27) or (32) is required while for M <0-1, more terms in the series will be 
necessary. The error introduced by restricting the series to four terms 
was determined by computing (38) since, by suitable normalization (see 
ordinate scale of fig. 3), this solution for an infinite body becomes identical 
to the early stage solution (small values of J) for a finite body; for 


1>M>0-01, 


the error was, in fact, negligible. 

Figure 1 shows the variation of the precipitated fraction with time 
(expression (32)), in the absence of any drift flow (L=0). The expression 
has also been computed assuming L/r,=0-2, which is a value consistent 
with the first-order perturbation solution. On the scale used in fig. 1, 
it'is found that inclusion of the drift terms leads to no significant increase 


in the precipitated fraction at any time for any of the values of a? adopted. 
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Fig. 1 


ail) 


1.0 
M=Dt/R? 


The precipitated fraction of solute as a function of reduced time M= Dt/R?*; 
pure diffusion flow (L=0). 


DAIFT 


DIFFUSION 


0.1 


0-01 
0.01 0.1 1-0 10 100 
M 


Ratio of drift contribution to diffusion contribution as a function of reduced 
time and «?(L/r,=0-2). 
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However, as one might expect, the relative contribution of drift to diffusion 
can be quite large in the early stages of precipitation as illustrated in fig. 2, 
where the ratio of the drift and diffusion contribution from eqn. (32) has 
been plotted against M. 

V(t)/V(co) from expression (31) has also been computed and shows no 
significant difference from W(t)/W(0co) shown in fig. 1. 


Fig. 3 


V(t) 
TICR°- 1 F )Co 


10" 1.0 10 102 10? 
m=Dt/rd 


Quantity of solute entering isolated dislocation, normalized for link up 
procedure. 


———— Diffusion flow only (L=0). 
— — — — Diffusion plus drift (2 =0-2). 


For very low dislocation densities, we require the solution (38) for an 
isolated dislocation and also for the small M values shown in fig. 1 (after 
suitable normalization). Figure 3 shows expression (38) with L|ro= 0 
and 0:2 respectively ; from these curves it is again apparent that the drift 
contribution is important during the early stages. 
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It is important to know the sensitivity of the curves in fig. 1 to possible 
variations in the dislocation density, specified by the ratio K. The three 
curves shown in fig. 1 for the low a? values may be closely approximated 
to by the first term in the infinite series (32) since most of the variation 
shown by these curves is over the range M>1. In fact, for a? 50-1 we 
may write 

W(t) 
W(c) 


=1—exp[—A,?1/] 
; 4] 
where A,” is given by a 
Ao? & 2[In (1/K) — $4 2/02]. 

It is clear that for the above small values of «2, A)? is not very sensitive to 
variations in K ; however, when a? is very large, variations of K will have 
a much greater effect on the values of A,” and hence on the value of M for 
which W(t)/W(co) attains a given value. This latter result, for infinite a, 
has been discussed previously by Ham (1959). It is also worth noting 
that if the dislocation density is reduced then because the value of W(c) 
is increased, the ratio of the drift to diffusion contributions, as shown in 
fig. 2, will also be reduced at_a given M value. 


§ 6. Discussion 


We now consider briefly some of the available experimental results for 
quench and strain ageing in relation to the above theoretical treatment. 

The quench ageing results obtained by Thomas and Leak (1955) have 
been replotted in figs. 4 and 5 where the vertical scale has been renormalized 
for all the curves such that the variation is between one and zero. The 
necessity for such renormalization indicates that there is a measurable 
non-zero equilibrium solubility of the impurity at these high temperatures. 
In these figures an arbitrary time scale has been used since a plot against 
reduced time (proportional to Dt where D has been given by Wert (1950) 
for carbon in iron, and Thomas and Leak (1954) for nitrogen in iron) does 
not bring the various curves into coincidence ; in fact, the relative separation 
of the curves as a function of temperature is much less than it would be 
if the ageing were a simple diffusion limited process. If the precipitation 
were controlled by a single mechanism, then the experimental results 
would indicate much reduced activation energies compared with those of 
the diffusion coefficients. A possible conclusion is that a finite rate 
limiting process occurs at the precipitate matrix-interface. That this is 
so can be seen from eqn (41) since for small values of a? the explicit 
dependence of W(t)/W(co) on D disappears. 

For nitrogen in iron (fig. 4) the results at 98°c lie on a very good straight 
line while the results at 84°o deviate somewhat from a straight line. For 
carbon in iron (fig. 5) the results at 109°, 300°c and 397°c all lie on 
straight lines; the latter two show a small zero error on the time scale of 
about 30sec. Such small time errors have been observed by Wert (1949) 
in similar experiments when high annealing temperatures are involved. 
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The curious behaviour at 200°c is completely inexplicable and is quite 
out of keeping with the results at the other temperatures. 

The linear variations found for the experimental results shown in figs. 
4 and 5 are consistent with the form of the theoretical curves in fig. 1 for 
small a”; the curves for a?=0-1 and 0:03 have been replotted in fig. 4. 
To compare theory and experiment, it is necessary to know either the 
dislocation density or a? to determine the other uniquely, since both these 
ie are directly involved in the evaluation of the controlling first root 
of (28). 


Fig. 4 


TIME (arbitrary units) 


Experimentally observed kinetics of formation of nitride precipitates in iron 
(Thomas and Leak 1955). 


In this discussion of quench ageing, we have assumed, by our adoption 
of cylindrical geometry, that all the precipitates are nucleated on the 
dislocations. However, if the dislocation density is very low, then it is 
probable that a significant proportion of precipitate nuclei form at random 
sites throughout the crystal. This could provide an explanation for the 
different precipitation kinetics, more characteristic of spherical rather 
than cylindrical geometry, observed in the quench ageing experiments 
of Wert (1949). 
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The kinetics of impurity depletion during strain ageing experiments do 
not have the form of the curves given in fig. 1 but are better represented 
by the semi-empirical formula proposed by Harper (1951): 


EO ee 9)1/3 2/3 

Weay 1—exp [ — 3p(n/2)"°(LZD1)""] 
where p is the dislocation density. Any explanation of this variation must 
still await a treatment which incorporates a large drift potential or a 
transient variation in the value of the rate-limiting parameter or possibly 
both these factors. Such a treatment requires a completely numerical 
approach and will be presented in a subsequent paper. 


Fig. 5 


W(t) 
Wee) 


(- 


0 \ 2 3 4 
TIME Carbitrary units) 


Experimentally observed kinetics of formation of carbide precipitates in iron 
(Thomas and Leak 1955). 
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The Hall Coefficient of Liquid Mercury +; 
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ABSTRACT 


A technique for measuring the Hall coefficient of liquid metals is described. 
Results for solid and liquid mercury are presented and discussed. The Hall 
coefficient for liquid mercury is found to be —75x10-5e.m.u. over the 
temperature range —38-9°o to +100°c. 


§ 1. IntRoDUCTION 


On both the experimental and theoretical sides, knowledge of electrical 
properties is far less complete for liquid than for solid metals. 
Observations of de Haas—van Alphen effect or other properties which 
have been significant for the understanding of electronic states in solids 
are not possible for liquids, while theoretical accounts of the change (if 
any) in electronic states at and above the melting point are still 
undeveloped. 

Although observations of resistivity, thermoelectric power, and 
galvanomagnetic coefficients cannot by themselves provide detailed 
information about electronic states, there are certain qualitative questions 
which should be answerable by observation of the Hall coefficient in 
liquid metals. These are: do the conduction electron states change 
significantly at the melting point? and do they vary with temperature 
in the liquid state? These questions are suggested by (i) the occurrence 
of the volume and structure change at the melting point; (ii) the possi- 
bility of progressive structural change in the liquid state as the tem- 
perature rises (Bernal 1959); (iii) the fact that the rate of increase of 
resitivity with temperature is frequently smaller in the liquid than in 
the solid and, though often constant, is much less than would correspond 
to pxT'°K; (iv) the fact that the absolute thermoelectric power often 
changes considerably at the melting point, e.g. in Li, Cs, Cu, Hg, Zn, 
Bi amongst others. 

Observation (iii), which has never been satisfactorily explained, can 
be exemplified as follows. Let a be the observed (dp/07')p and let a’ 
be the hypothetical value of « if, above the melting point, p were 
proportional to 7'°x. Then some values of a’/a are: Na, K, Rb, 0-7; 
Cu, Ag, Au, 1:6; Al, 1-8; In, Sn, Pb, 3 to 4; Hg, Ga, TI, Sb, Bi, € 10.5: 
The proportionality of p to 7 in the solid state, for 7’'> 0, is inferred 


+ Communicated by the Authors. 
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from the theory of electron-phonon interaction on the assumption of 
constant volume, and some observed departures from it have been 
attributed to thermal expansion. It is therefore reasonable to consider 
(dp/0T)y where that is possible. For a few liquid metals the variation 
of p with pressure is known and this enables (dp/0T)y to be calculated. 
At constant volume, o’/« is greater than stated above. In particular 
for the alkali metals it becomes approximately 1-0. (Lithium is excep- 
tional in that p increases with increased pressure in the liquid and solid 
states.) The most striking example of observation (iii) is in the Zn, 
Cd, Hg group in which (@p/07’)p is negative for some way above the melting 
point in Zn and Cd (Scala and Robertson 1953). This is also true of 
(dp/07')y in Hg (see fig. 1). 


Fig. | 


Relative Resista 


0.950 
| 


° 20 30 40 50 60 70 80 90 100 
° 
Temperature C 


The variation of resistivity of liquid Hg with temperature at constant volume 


relative to that at 0°c, zero pressure, calculated from th 
Birch (1932) and Bett et al. (1954) ed from the results of 


When dp/07' is abnormally low it has been supposed that the effective 
number of free electrons is increasing, possibly as a result of some change 
in bonding, to counteract the simultaneous increase of scattering with 
temperature (Mott 1934, Regel 1958). It would therefore be of interest 
to measure the Hall coefficient of Zn, Cd, and Hg at and above the meltin 
point to see whether this property reflects any such change The 
present paper is concerned with the method of measurement bad its 


application to Hg. A preliminary statement of results h 
appeared (Kendall and Cusack 1960). bab mee es 
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§ 2. MerHop or MEASUREMENT 


The measurement of the Hall coefficient in solids is so well known that 
it is only worth mentioning here the additional difficulties presented by 
liquids. These are : (i) definition of the shape and size of the specimen ; 
(ii) the possible occurrence of magnetohydrodynamic effects. 

Since a very thin parallel-sided layer of metal is required, several 
methods of constructing glass cells were tried. For Hg, a suitable 
method was to bond flat pyrex plates with hot-setting epoxy-resin. 
During the setting, the separation of the plates was kept under obser- 
vation by forming interference fringes in the air-film between them and 
corrections made to the parellelism when necessary. After setting, and 
before and after the electrical measurements, the separation of the 
plates was measured by means of Edser—Butler fringes. In this way, 
rectangular cells of well-defined thickness in the range 20 to 100 microns 
were obtained (see fig. 2). 


Seale 


inches Cisc 


Glass Hall cells. Ecrmiienas 
A. Two-probe cell. B. Three-probe cell. (a) Mercury specimen ; all probe ; 
Es () resistivity probe ; (d) current lead ; (¢) thermocouple ; (/) filling 
tube. 


Current leads to the cells were of tungsten wire or foil and the Hall 
probes were columns of mercury in glass capillaries cemented into the 
cells. Both 2- and 3-probe methods were used. The ends of the Hall 
probes remote from the specimen were sealed off with tungsten electrodes, 
brought close together, and kept in a Dewar vessel; this is desirable to 
eliminate voltages from the Ettingshausen and thermoelectric effects, 
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The cells were filled under vacuum with purified Hg and encased in a 
thick copper box through whose walls passed tubes for circulating 
heating or cooling fluids. The whole assembly was mounted between 
the pole pieces of an electromagnet. 

Careful filling is of considerable importance because voids in the Hg 
film will permit irreversible movements under electro-dynamic forces 
thus changing the current flow lines. The absence of voids was indicated 
by visual inspection and by reproducibility of the readings after reversals 
of current or field. 

The Hall coefficient was obtained by the standard d.c. method described 
by Lindberg (1952), the voltages being measured to about +0-02yuVv 
with a constant resistance potentiometer based on the design of Teele 
and Schulman (1939). The whole measuring apparatus was tested by 
measuring known voltages from calibrated standards and by measuring 
the Hall coefficient of pure solid copper with good agreement with cur- 
rently accepted values. 


Fig. 3 


\ 
o 


Room Temp. 


Hall E.M.E pV 


-0.5 Specimen Thickness 
— 5.51x10 37cm 
-~=-2.03x10 “cm. 


° ' 2 3 4 5 6 7 


Field kilo-Oersted 


Liquid mercury : Hall e.m.f. vs magnetic field. 
Ihe dots indicate independent measurements. 
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It is possible that the voltages at the Hall probes of a cell contain a 
component due to magnetohydrodynamic effects. It was found too 
difficult to estimate its value from magnetohydrodynamic theory and 
some empirical checks on its importance were made. The general 
problem was first discussed by Williams (1925) from whose work it can 
be inferred that the magnetohydrodynamic effect will be (i) present 
only in so far as the current density and magnetic field vary over the 
specimen; (ii) proportional to H? for constant specimen current, J; 
(ii) proportional to the linear dimensions of the specimen. The Hall 
effect on the other hand is proportional to H and inversely proportional 
to the linear dimensions. JH here refers to the total field which means 
in practice the applied field since the field of the specimen current is 
negligible by comparison. Since it is only inhomogeneities in H and I 
that produce the undesired effects, it was thought that the best check 
would be to perform the experiment with several cells of different thick- 
nesses. Other things being equal, the ratio of Hall voltage to spurious 
voltage should have varied by a factor of 25 from the thinnest to the 
thickest cell. In one experiment the pole separation and therefore, 
probably, the field inhomogeneity was varied, and different geometries of 
the current electrodes were also tried. Inhomogeneities of current 
density, though bound to exist to some extent, are certain to vary from 
cell to cell and the observed agreement of the results from cell to cell 
was regarded as good evidence that magnetohydrodynamic effects, if 
present, were not serious. It was also confirmed that the Hall voltage 
varied linearly with H for constant J and vice versa (see fig. 3). 


§ 3. ResuLTs AND COMPARISON WITH OTHER EXPERIMENTS 


Figure 3 shows a typical result for liquid Hg. The average value of 
the Hall coefficient from many measurements on four cells of thicknesses 
2-03, 4:17, 5°51, 10-5 x 10-8cm was 74:°6x 10-°+4%em.u. The variation 
with temperature is shown in fig. 4. 

Some measurements on solid Hg were also made. After each freezing, 
the readings were internally very consistent, somewhat more so than 
on a liquid specimen, but the Hall coefficient changed when the specimen 
was melted and refrozen. In this way, values ranging from —87 to 
—101x10-e.m.u. were observed. This is probably due to the lack of 
control over the crystallization; very likely the orientation of the 
crystals varied from specimen to specimen and it is highly probable 
that the Hall coefficient of solid Hg is anisotropic. This deserves a 
separate study. The suspected anisotropy makes it difficult to compare 
the present results with those of Jaggi and Sommerfelder (1959; table) 
especially as they worked at 4:17°K. The Hall coefficient may very well 
vary with temperature between the freezing point and 4-17° and the 
variation will not necessarily be the same for all orientations; such 
phenomena have been demonstrated’ in gallium (Yahia and Marcus 


1959). 
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Only a few measurements of the Hall coefficient of liquids have been 


made before but naturally Hg features in these. The table shows the 
unsatisfactory state of the subjecty. 


Fig. 4 
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Hall Coefficient CGSM x 10° 


Specimen Thickness 
+ 5.51 x 1O-cm. 
x 10.5x IO *cm. 


-50 ° 50 100 ° 
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Temperature '& 


Hall coefficient of mercury as a function of temperature. Its sign is negative. 


Value of 10°R (in e.m.u.) for mercury 


Author R solid R liquid 
des Coudres 1901 
Fenninger 1914 
Serduke and Fisher 1932 
Kikoin and Fakidow 1931 


Jaggi and Sommerfelder 1959 
Present work 


0-0 
—72 to —88 (at m.p.) ’ 


— 200 (at 4:17°K) 
—87 to —101 (near m.p.) 


Measurements on solids are less subject to error than those on liquids 
and the agreement between the present results and those of Serduke 
and Fisher is satisfactory for solid Hg. The authors therefore agree 
with Serduke and Fisher in regarding the results of Fenninger with some 


} Note added in proof.—Satisfactory confirmation of the values in the present 


paper has now been given by G. Wilson (private communication) and Y. Tiéche, 
Helv. phys. Acta, 38, 963, 1960, 
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suspicion. Kikoin and Fakidow did not measure solid Hg and could 
not detect a Hall effect in the liquid. It is difficult to discuss this result 
without detailed knowledge of the apparatus, but it seems likely that 
their specimen thickness was greater than those of the present 
experiments—a circumstance tending to greater inaccuracy. Their 
stated error in the Hall coefficient of liquid Na-K alloy was +10% 
and the voltage from a mercury specimen of similar size would be only 
17% of that from Na~K—if the present result for Hg is correct. 

The early work by des Coudres (1901) agrees with this paper in that a 
definite negative effect was found. There is no point in discussing the 
magnitude because other of des Coudres’ results, on solid materials, 
are now known to be very inaccurate quantitatively. 

The present experiments have the advantages that the specimens were 
thinner and of more accurately determined thickness, that more than 
one cell was used, and that the temperature was varied. In view of 
the disagreements in the table, it would be desirable to repeat the 
measurements using an a.c. technique which, though involving difficulties 
of its own, offers another way of minimizing the specific disadvantages 
of working with liquids. 


§ 4, DISCUSSION 


The sign of the Hall coefficient in solid Hg is negative. In the course 
of discussing the bonding of Group IIB elements, Wallace (1955) put 
the hypothesis that 0-67 valence electrons per atom participate in 
covalent bonds in the basal planes, leaving 1-33 electrons per atom as 
conduction electrons contributing to metallic binding. It was sub- 
sequently inferred that the Hall coefficient would be positive, in agree- 
ment with experiment for Zn and Cd. The negative value for Hg, 
while by no means sufficient in itself to disprove this theory of bonding, 
certainly casts some doubt on it. 

The small change in the Hall coefficient at the melting point is con- 
sistent with the constancy of the Knight shift across the melting point, 
since the latter argues that the electron states cannot alter much on 
melting (Knight et al. 1959). 

In a solid metal the ‘free electron value’, —1/nle|c, of the Hall 
coefficient obtains if (i) the Fermi surface is spherical and (ii) the electron 
scattering is independent of k for a given energy. Naive considerations 
of the structural disorder in a liquid metal might encourage the expecta- 
tion that these conditions should be fulfilled, since their contravention 
in solids arises from the crystal structure. However, it is not yet clear 
whether liquids are sufficiently like solids for similar electronic terminology 
to be used, e.g. it is not certain whether it is permissible to speak of a 
Fermi surface. Assuming that it is, one may ask to what extent the 
electronic properties of liquid mercury are accounted for by a free electron 
model with two electrons per atom. 

On this model, the Hall coefficient would be —77-4x 10-°e.m.u. 
which agrees within experimental error with the observed one, The 
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electron mean free path at the melting point would have to be 7:44 
which is smaller than in many other metals but not of a different order. 
The Lorenz number would be 2-25 x 10-8 watt ohm deg! and independent 
of temperature, whereas the observed value is 2:62+0-03 x10-® from 
100 to 300°c (Ewing et al. 1955). The absolute thermoelectric power 
should be negative and proportional to the absolute temperature with 
an approximate value of —2-9uvdeg-! at 0°c; the observed values 
vary approximately linearly with, but not proportionally to, the absolute 
temperature and the 0°c value is — 4-4yv deg? (Gehlhoff and Neumeier 
1919). According to Schulz (1957), though not according to Hodgson 
(1959), the optical properties of liquid Hg are also very well accounted 
for by the Drude free electron theory with two electrons per atom. It 
seems therefore that liquid Hg shows free electron behaviour to a 
considerable degree. 

However, between the melting point and 100°C, the observed 
resistivity of Hg increases by a factor 1-14; if the resistivity were 
roughly proportional to absolute temperature, as occurs commonly in solid 
metals (including Hg) and also in liquid Na and K, the factor would be 
about 1:59. To attribute the difference to a rise in the effective number 
of free electrons requires the latter to increase by 40% in 140°c. As is 
clear from fig. 4, the present experiments fail to reveal any comparable 
variation in the Hall coefficient which remains at about the ‘free electron ’ 
value throughout. It would be of interest to know if similar constancy 
characterizes the Knight shift. Thus either the effective number of 
free electrons must increase considerably without alteration in the Hall 
coefficient or the scattering of electrons by the disordered structure is 
sufficiently different from ordinary electron—phonon interaction for a 
much smaller rate of increase with temperature to be possible. The 
first alternative, though not impossible, is less attractive because both 
quantities depend on the Fermi surface. The second is in line with the 
view that liquid metal resistivity has much in common with a residual 
resistance created by lattice defects, this disorder entering at the melting 
point (MacDonald 1959, Cusack and Enderby 1960). It remains of 
interest to find whether the Hall coefficient is also constant when the 
resistivity actually decreases with rising temperature as in Zn. Such 
work is in progress but the higher melting point considerably increases 
the technical difficulty. 

If the resistivity increment, Ap, on melting of Hg is of the ‘residual’ 
type, the abnormally large ratio of resistivities, p,/p,, at the melting 
point (~4 instead of ~2) is a simple consequence of the abnormally 
low melting point. For if the solid state were to persist to ~500°c 
p, would be trebled and Ap remain about the same. However, other 
explanations of the high p,/p, have been offered, e.g. a decrease in the 
effective number of free electrons on melting (Mott and Jones 1936) 
extra scattering mechanisms in the liquid state (Wallace 1955) and it is 
not possible to decide certainly between them. 
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ABSTRACT 


Etching experiments are described in which diamonds showing octahedral, 
cubic and dodecahedral faces are heated in air within the temperature range 
of 800°c to 1400°c. A change in the etch pit orientations of the characteristic 
etch pits for these faces has been found as the temperature of etching is 
raised. Etching with mixtures of gases has led to the conclusion that oxygen 
and water vapour are the relevant constituents of air in the formation of well- 
defined etch pits. Surface graphite formation at temperatures below 1400°c 
is considered and it is proposed that the mechanism of formation is essentially 
a chemical reaction and not connected with a purely physical phase change. 


§ 1. IntTRODUCTION 


Most of the previous work on the etching of diamond surfaces has been 
done with fused potassium nitrate or sodium carbonate as the etchant 
(see, for instance, Williams 1932, Omar et al. 1954, Patel and Tolansky 
1957) in which characteristic etch pit shapes for the octahedral, cubic and 
dodecahedral faces have been produced. Omar and Kenawi (1957) how- 
ever formed well-defined pits of triangular outline by heating an octa- 
hedron on a molybdenum strip to about 1000°c in the presence of air at 
2-5 microns Hg pressure. Using the nomenclature of Frank et al. (1958), 
these pits were in the positive orientation, that is, in the same orientation 
as the octahedral face which is opposite to the negatively oriented natural 
‘trigons’. On an attempt to produce negatively oriented etch pits, 
Frank and Puttick (1958) heated octahedra in steam, chlorine and wet 
hydrogen at 800°c but did not succeed in doing so. They did however, 
produce such pits by heating in fused kimberlite at about 1450°c. 

In the experiments to be described, the above work has been extended 
by etching diamonds in air at reduced and at atmospheric pressure in the 
temperature range of 800°c to 1400°c. As the temperature of etching is 
raised, the etch pit orientations have been found to change on the octa- 
hedral, cubic and dodecahedral faces. In an attempt to determine which 
of the constituents of the air are responsible for the observed etch pits, 
octahedra have been etched by heating in various combinations of the 
gases found in the atmosphere. 
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§ 2. EXPERIMENTAL PROCEDURE 

The diamonds primarily used in the etching experiments were 0-1 carat 
octahedra which were selected for their surface perfection. The cleaning 
procedure prior to etching consisted of immersing the diamond in boiling 
concentrated nitric acid, washing in distilled water and refluxing in ethyl 
alcohol. The stones were handled throughout with glass tweezers to 
avoid contamination. For etching experiments at temperatures less than 
1000°c the specimen was placed in a silica tube, closed at one end, which 
was evacuated to 10-*mm Hg after flushing with argon. At temperatures 
greater than 1000°c alumina tubes were used. The stone was heated to 
the etching temperature in this vacuum which served to complete the 
cleaning procedure. Control experiments showed that this heating 
operation did not result in any detectable etching of the surfaces. The 
air was then bled into the tube at a controlled rate so as to maintain the 
required etching pressure whilst continuously pumping. On completion 
of the run the tube was again evacuated to 10-*mm Hg and the specimen 
allowed to cool. An optical microscope examination was then made of 
the surfaces. 

The main constituents of air are approximately 78% by volume of 
nitrogen and 21% by volume of oxygen. Water vapour is also present 
at concentrations usually less than 1%. The inert gases and the trace 
quantities of hydrogen and carbon dioxide have been found not to influence 
the etching significantly. In experiments on etching with mixtures of 
gases, two flasks were evacuated to 10-*mm Hg, one filled with nitrogen 
to atmospheric pressure, and the other with oxygen to the same pressure. 
The volumes of the flasks were in the ratio of 4:1 so as to simulate the 
nitrogen/oxygen ratio in air. The flasks were interconnected and the 
mixture fed into the heated silica or alumina tube, which had previously 
been evacuated and isolated from the pumping unit. After the required 
etching time, the tube was again evacuated to 10-*mm Hg and allowed to 
cool. In this case, the mixture in the tube was at approximately atmo- 
spheric pressure during the etching. With this system it was possible to 
etch separately with oxygen, nitrogen, water vapour, oxygen and nitrogen 
mixtures, oxygen, nitrogen and water vapour mixtures at the partial 
pressures at which these gases were found in air. In some experiments 
it was desirable to substitute an inert gas in place of the nitrogen, and 
argon was used. The nitrogen, argon and oxygen were obtained from 
cylinders, and these gases were separately dried and impurities removed 
before admission into the flasks by passing the gases through: 

(1) a horizontal column designed for turbulent flow containing soda 

lime to remove carbon dioxide ; 

(2) heated pailadized asbestos for removal of hydrogen; 

(3) a horizontal tube containing coarse silica gel which removed most 

of the water vapour ; 

(4) a heat exchanger type of cooling trap at —80°o; 

(5) a cold trap at —80°c filled with fine silica gel to remove organic 

matter and water vapour ; 
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(6) a turbulent flow tube containing magnesium perchlorate for the 
final drying. 
When the gases were required to be saturated with water vapour, they 
were subjected to the above sequence and then finally passed through two 
water bubblers. 

At etching temperatures of 900°c and above, the diamond surfaces were 
covered with a layer of graphite which was removed by heating for 10 min 
in a mixture of perchloric, nitric and sulphuric acids before the surfaces 
were examined. 

Although most of the experiments were done with octahedral stones, 
cubes with {100} surfaces and also, cubes with a mixture of {100} and {110} 
surfaces were investigated. 


§ 3. RESULTS 


In the first series of experiments octahedra were etched at temperatures 
between 800°c and 900°c in air at pressure ranging from 10-?mm to 
atmospheric pressure. The severity of the attack increased with 
increasing pressure, temperature and time, but the sequence of etching 
appeared to be the same at these ranges of pressure and temperature. 
First of all, the microcracks and edges were attacked more vigorously than 
the rest of the surface, but at a later stage the surface became covered with 
etch pits of triangular outline in the positive orientation. Figure 1+ shows 
such pits on an octahedron which had been etched for 15min at 850°c at 
an air pressure of 760mm Hg. 

When the same experiments were repeated at temperatures in excess of | 
1000°o, the stones were covered with a layer of graphite. Removal of the 
graphite and examination showed that the same sequence of etching 

_occurred as before except that the pits were now in the negative orientation, 
i.e. the same as for natural ‘trigons’. Figure 2 shows these pits on an 
octahedron etched for 5min at 1250°c at atmospheric pressure. 

It was thus interesting to investigate the region between 900°c and 
1000°c which would be the region where reversal of pit orientation occurs. 
An octahedron was etched for 1 hour at 950°c in air at a pressure of 10-'mm 
of Hg, and the resultant etch pits are shown in fig. 3. They were primarily 
circular pits with a tendency towards the hexagonal outline which would 
be expected in an intermediate region where reversal of orientation of 
triangular etch pits takes place. 


3.1. {100} Surfaces 


The etching of these surfaces could not be investigated at various 
pressures and temperatures to the same extent as {111} surfaces, owing to 
the comparative scarcity of the polished cubes which would be necessary. 
However, a cube with {100} surfaces was heated at 775°c for 2 hours at an 
air pressure of 5x 10-'mm of Hg. The etch pits formed were of square 


+ Figures 1-7, 9, 10 and 11 are shown as plates. 
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outline such that the sides were in (100) directions. Figure 4 shows these 
pits. A second cube was heated at 1400°c in air for 6 min at atmospheric 
pressure and the surface graphite removed. Again there were pits of 
square outline on the surface but now the sides were parallel to ¢110), 
i.e. the high temperature square pits were rotated through 45° relative to 
the low temperature pits. The orientation of the high temperature pits 
was the same as those found on the untreated natural cubic crystals which 
occur in nature. Figure 5 shows the high temperature pits on a (100) 
surface. 


3.2. {110} Surfaces 


In the case of dodecahedral faces, difficulty was experienced in forming 
etch pits and prolonged etching times were required to obtain recognizable 
pits. However, by etching in air at atmospheric pressure for 2} hours at 


Fig. 8 
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800°c ‘boat-shaped’ pits were formed such that on a (110) face, the keel 
of the boat was in the [001] direction (fig. 6). In the corresponding 
experiment at 1400°c, etching in air at atmospheric pressure for 10min 
produced ‘boat-shaped’ pits with the keel in the [110] direction on the 
(110) face (fig. 7). A graphite layer had to be removed before the diamond 
surface was revealed for examination. Thus for the dodecahedral face as 
well, the high temperature pits were rotated relative to the low temperature 
pits: through 90° in this instance. : ee: 
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Summarizing these results, fig. 8 shows the etch pit orientations on the 
three faces which have been considered. 


3.3. Graphite Formation 


The etching experiments at 900°c and above with octahedral stones 
showed that the formation of graphite on the surface depended, not only 
on the temperature, but also upon the pressure of the air surrounding the 
diamond. Heating for 16 hours at an air pressure of 10-*mm of Hg 
resulted in a faint grey graphite layer with slight attack at microcracks 
and edges of the underlying surface. However, heating for 5min in air 
at atmospheric pressure provided an appreciable opaque black layer of 
surface graphite at 900°c and above, with a corresponding heavy attack 
of the diamond surface. Reflection electron diffraction revealed that the 
crystallite size varied from less than 504 in the case of rapid graphite 
formation, to greater than 200A when formed at 1400°c at 10-*mm air 
pressure. In the cases examined, namely octahedral surfaces, oriented 
growth of the graphite had occurred with the graphite layers lying parallel 
to the puckered (111) surface. 

Heating octahedra in an atmosphere of carbon dioxide also produced 
surface graphite at temperatures in excess of 900°C with heavy attack of 
the underlying diamond surface. Carbon monoxide however, did not 
alter the surface of the diamond appreciably, even after heating for 45 min 
at 1400°c at atmospheric pressure. 


3.4. Etching with Mixtures of Gases 


The use of air as the etchant has been seen to produce well defined etch 
pits on the three surfaces considered. It cannot be assumed that the 
oxygen in the air was solely responsible for the etch pit formation as 
current theories (Frank 1958, Cabrera and Vermilyea 1958, and Gilman 
et al. 1958) postulate that an etchant which produces well-defined etch 
pits has at least two consituents, namely the substance which attacks the 
surface and an inhibitor which modifies the attack so that pits are formed 
instead of, for instance, surface polishing occurring. In order to separate 
the roles of the constituents of air in the etching process on octahedral 
surfaces of diamond, stones were etched in different mixtures of the gases 
found in air. The pressures of the constituents in these mixtures were 
the same as the partial pressures of the same gases in air. 


3.4.1. Low-temperature etching with gas mixtures 
In all these experiments, the etching time was 12min at 850°c. 


(a) Etching at atmospheric pressure in a mixture of 4 parts nitrogen 
plus 1 part oxygen plus water vapour, produced positively oriented tri- 
angular pits over the surface as in fig. 9. These showed all the characteri- 
stics of etch pits formed in the previous experiments with air. 
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(b) Etching with dry oxygen at 170mm of Hg pressure resulted in a 
rough matt-type surface with no triangular pits (fig. 10). An important 
and significant observation however, was that the microcracks had been 
widened by etching to form large deep positively oriented pits of tri- 
angular outline. Etching at atmospheric pressure in a dry mixture of 
4 parts argon plus 1 part oxygen gave essentially the same result. Thus 
the pits formed by air etching could not have been due to the action of 
oxygen alone. 

(c) Etching in wet oxygen at 170mm of Hg pressure produced well- 
defined positively oriented etch pits and these could not be differentiated 
from those produced in air. The same type of etch pits were obtained by 
etching in 4 parts argon plus | part oxygen plus water vapour at atmo- 
spheric pressure. 

(d) To determine the effect of nitrogen, an octahedron was etched in a 
dry mixture of 4 parts nitrogen plus 1 part oxygen. No etch pits were 
formed and the surface was roughened in the same manner as when dry 
oxygen was used as the etchant. 

Control experiments were done on the etching effect of nitrogen and 
water vapour separately. Treatment with nitrogen alone at 590mm of 
Hg pressure did not produce any detectable attack on the diamond surface ; 
neither did water vapour at 18mm of Hg pressure under similar con- 
ditions to the etching experiments described in this section. 


3.4.2. High-temperature etching with gas mixtures 


The above experiments in the previous section were repeated at 1250°c 
with an etching time of 5min. Surface graphite had to be removed before . 
the surface could be examined. 


(a) The wet nitrogen and oxygen mixtures produced well-defined 
triangular etch pits in the negative orientation (fig. 11), which were 
indistinguishable from those obtained by air etching. 

(b) The dry oxygen and also the dry nitrogen and oxygen mixture 
resulted in rough pitting with no crystallographically oriented etch pits. 

(c) Wet oxygen produced well-defined etch pits in the negative orien- 
tation, similar to those formed in (a) above. 

The experiments in the previous two sections showed that the oxygen 
and the water vapour were the constituents responsible for the etch pit 
formation observed when diamonds were heated in air. 


3.5. Htch Rates 


The effect of water vapour on the etch rate of oxygen on octahedral 
diamond surfaces was determined by measuring the rate of attack in 
both dry and wet oxygen in the closed system used in the mixture experi - 
ments. ‘The presence of water vapour reduced the rate of oxygen attack 
by approximately four times. This result is at variance with the work of 
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Hennig (1960) who found that water vapour acted as a catalyst to the 
reaction of oxygen with the surface of single crystals of graphite. Work 
done on the etch rates of diamonds in a fast-flowing stream of oxygen also: 
showed that under these conditions water vapour acted in such a way as 
to double the etch rate with wet oxygen over that of dry oxygen (Evans 
and Phaal, to be published). This apparent contradiction in results can 
be explained by assuming that water vapour accelerates the oxygen— 
carbon reaction, but that one or more of the products of the reaction acts 
as an inhibitor. In the case of flow experiments, the gaseous products of 
the primary reaction tend to be swept from the surface as they are formed, 
whereas in the closed system, secondary surface reactions take place 
resulting in inhibition to the oxygen—carbon reaction. 


§ 4. Discussion 


Frank (1958) and Cabrera and Vermilyea (1958) ascribe the growth of 
crystals from solution to the movement of atomic steps over the surface. 
They also apply the theory to dissolution and propose that in formation of 
well-defined etch pits, the etching solution contains a solvent and also an 
impurity which inhibits the movement of steps, causing ‘bunching’ of 
steps to occur. The results of Gilman ef al. (1958) confirm this theory for 
the etching of lithium fluoride in which they used water as the solvent 
and ferric ions as the inhibitor in the formation of well-defined etch pits. 
The experiments described have shown that the theory can also be 
applied to etching by gaseous mixtures. In the case of air etching on 
diamond surfaces, the carbon atoms are removed by the oxygen—carbon 
reaction which is accelerated by the presence of water vapour on the 
surface. However, one or more of the primary products of the water— 
carbon reaction, if it is not removed, can accumulate on the surface and 
act as inhibitor to the oxygen—carbon reaction, resulting in well-defined 
etch pits. This also explains the difficulty which we have found in forming 
good etch pits in a flow type experiment. Extensive work has been done 
on the steam—carbon reaction, and most investigators agree that the 
primary reaction is H,O+C——>H,+ CO (see for instance Binford and 
Eyring 1956); Long and Sykes (1948) consider that hydrogen molecules 
retard the steam—carbon reaction by being strongly adsorbed on the surface. 
Key, quoted by Strickland-Constable (1950), on the other hand, considers 
the retardation as due to the reaction of hydrogen gas with surface oxide 
to give H,O gas. Whatever the details of the surface reactions, it is 
probable that in the air etching experiments on diamond surfaces, the 
primary reaction product which acts as an inhibitor to the oxygen—carbon 
reaction is associated with the hydrogen. Experiments on etching with 
dry hydrogen and oxygen mixtures have failed to produce well-defined 
etch-pits on octahedral surfaces, but it was not possible to vary the 
hydrogen concentration widely owing to the explosive nature of the 
mixture at elevated temperatures. 


2H2 
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4.1. Htch Pit Formation 


If an etch pit is considered as being formed by the movement of atomic 
steps over the surface, the resultant pit is outlined by the steps moving at 
the slowest rate. Thus the observed change in etch pit orientation can be 
explained by a change in the relative velocities of different types of steps, 
under the varying experimental conditions. Consider the etching of an 
octahedral face where, of course, the triangular {111} faces are outlined by 
(110) directions. Frank, Puttick and Wilks (1958) describe the atomic 
structure of steps formed in the (110) directions. Take, for the sake of 
illustration, an etch pit as in fig. 12—which is initially hexagonal in outline 
and bounded by <110) directions. Steps 1 are atomically identical as are 
steps 2, and Frank ef al. (1958) point out that the edges of steps 1 consist 
of atoms which are trebly bonded to the lattice, whereas atoms at steps 2 


Fig. 12 
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are doubly bonded to the lattice. When steps 1 move at a higher velocity 
than steps 2, positive etch pits result, whereas when the reverse is the case, 
negative pits are formed. The observation that oxygen attack alone 
opened microcracks to form deep pits of positive orientation at 850°c, 
indicates that the relative step velocities are determined by the relative 
stabilities of the different oxygen—carbon complexes at the steps and is 
distinct from the action of the inhibitor. Information is lacking about 
the type of complexes formed between oxygen and carbon on diamond 
surfaces, but steps 2 must be more strongly stabilized by oxygen than 
steps 1 at temperatures below 900°c, The tentative explanation by Frank 
et al. that steps 2 are strongly stabilized by oxygen atoms forming bridges 
between carbon atoms as in a lactone, appears reasonable by our experi- 
ments. No similar stabilization is possible on steps 1 as the free bonds are 
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not favourably oriented for bridge formation. At temperatures exceeding 
1000°c, the observation of negative pit formation can be explained by a 
reduction in the stabilization by bridge formation at these temperatures 
such that the lifetime of the oxygen—-carbon complex on steps 2 is less 
than the different type of complex formed on steps 1. At temperatures 
between 900°c and 1000°c the lifetimes of the complexes at the two types 
of step are the same and circular pits result. 

Now the (100) surface will be considered. It is not possible to describe 
the results by considering that the kinetics of dissolution is determined by 
the number of bonds per carbon atom which has to be broken during 
dissolution. If the diamond structure is considered in terms of (400) 
planes, the (100) surface consists of atoms which are doubly bonded to the 


Fig. 13 


lattice. Removal of these atoms reveals a plane at distance ja[100] 
below the original surface which again consists of doubly bonded atoms 
and so on. Step formation and also kink formation at steps on such a 
surface do not produce favoured sites for dissolution as the terminal 
atoms are also doubly bonded to the lattice. To account for the fact that 
etch pit formation is readily obtainable on {100} diamond surfaces, it is 
necessary to postulate that the flat surface is strongly stabilized by adsorbed 
atoms so that sites which are kinetically favoured can be formed at steps. 
Examination of the free bonds on a flat (100) surface, shows that bridge 
formation by adsorbed oxygen atoms between carbon atoms is possible 
because the free carbon bonds are favourably oriented and spaced for the 
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formation of such a complex. Thus a flat (100) surface can be stabilized 
in the same way as steps of type 2 on octahedral surfaces since the free 
carbon bond configurations are similar in both cases. If the flat (100) 
surface is assumed to be strongly stabilized by bridge formation of oxygen 
atoms between carbon atoms, steps can again be considered as positions 
where it is kinetically favourable for etching to occur. Consider the 
initial etch pit to be bounded by steps of height }a[100] and in directions 
shown in fig. 13. Steps 3 are identical atomically as are steps 4. Low 
temperature pits will form if steps 3 have a lower velocity over the surface 
than steps 4, and high temperature pits will be produced when the reverse is 
the case. Examination of the atomic configurations at these steps shows 
that if the rates of step movement are determined by the number of bonds 
per carbon atom which have to be broken during dissolution, then steps 
3 and 4 would move with the same velocity. This is clearly not so by 
the experimental evidence. The relevant asymmetry in the two types of 
steps is the distribution of free bonds at the terminal atoms. These free 
bonds determine the type of oxygen—carbon complexes which form at 
the steps. It is postulated that the complexes on steps 3 are more stable 
than those formed on steps 4 at low temperatures (i.e. less than 900°c) 
and the reverse occurs at temperatures greater than 1000°c. This could 
then account for the change in etch pit orientation as the temperature is 
raised. 

As mentioned before, difficulty was found in forming boat-shaped pits 
on {110} surfaces, and considerable roughening of the surface was necessary 
before they developed. Williams (1932) in his book also states that the 
boat-shaped etch features are sometimes raised hillocks on the dodeca-_ 
hedral face. It is considered that the boat-shaped pits are not formed 
as a result of step movement over the surface, but are due to the tendency 
for depressions formed during roughening by etching to stabilize themselves 
by assuming a shape bounded by surfaces with the slowest etch rates. 
This necessitates that there is a change in the relative etch rates of 
different diamond surfaces as the temperature is raised, and experiments 
are being done to test this. 


4.2. Graphite Formation 

As stated in the results, diamond octahedra which had been heated in 
air between 900°c and 1400°c had a carbon layer on the surface. The 
structure of the layer varied from effectively amorphous carbon to an 
oriented layer of graphite with a crystallite size of greater than 2004. 
The presence of this layer depended markedly upon the atmosphere 
surrounding the stone during the heating operation. In a high vacuum 
(10-* mm of Hg) very little formation of surface carbon resulted, whereas 
heating in air at atmospheric pressure caused the stone to be covered with 
an opaque carbon layer at 900°c. This has led to the conclusion that in 
the temperature range of 900°c to 1400°o the formation of surface graphite 
is due to the deposition of carbon on the diamond surface as a result of 


Etching of Diamond Surfaces with Gases 439 


surface chemical reactions and the effect of a purely physical phase change 
is not evident. Work on graphite, summarized by Nelson Smith (1959) 
has shown it is possible to deposit carbon on the surface by a secondary 
surface reaction in which carbon monoxide reacts with a surface oxide 
complex to form carbon dioxide and carbon, i.e. 


CO +(CO)>CO,+C. 


This is consistent with our results, as the primary reaction product of the 
oxygen—carbon reaction is presumably CO which at a sufficiently high 
temperature can react with the (CO) surface complexes (Nelson Smith 
1959). Also no graphite formation would be expected by heating in carbon 
monoxide alone as no surface oxide complexes would be formed with which 
it could react. The graphite formation as a result of heating diamonds 
in carbon dioxide can also be explained by the same mechanism following 
the primary reaction 


Ce CeO 100). 


The thickness of the carbon layer under particular experimental conditions 
would be an equilibrium thickness determined by the rate of formation 
of the deposited carbon as a result of these secondary surface reactions and 
the rate of removal of the layer by oxidation. It would be expected that 

the structure of the layer would play an important role in the equilibrium 
thickness as this would determine the effectiveness of the layer in protecting 
the underlying diamond surface from the surrounding oxidizing atmosphere. 
It must be emphasized that the experiments were done at 1400°c and 
below and the interpretation does not necessarily hold for temperatures 
in excess of 1400°c. 
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ABSTRACT 


A study has been made of the development of etch pits and dehydration 
nuclei on the (010) face of crystals of gypsum. The etch pits marked the 
sites of features that extended into the crystals. Dehydration nuclei in 
general did not occur at the same sites as etch pits. The nature of the sites 
favourable for the formation of etch pits and of dehydration nuclei is discussed. 


CHEMICAL reactions and phase changes in solids are likely to be initiated 
at particular sites such as lattice imperfections. This applies also to the 
dissolution of crystals, which probably begins at unit pits nucleated 
preferentially at the points of emergence of dislocation lines, with removal 
of the matter of the crystal proceeding by movement of monomolecular 
steps across the surface. 

This note presents the results of a study of dissolution and thermal 
decomposition at the (010) face of naturally occurring, high quality 
crystals of gypsum. When the faces, which were produced by cleaving, 
were washed with distilled water etch pits were formed on parts of the 
surface. Some of these pits occurred in rows (cf. Grengg 1914, Viola 
1897). Washing the crystal with a 1% solution of sodium hexameta- 
phosphate, which is known to have a marked effect on the crystallization 
of gypsum (Ridge and Surkevicius 1960), made the pits more pronounced 
and somewhat increased their number. It also brought about the 
formation of geometric features on the parts of the surface free from 
etch pits. In general, the full crop of pits appeared within a minute or 
two and was not greatly increased by further washing, showing that the 
probability of forming a pit was high over the population of sites favourable 
for the formation of pits. 

It was found that if a thickness of 10, the approximate depth of the 
largest pits, was removed from the surface by grinding, rewashing 
frequently produced the original pattern of pits. Indeed, even after the 
removal of 60 1 some of the features of the original pattern of pits could 
often be distinguished. This shows that the pits occur at fixed sites 
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which extend into the crystal. In a particular case the x-ray technique 
devised by Schulz (1954) revealed a sub-grain boundary at which the 
change in orientation was about 1°; on etching this boundary was marked 
by a dense line of pits. These facts suggest that in the crystals examined 
the etch pits do mark the points of emergence of dislocation lines, and 
the rows of pits the positions of sub-grain boundaries, although some of 
the observations are difficult to reconcile with this view. The pits show a . 
density of about 105 per em2, which is probably several orders of magnitude 
less than the density of dislocations in such crystals. Furthermore, 
examination of the crystals by the Schulz x-ray technique revealed sub- 
grain boundaries in parts of the crystals where the etch pits did not form 
rows. Some of the rows of etch pits terminated at the edge of the crystal 
or at junctions with other rows of pits. On the other hand some rows 
terminate in the crystal, as is shown in figs. 1 to 47. These facts show 
that etching did not reveal all dislocations cut by the surface. Breaks in 
rows of etch pits often occur in the region of cleavage steps (see fig. 5), 
and demonstrate that some parts of the features promoting the formation 
of rows of pits are immune to etching. It is possible that such effects 
may be due to variations in the distribution of impurities in the crystal, 
or to the capacity of dislocation lines to produce pits depending upon their 
direction. 

Focusing below the surface of the crystal showed that some, but by no 
means all, of the rows of etch pits were associated with walls of erystal-like 
inclusions. No tubular features suggesting hollow dislocations (Verma 
1953) were seen. 

Substances promoting the formation of etch pits probably act by 
being adsorbed on the edges of the receding sheets of ions, and hinder 
horizontal dissolution (Gilman et al. 1958). For this reason substances 
that modify habit seem likely to promote the formation of etch pits. 
Other habit modifiers for gypsum were tried, including acetates, citrates 
and borates, but sodium hexametaphosphate was the most effective. No 
reagent tried was effective in filling in the gaps in the rows of pits. The 
relative effectiveness of distilled water suggests that the crystals had 
grown in an environment in which substances capable of adsorption on 
the growing surface were present. 

Holding crystals at 98°c for a few hours in an air oven resulted in the 
formation of areas of slow decomposition which began as rectangular 
nuclei, the sites of which did not generally coincide with those of etch 
pits, although some cases of coincidence were observed. Unlike the 
etch pits, further treatment (heating) caused the dehydration nuclei to 
grow and new ones to appear continuously in the empty parts of the 
crystal. When the nuclei had grown sufficiently to begin to overlap they 
were about two orders of magnitude less numerous than the etch pits. 
The nuclei were not of quite the same shape as those obtained by 
ee eee ive Sete es 

{ All figures are shown as plates. 
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Zolotov (1958) on cleaved, unwashed surfaces, and the appearance varied 
a little in crystals from different sources. Re-etching the surface after 
formation of the dehydration nuclei did not result in the formation of a 
new pattern of etch pits. If 60. were ground off the surface to remove 
the nuclei and the crystal reheated, new dehydration nuclei were formed 
but not in the same position as the previous ones. 

The behaviour of the surfaces on etching and on dehydration is illus- 
trated in figs. 1 to 4. Figure 1 shows the surface after washing with a 
solution of sodium hexametaphosphate, the field containing a row of pits 
which ends abruptly. Figure 2 shows the same field after heating at 98°c 
for several hours. A row of dehydration nuclei have formed and do not 
correspond with the row of etch pits. Figure 3 shows the same field 
after 60 » has been ground off and the surface washed with distilled water. 
The etch pits are fewer in number than in fig. 1, but the rows of pits can 
still be distinguished. 

If dehydration tends to occur at the points of emergence of dislocation 
lines, the facts covering dehydration could be explained on the ground 
that the probability of the formation of a dehydration nucleus at a 
dislocation line at 98°c is low. In this way the dehydration nuclei would 
be unlikely to coincide with etch pits, since the pits mark the site of no 
more than some of the dislocations. It is also possible that the dislocations 
that readily nucleate etch pits are not favourable sites for dehydration. 
This may be due to the type or direction of the dislocation or to the 
segregation of impurities. The facts are also in accord with the view that 
dehydration begins at sites other than dislocations. 

The shapes of dehydration nuclei like those of etch pits are affected by 
the presence of materials adsorbed on the surface. When the surface had 
been washed with water well defined rectangular nuclei were formed, as is 
seen from fig. 4, which shows the field of fig. 3 after heating at 98°c for 
several hours. When the surface had been treated with sodium hexameta- 
phosphate the nuclei were less clearly defined, and were not of the same 
apparent orientation as those formed on water-washed surfaces (see figs. 2 
and 4). The nuclei formed on water-washed surfaces showed the * hour 
glass’ structure observed by previous authors (Zolotov 1958), whereas 
those formed on surfaces treated with sodium hexametaphosphate did 
not. This effect is not very clear in figs. 2 and 4, but was quite distinct 
when lower magnifications were used. ‘The nuclei on surfaces treated 
with sodium hexametaphosphate showed a marked tendency to occur in 


rows. 
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Observation of de Haas—van Alphen Type Oscillation in the 
Acoustomagnetic Attenuation of Zinc+ 


By D. F. Grszons 
Bell Telephone Laboratories, Incorporated, Murray Hill, New Jersey 


[Received August 22, 1960] 


ABSTRACT 


de Haas—van Alphen type oscillations have been observed in the acousto- 
magnetic attenuation of zinc. The oscillations were observed only for 
specific modes and direction of propagation with a field strength greater than 
20-0 k oersteds. Periods in A(1/H) of 2:1+ 0-05, 0-3+0-1 and 7-:0+1-0 x 10-5 
oersteds-+ were measured; these data are in excellent agreement with the 
magnetic data of Verkin and Dmitrenko (1955). 


Tue purpose of this note is to report upon the observation of de Hass—van 
Alphen (dHvA) type oscillations in zinc single crystals in magnetic fields of 
20-70 koersteds. The ultrasonic losses in metals when g/>10, where q is 
the sound wave vector (27/A) and /.the mean free path of the electrons, have 
been shown (Cohen et al. 1960) to originate in the interaction of the acoustic 
wave with the conduction electrons. Any changes in the conductivity 
tensor produced by the application of a magnetic field should in principle 
result in changes in the ultrasonic attenuation. Oscillations in components 
of the conductivity tensor with applied field have been seen to give rise to 
oscillations in the acoustic attenuation, and the maxima in the attenuation 
can be shown to occur where gR=nz for small magnetic fields; R is the 
cyclotron orbit radius. These oscillations are the so-called geometric 
resonances which have been observed by many workers (Morse and Gavenda 
1959, Reneker 1959). The period of these resonances will scale with 4, 
the acoustic wavelength. However, when &<4/2 and the separation 
between the Landau levels iw,, where w, is the cyclotron frequency, is com- 
parable to k7', the semi-classical arguments used to describe geometrical 
resonances breakdown and dHvA type oscillations may occur just as in the 
case of magnetic susceptibility (Shoenberg 1959) and Hall coefficient 
(Bergeran et al. 1960). This assumes that there is sufficient coupling 
between the acoustic field and the electrons for the phenomena to be 
observed. The period of these oscillations is independent of A. This type of 
oscillation in the acoustic attenuation has until now only been observed in 


the semi-metal Bi (Reneker 1959). 
+ Communicated by the Author. 
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Many geometrical resonances have been observed in zinc single crystals 
grown from zone refined material, these crystals have a residual resistance 
ratio Rpp|Ryo_%~ 30x 10? which gives a mean free path ~1-2mm 
(Gibbons 1959). These geometrical oscillations occur in magnetic fields 
~ 100-1000 oersteds. Using a Bitter type solenoid magnet dHVA oscilla- 
tions have been observed in some directions of these crystals in fields 
~20-70k oersteds. The figure shows the data from a zinc crystal for a 
shear wave propagated along (1010) with the particle motion u along 
{0001 ) and the field perpendicular to u (i.e. along (1120) at a temperature 
of 1:0°K. Less pronounced oscillations have also been observed for the 
orientation, q along (1120), u, (0001) and H, (1010). So far,no dHvA 
oscillations have been observed for other directions or modes of propagation 
of the acoustic wave. The temperature dependence of the amplitude of the 
oscillations was measured for the orientation of the figure, and found to be 
0-25 dB/°K between 4:2°K and 1:0°K. At the present time this temperature 
dependence data cannot be used to determine the effective mass of the 
carriers as the coupling coefficient between the acoustic wave and the carriers 
is unknown. 


CHANGE IN ATTENUATION 
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The change in attenuation of a shear wave, propagated along {1010) with 
particle motion along (0001), as a function of 1/H, where H is along 
(1120), at a temperature 1-0°K (27=1-573 cm). 


In the data of the figure the main period exhibited A(1/H) is 
2-1 + 0:05 x 10-* oersteds—1. 


This value is in excellent agreement with the magnetic susceptibility data 
of Verkin and Dmitrenko (1955), (their orientation No. 1 with Gie= D0): 
Below approximately 28koersteds the amplitude of these oscillations 
becomes too small to resolve accurately with the techniques availabe at 
present. In fields greater than 59koersteds a fine structure becomes 
discernible with a period of 0:3 + 0-1 x 10-* oersteds which is in agreement 
with Verkin and Dmitrenko (their curve No. 3); however, they did not 
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observe this period for > 50°, whereas these present data are for 6=90°. 
Both of these components appear to be superimposed upon a considerably 
smaller amplitude low-frequency component with a period 


~7:0+ 1-0 x 10-* oersteds—. 


The small amplitude of this component results in the low precision with 
which this period can be measured. However, this value is in good 
agreement with the value of 6-0 x 10-®oersteds observed by Verkin 
and Dmitrenko. This latter period is the main oscillation seen in their 
work. This branch of carriers, but with @=0° where the period becomes 
60-0 x 10-® oerstedst, has also been seen in Hall coefficient measurements 
(Bergeran et al. 1960). 

Harrison (1960) has attributed the main period observed in this work 
(2-1 x 10-* oersteds) to the horizontal arms (i.e. normal to the C axis) 
of the Fermi surface in the first and second bands of zinc. The geometry 
used in this investigation would be the optimum to observe such segments 
of the Fermi surface. Harrison attributes the smallest period observed to 
the diagonal arms (lying in {1010} planes) of the Fermi surface. The 
geometrical resonances observed for this orientation indicate that the 
wave vector of the carriers in this orientation is k’=0-25x 108cm-!. The 
effective mass is 0-03m,, which would correspond to the horizontal arms. 
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On Dislocation Loops Formed in Zinc Crystals During 
Low Temperature Pyramidal Glide+ 


By P. B. Pricnt 
Physics and Chemistry of Solids, Cavendish Laboratory, Cambridge 


[Received August 2, 1960] 


RECENTLY electron transmission observations have been reported (Price 
1960) of (1122)[1123] pyramidal glide in thin vapour-grown zine crystals 
which were strained in tension inside an electron microscope. The 
orientation of these platelet crystals was such that the dislocations with 
a 3[1123] Burgers vector were of predominantly screw orientation, and 
it was frequently observed that in the path of one of these gliding 
dislocations a long, narrow dislocation loop was left behind on a basal 
plane, presumably as a result of the formation of a multiple jog on the 
screw by cross-glide to enable it to avoid an obstacle. This long loop 
was unstable and was able to lower its energy by splitting up into a row 
of circular loops. It was shown that some of these loops had a 3[{1123] 
Burgers vector and others had a [0001] Burgers vector. Both types 
were sessile but were observed to shrink slowly by climb at room 
temperature and eventually to disappear. 

The intermediate stage of splitting up of a long, narrow loop was 
rapid (1 to 5sec) and it was impossible to take photographs, but it 
appeared that the total area of the circular loops was about the same as 
the area of the elongated loop and therefore that the process took place 
without a net flow of vacancies to or from the loop. The redistribution 
of the area inside the long loop clearly could only take place at a 
temperature at which jogs could easily be formed. It was therefore 
predicted that if the deformation were done at low temperature, this 
process should not be able to occur and the crystal should contain 
elongated loops instead of circular loops. 

The purpose of this note is to describe some electron microscope 
observations, made during straining at low temperature, which confirmed 
this prediction. Figure 1 shows a low temperature micro-tensile device 
which was designed for use inside the Siemens Elmiskop 1. It was similar 
to one designed by Forsyth and Wilson (1960) except that the bimetal 
strips were arranged to stretch the specimen when the temperature was 
lowered instead of raised. The strips, which were made of aluminium 
and molybdenum, were fixed at one end to a phosphor—bronze body 
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which was thermally insulated from the rest of the microscope by a 
teflon collar but electrically connected by a small wire to prevent the 
specimen from charging up during examination. The strips and specimen 
were cooled by bringing into contact with the body of the device a 
thermally insulated finger, which was connected through the microscope 
column to a low temperature bath. By changing the level of liquid 
nitrogen in the bath the temperature of the bimetal strips, as determined 
by a contact thermocouple, could be lowered from 20°c to —150°c and 
the tensile strain increased to a maximum of 8%. A specimen was 
always mounted in a sagging position so that the slack was not taken up 
and deformation started until it was considerably below room temperature. 


TEFLON 
COLLAR 


TO COLD 
FINGER 


Low temperature straining stage for Siemens Elmiskop 1. 


As was expected, it was found that loops produced at low temperature 
were elongated. Figure 2+ is a series of six transmission electron 
micrographs of the same set of loops taken at intervals of about 20 sec 
when the recorded temperature was —60°c. As long as the electron 


} Figure 2 is shown as a plate. 
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beam was off, the loops were stable and showed no tendency to split 
up at this temperature. During observation, the local temperature was 
raised enough (probably by at least 10 to 20°c) to allow the middle one 
of the three long loops to split into a row of circular loops. By proper 
adjustment of the beam current, the process could be made to occur 
slowly enough to photograph, and fig. 2 shows the various stages in the 
formation of the second circular loop in a row which eventually comprised 
six loops. These stages are exactly the same as those which occurred at 
room temperature (but at a different rate) and which are depicted in the 
drawing, fig. 8, of the paper by Price (1960). The elongated loops to 
either side of the centre one (fig. 2) did not change their shape, presumably 
because they were not in the centre of the illuminated region and were 
therefore at a lower temperature. 

It was verified by direct measurement that the total area of the six 
circular loops was the same, to within 10%, as that of the parent loop 
and therefore that there was no net flow of vacancies to or from the 
loop. Measurements also showed that the criterion for loop splitting 
(Price 1960) was satisfied and that the energy, as calculated from the 
geometry of the loops, was lowered by ~ 25% by splitting. 

Further studies of low temperature deformation are in progress ; the 
results will be reported later. 
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The Cooperative Movement of Dislocations in Very Pure 
Aluminium Foil 


By R. N. Witson and P. J. E. Forsyrsa 


Metallurgy and Physics Department, Royal Aircraft Establishment, 
Farnborough, Hants 


[Received September 29, 1960] 


DisLocatTion sub-boundaries in aluminium have been shown to be 
essentially of two types: (a) tilt boundaries consisting of arrays of edge 
dislocations formed by thermally activated dislocation climb and elastic 
interaction with like dislocations on parallel planes, and (b) twist 
boundaries formed by the cross slipping and arraying of screw dislocations. 

Tilt boundaries in aluminium have been studied by Cahn (1949) who 
showed they could be produced in a single crystal by light deformation 
followed by an annealing treatment. He called this phenomenon poly- 
gonization. Hirsch ef al. (1956) and Hirsch et al. (1958) have studied 
dislocation sub-boundaries formed by heavy cold-working. They have 
used transmission electron metallography for direct observation and have 
shown many of these boundaries to be chiefly composed of screw disloca- 
tions. 

The authors have studied the movement of dislocations in foils of 
aluminium using transmission electron microscopy with a straining 
device described elsewhere (Forsyth and Wilson 1960). As part of this. 
work, a study was made of the influence of impurity level upon dislocation 
motion and the mode of fracture of the foils. Thin foils were prepared 
from super purity aluminium (99-99°% Al) and zone refined aluminium 
(99:9999% Al). It was found that whereas deformation in the super 
pure aluminium was by the independent movement of dislocations, those 
in the zone refined metal moved in a cooperative manner as shown in the 
plate. The axis of tension in the foil is vertical in the micrograph and 
dislocations, believed to be essentially screw in character, can be seen to 
have been moving from right to left before the photograph was taken. 
It can be seen that many dislocations are in the process of moving in 
formation as well defined arrays of several dislocations, i.e. they are not 
only being influenced by the repulsive forces from like dislocations on 
their slip planes but also by the attractive forces of dislocations on parallel 
slip planes. Wilsdorf (1959), discussing the frictional force influencing 
dislocation movement in thin foils of stainless steel, has suggested two 
factors may be operative. The first is the presence of a surface oxide 
film anchoring the ends of the dislocation: This could account for the 
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convex bowing of dislocations under stress. The second factor is the 
pinning of a dislocation by atmospheres of point defects or impurity 
atoms along its length. This he suggests could account for the great 
variation in width of stacking fault in stainless steel and also be one cause 
of the jerky dislocations observed in metals of higher stacking fault energy 
suchasaluminium. These frictional forces may also retard one dislocation 
with respect to its neighbours due to the variation in degree of locking 
of the dislocations as they move through the lattice. 

In aluminium of very high purity the latter frictional force might be 
expected to be low. Also the absence from solid solution of such foreign 
atoms as magnesium and copper, which have a tendency to associate with 
vacancies, might reduce the number of vacancies present and increase the 
mobility of those remaining. It is suggested that it is this reduction of 
lattice friction which increases the effectiveness of the attractive forces 
between like dislocations on neighbouring planes causing this cooperative 
movement of dislocations. By this process small angle arrays can move 
unbroken through the lattice. 
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Correlation between Magnetic Susceptibility and Hydrogen 
Solubility in Alloys of Early Transition Elements 


By D. W. Jonzs, N. Pessaut and A. D. McQuitLan 
Department of Physical Metallurgy, University of Birmingham 


[Received December 1, 1960] 


Tue elements of groups IVA and VA of the periodic table share with 
palladium the ability to take into solid-solution large quantities of 
hydrogen, but whereas in the case of palladium this type of behaviour 
can be explained in terms of an almost filled d-band, no explanation has 
yet been offered for the case of the early transition elements. As part of 
a general investigation of the physical properties of alloys of the early 
transition elements with each other, it was decided, therefore, to include 
hydrogen solubility as one of the properties to be examined. In the 
preliminary survey the elements chosen for study were titanium, niobium 
and molybdenum because these metals are from successive groups of the 
periodic table, and because they are completely miscible with each other 
in the body-centred cubic structure. It was possible, therefore, to follow 
continuously changes in physical properties in passing from one metal 
to the next. Since molybdenum, being in Group VIA, does not dissolve 
hydrogen, the systems containing molybdenum are of particular interest 
because in these systems hydrogen solubility must fall to zero with 
increasing molybdenum content. 

Figure 1 shows the maximum amount of hydrogen which can be 
absorbed by alloys in the binary systems, titanium—niobium, and 
titanium—molybdenum. The alloy compositions have been plotted in 
terms of electron to atom ratio assuming the values 4, 5 and 6 for titanium, 
niobium and molybdenum respectively. With increasing electron/atom 
ratio, the capacity of the alloys to absorb hydrogen decreases continuously, 
and is effectively zero at about 5-6 electrons per atom. Because titanium- 
rich alloys form hydride phases with face-centred cubic structures, the 
crystal structure does not remain the same over the length of the curve 
shown in fig. 1. It is preferable, therefore, to measure not the total 
solubility of hydrogen in the alloys, but the heat of solution of a very small 
quantity of hydrogen in the alloys having a body-centred cubic form. 
Figure 2 shows the results of this type of investigation. As the electron/ 
atom ratio approaches 5-6, the heat of solution decreases rapidly, passing 
through zero at 5-5 electrons per atom, and thereafter becoming endo- 
thermic. The zero value of heat of solution has little significance, being 
merely the point where the reaction with the metal changes from 
exothermic to endothermic with respect to gaseous molecular hydrogen. 
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The magnetic susceptibility at room temperature of alloys in the three 
binary systems titanium—niobium, niobium—molybdenum and titanium— 
molybdenum are plotted in fig. 3. The most characteristic features of 
the curves are the minimum points which occur in the systems titanium— 
molybdenum, and niobium—molybdenum at electron/atom ratios between. 


Fig. 1 
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Total amount of hydrogen absorbed by alloys of titanium with niobium 4A, 
and titanium with molybdenum @, as a function of alloy composition 
expressed in terms of electron/atom ratio. 


5-6 and 5-7 and which correlate very well with the disappearance of 
hydrogen solubility in the alloys. There exists, furthermore, one other 
piece of evidence that a change in electronic structure occurs at this 
critical composition. Figure 4 gives results (Grum-Grzhimailo and 
Gromova 1956) on the Hall coefficient of body-centred cubic titanium— 
molybdenum alloys at room temperature. The Hall coefficient is positive 
throughout and reaches a maximum value at about 5-6 electrons/atom. 
The magnetic susceptibility and Hall coefficient curves taken together 
strongly suggest that, as the critical composition of 5-6 is approached, a 
sub-band is being filled, since the Pauli paramagnetism of the alloys, 
which is proportional to the density of states at the Fermi level is decreasing 
rapidly immediately before the critical composition is reached, and the 
Hall coefficient, as would be expected if the top of a band were being 
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Fig. 2 
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Heat of solution of molecular hydrogen in body-centred cubic alloys in sthe 
systems titanium—niobium 4A, and titanium—molybdenum @ in terms 
of electron/atom ratio. 
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- Magnetic susceptibility of body-centred cubic alloys of titanium—niobium A, 
niobium—molybdenum , and titanium—molybdenum @ at room temperature. 
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approached, is becoming increasingly positive in the same composition 
range. A similar picture of a sub-band filling at an electron/atom ratio 
of approximately 5-6 electrons emerges from results obtained from the 
measurement of the low-temperature specific heats of body-centred 
cubic alloys of the early transition metals of the first long period (Cheng 
et al. 1960). 


Fig. 4 
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Hall coefficient of body-centred cubic alloys of titanium—molybdenum @ at 
room temperature, after Grum-Grzhimailo and Gromova 1956. 


The form of the curves of magnetic susceptibility and Hall coefficient 
shown in figs. 3 and 4, together with the results of Cheng et al. strongly 
suggest that, beyond the critical composition of 5-6 electrons per atom, 
there exists another sub-band, and that the top of the lower sub-band 
overlaps the bottom of this upper band. In a recent paper (Goodenough 
1960), it is proposed that early transition metals with the body-centred 
cubic structure have a band structure consisting of a partially filled 
s-p bonding band overlapping a sub-band composed of d-electrons in 
f,, States. On this model, the alloy having the critical composition of 
5-6 electrons per atom would have three electrons per atom in the is 
sub-band and 2-6 electrons per atom in the s-p bonding band. Further- 
more, an increase in electron/atom ratio would lead to the filling of 
localized e, states. While there are difficulties in reconciling the details 
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of Goodenough’s model with our results, it is evident that this type of 
approach, in which the d-band of the early transition metals is split 
into sub-bands, is likely to be extremely productive in attempting to find 
an explanation for the unusual physical properties of these materials. 
The question will be discussed in greater detail in later papers when the 
results of further experiments now in progress become available. 
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Density and Expansivity of Solid Xenon 


By A. J. Eatwewu and B. L. Smrrx 


Department of Physics, Queen Mary College (University of London), 
London, E.1 


[Received October 17, 1960] 


Ty previous communications, measurements of the densities and expan- 
sivities of solid argon (Dobbs et al. 1956) and solid krypton (Figgins and 
Smith 1960) have been reported. We have now determined the lattice 
parameter of solid xenon between 20° and 120°k using the x-ray method of 
Debye-Scherrer powder photographs and a bulk density method and 
have thus derived values of the density and the volume coefficient of 
expansion. 

In the x-ray measurements, polycrystalline xenon was deposited from 
the vapour in a Perspex tube of 1 mm internal diameter, and all specimens 
used were free from visible defects. The temperature was kept constant 
to within 0-04°K at 20°K and to 0-01°K at 90°k. Thus temperature 
fluctuations could only contribute a negligible error (5x 10-7 A) in the 
values of lattice parameter. ‘ Spottiness ’ in the films due to large grain 
size at the higher temperatures has been much reduced by the introduction 
of vertical and rotational motion of the cylindrical specimen. Thus, the 
accuracy quoted for x-ray values at 20°K by Figgins and Smith (+ 0-001 A) 
now applies to the whole range of working temperatures. The lattice 
parameter values are not corrected for refraction, this correction being 
small (Jette and Foote 1935). The xenon used in both experiments was 
stated by the British Oxygen Company Ltd to be of 99-100 purity, 
the remainder being krypton. The estimates of lattice parameter made 
from the bulk density measurements are accurate to within + 0-005 A. 

The results are shown in the figure. It will be seen that there is agree- 
ment between the results obtained by the two methods. Our results are 
in fairly good agreement with the x-ray measurements of Cheesman and 
Soane (1957), are within the limits of the experimental error of the single 
value of Simon and Ruhemann (1932) but outside those of the single value 
of Natta and Nasini (1930). It is of interest to note that a separate 
sample of gas, also stated to be of 99-100% purity, was used for the x-ray 
measurement at 60°K. 

In the tables, we summarize our results for the density and estimates 
of the volume coefficient of expansion «=1/V (dV/0T), at rounded 
temperatures. For convenience the corresponding results previously 
obtained for argon and krypton are repeated. 
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Table 1. Density (g cm) 


ae I) 20 40 60 80 100 120 


Argon 1-764 List 1-691 1-636 a= = 
Krypton 3-078 3-040 2-988 2-926 — 
Xenon 3°772 3-736 3-689 | 3-642 3-593 


3-542 


80 


18 
11-0 
6-65 6-77 


An extrapolation of our curve leads to an estimate at absolute zero of 
a lattice parameter of 6-13 A for xenon and a corresponding density of 
3-78 g cm->. 
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Dislocation Loops and Hardening in Neutron Irradiated Copper 


By M. J. Maxty, A. D. WuHaruam and F. J. Minter 
Metallurgy Division, A.E.R.E., Harwell 


[Received October 17, 1960] 


Silcox and Hirsch (1959) using the transmission electron microscope 
technique have shown that in copper prismatic dislocation loops are 
formed during neutron irradiation at 35-100°c. The loops increased in 
diameter with dose and their density could be substantially reduced by 
annealing for a few minutes at 350°c. They suggested that the loops 
were vacancy loops similar to those observed in quenched metals, but 
that they were nucleated by the collapse of vacancy-rich regions in the 
centre of Brinkman displacement spikes. They further concluded that 
the lattice hardening formed during irradiation was the stress required 
for a dislocation to cut through the ‘ forest’ of loops threading the glide 
plane. In considering the details of this model certain objections become 
apparent, notably the difficulty in simultaneously accounting for the 
very large temperature dependence and the absolute magnitude of the 
lattice hardening. For these reasons experiments were undertaken to 
investigate (a) the nature of the loops and (b) the relation between the 
loops and the irradiation hardening. 

Foils, 0-001 in. thick, of 99:998°% copper were irradiated at ~27°c 
in an air-cooled hollow uranium cylinder in the BEPO reactor at Harwell, 
and were subsequently electropolished in an orthophosphoric acid—water 
mixture and examined in a Siemens Elmiskop I microscope. The 
distribution of the number of loops as a function of diameter after 
various neutron doses is shown in fig. 1. The average diameter of the 
loops increased during irradiation, showing that individual loops must 
grow by acquiring mobile defects. Loops can grow therefore only in the 
temperature range in- which the constituent point defects are mobile. 
In principle it is possible to distinguish whether the loops are formed 
from platelets of interstitial atoms or vacancies by irradiating at a 
sufficiently low temperature to immobilize one of these defects. At 27°c 
it is believed that both interstitial atoms and vacancies are mobile, 
-assuming that the energies of migration are ~0-l ev for interstitials 
and ~0:7ev for vacancies (Holmes et al. 1958). Irradiations were 
therefore carried out at —195°c and the specimens examined on the 

-eold' stage of a Siemens microscope keeping the temperature low 
throughout. Two procedures were adopted. In the first, bulk foils 
' 0/001 in. thick were neutron irradiated (3-1 x 10" fast n cm~?) at — 195°c 
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and then polished in a methyl alcohol-nitric acid bath at — 70°c. They 
were washed in methyl alcohol and then ether at —78°c and dried 
in vacuo at —50°c. Transfer into the microscope specimen holder was 
performed under liquid nitrogen. In the second procedure previously 
polished specimens were irradiated at —195°c with pile neutrons or 
with 100 kev alpha particles in an accelerator, care being taken to ensure 
good thermal contact between the specimen and the cooled support. 


Fig. 1 
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The density of dislocation loops as a function of diameter for neutron doses 
up to 2-4 x 1018 fast n em-, 


As before these specimens were inserted into the cold stage at —195°c 
This technique eliminates the polishing and washing at —70°c and he 
drying at —50°o entailed in the first procedure and enables the specimens 
to be kept within a few degrees of —195°c until they are inside the 
microscope. The cold stage operates at ~ — 120°c but at this temperature 
ice condenses from the microscope vacuum onto the specimen. This 
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can be removed by the electron beam at a local temperature of — 50°C 
and areas could be photographed within one minute of warming to this 
temperature. In all cases loops were observed immediately the ice 
cleared (fig. 2+) and no further change occurred even on warming to 
room temperature. This result virtually eliminates the possibility that 
the loops are formed from vacancies since these are capable of only 
0-1 Jumps per minute at —50°c assuming that the energy for migration 
is 0-7ev. Hence loops are formed below the temperature at which 
vacancies are mobile and it is inferred therefore that the loops must be 
formed from interstitial atoms. This conclusion is supported by the 
observation of Barnes and Mazey (1960) of the growth of the loops in 
the helium band of alpha irradiated copper on annealing. In this case 
it is concluded that the loops expand by emitting vacancies (which are 
absorbed by the helium bubbles) and must therefore be interstitial in 
character. 

The role that the interstitial loops play in the irradiation hardening 
is shown by experiments in which electron micrographs are compared 
with the mechanical properties during annealing at 306°c. A large 
number of polycrystalline tensile test specimens and 0-001 in. foils were 
irradiated with 2 x 1018 fast nem? and annealed at 306°c for various 
times. The electron micrographs from the annealed foils show two 
effects : (a) the density of loops decreases rapidly with annealing and 
(0) a high density (1:3 x 1016 em~$) of very small (<25 A) defects become 
readily visible (fig. 3). That these defects are different in kind from the 
loops is indicated by the very large difference in size and also the different 
response to the grain boundaries. The loops form a wide denuded zone 
on annealing but the small defects show no appreciable denudation 
(fig. 4). The small defects are thought to be vacancy clusters formed 
by the redistribution of the multiple vacancies created in the displacement 
spikes (Seeger 1958). If the defects observed in the electron microscope 
are responsible for the hardening the yield stress should be proportional 
to the square root of the loop density. A plot of the square root of the 
densities of both loops and small defects versus the yield stress after 
various anneals (fig. 5) shows that whereas there is no linear relationship 
between the hardening and the number of large loops, there may well 
be such a relation between the small defects and the hardening. This 
suggests that the small defects may be responsible for a large part of 
the hardening. This is supported by the observation that on mild 
annealing (i.e. 10-60 min at 306°C) simultaneously with the decrease 
in density of the large loops and the formation of the small defects there 
is a large reduction in the temperature sensitivity of the hardening. 
It is unlikely that this is associated with the decrease in density of the 
large loops since theoretically these are expected to produce a hardening 
with a low temperature dependence. It is more reasonable to associate 


+ Figures 2, 3 and 4 are shown as plates. 
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the temperature dependence with submicroscopic defects which anneal 
away to produce the observable small defects. Furthermore, nickel 
does not show any large loops on irradiation at 35°C (Silcox 1960) although 
the great similarity in the characteristics of the hardening in this metal 
to those in copper (Makin and Minter 1960) suggests that the same types 
of defects are responsible in both. 
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It is concluded that in copper although the interstitial loops must 
play some part in the hardening they are unlikely to be the sole cause. 
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REVIEWS OF BOOKS 


Translation from Russian for Scientists. By C. R. Buxton and H. SHenpon 
Jackson. (Blackie & Son, 1960.) [Pp. xvii+299.] 30s. 

A concise survey of Russian grammar is followed by an extensive series of 

texts drawn from physics, chemistry and engineering. The book provides 

a valuable introduction to scientific Russian and should enable scientific 

journals to be read without difficulty. MaGye: 


Magnetic Materials. Third edition. By F. Brattsrorp. (London : Methuen 
& Co; New York : John Wiley & Sons, 1960). [Pp. vii+188.] 16s. net. 


THE study of magnetism and magnetic materials has expanded so enormously 
since the publication in 1948 of the first edition of this well-known Methuen 
Monograph that at a recent conference no less than 164 papers on these topics 
were presented. This new edition nevertheless encompasses many of .the 
advances of the past 12 years in a volume not exceeding unduly that of the earlier 
one. It should hence continue to be a useful companion to those engaged in the 
art and practice of magnetic materials. The major addition is a chapter on 
ferrites, both cubic and’ hexagonal, appearing here with chapters on iron- 
silicon, iron—nickel and related alloys as well as on permanent magnet materials. 
One looks forward to the fourth edition, where, it is hoped, an equally coherent 
account will be given of the properties of thin films for computer applications, of 
garnets and of the novel permanent magnet materials, as well as of other sub- 
stances as yet undreamt of. One despairs, however, of magnetic recording tapes 
being regarded as respectable enough to be included even in the nth edition. 
E. P. W. 


Cartesian Tensors. By G. Tempte. (London: Methuen & Co. Ltd., 1960.) 
[Pp. 92.] 12s. 6d. 


Every physicist needs to know just a little about tensors, and it is not easy to 
find a good simple account of the elements. Professor Temple has attempted 
this; for its very existence his book is welcome. Yet it is really a trifle dry— 
much more suitable for applied mathematicians than ordinary physicists. It is 
a pity that he has concentrated so much on the representation of vectors and 
tensors by their components and ignored the advantages, at least for some ways 
of thought, of conceiving them as geometrical objects which can be symbolized 
as a whole, in Clarendon Bold type as you might say. The whole theory of the 
eigenvalues, for example can be flooded with light by the simple analogy of the 
reduction of a quadric surface to its principal axes. The chapter on spinors 
also seems very formal and difficult for a beginner. However, it will be useful 
to know where to refer students for an elegant deviation of the main theorems and 
their application. ice J.M. Z.. 


Progress in Solid Mechanics. Edited by I. N. Syeppon and R. Hint. (The 
North-Holland Publishing Company, 1960,) [Pp. xii+448.] Price £5 0s. 0d. 


Tus book is Volume I of a new series of which the aim is to present “ review 
articles of the highest standard which from time to time summarize and unify 
the most recent work in a particular field or group of fields” in the mechanics 
of solid continua. It contains the following articles: Viscoelastic waves (8. C. 
Hunter), Matrices of transmission in beam problems (K. Marguerre), Dynamic 
expansion of spherical cavities in metals (H. G. Hopkins), General theorems 
for elastic-plastic solids (W. T. Koiter), Dispersion relations for elastic waves 
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in bars (W. A. Green), Thermoelasticity: the dynamical theory (P. Chadwick), 
Continuous distributions of dislocations (B. A. Bilby), Asymmetric problems 
of the theory of elasticity (R. Muki). 

Each writer reviews the work done on his topic up to 1958 and gives a useful 
list of references. In general the standard of exposition is high, being concise 
without assuming too much previous knowledge on the part of the reader. 
The layout of the mathematics is pleasing and easy to read. 

This series promises to be of very great value to workers on the mechanics of 
solids: the North-Holland Publishing Company is to be congratulated on its 
initiative, and the editors thanked for the excellent start they have made on 
the task they have set themselves. K. R. L. 


Mechanics. Second edition. By K. R. Symon. (Addison-Wesley, 1960.) 
[Pp. 557.] Price £3 Os. Od. 


Tuts text-book can be warmly recommended. It is a treatment of classical 
mechanics with a choice of topics intended to emphasize present-day points 
of view. The inadequacies of classical mechanics are pointed out, as are the 
points of connection with relativity and wave mechanics. Only a knowledge 
of the calculus is assumed; there are adequate accounts of other mathematical 
techniques such as vector analysis and tensors. The first eight chapters deal 
with such topics as Newtonian concepts, motion of a particle in three dimensions, 
rotation about a fixed axis, gravitation, moving coordinate systems and wave 
motion. The treatment is on the whole too sophisticated for students attending 
a first University course, although the brighter members of a class would find it 
very stimulating. The last four chapters deal with a Lagrangian mechanics, 
inertia and stress tensors, rotation about a point, and small vibrations of a 
system. The book in fact covers all aspects of classical mechanics for which 
space is likely to be found in undergraduate courses up to and including the 
Honours level. A student who, with some guidance, worked through even an 
appreciable part of this book, would have laid a very solid foundation for study 
in other branches of physics. There are approximately 360 problems, and 
answers are given to half of these. The book is clearly written and produced, 
and for its length it is very reasonably priced. Wea ue 


Crystal-Structure Analysis. By M. J. Burreer. (New York: John Wiley, 
Inc., 1960.) [Pp. 668.] Price £7 8s. Od. 


Tus book by a distinguished American crystallographer will be in continual 
use in all laboratories concerned with crystal-structure analysis. It contains a 
wealth of material, set out clearly at considerable length, and there are many 
references for further reading. The present volume is intended as a companion 
volume to the author’s X-ray Crystallography and Vector Space and its appli- 
cations in Crystal Structure Analysis. The former dealt largely with experi- 
mental methods; to some extent the ground is gone over again in the present 
volume, which does not however overlap to any extent with ‘ Vector Space ’. 
The book opens with an interesting historical survey, and a discussion of the 
theory of diffraction by crystals. Stibsequent chapters deal with, for example 
the measurement and correction of X-ray intensities, examples of crystal 
structure determinations, the theory of Fourier methods, the reciprocal lattice 
methods of phase determination, and accurate structure analysis. The very 
complete treatment of most topics serves to emphasize some surprising omissions. 
The author’s stated reason for omitting optical analogue methods may or may 
not be sound, but the fact that intensity statistics is a self-contained topic on 
which a monograph could be written does not justify its omission from a com- 
prehensive textbook. Most surprising of all for a textbook which, according 
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to the publishers, sets out to review all the available methods for crystal 
structure analysis, the word ‘neutron’ does not even appear in the index. 
X-ray crystallographers should not be encouraged in their insularity by ignoring 
the fact that there is now a large and increasing effort devoted to structure 
determination, particularly of magnetic materials, by neutron diffraction, 
This is a good book; it could so easily have been a better one. W.C. 


Mechanical Waveguides. By Martin Repwoop. (London and New York: 
Pergamon Press, 1960.) [Pp. 300.] Price £2 10s. Od. 


Many measurements of the velocity and attenuation of ultrasonic waves in 
fluids and solids have been made in recent years and this book is primarily 
addressed to researches in this field. It is a mathematical text on elastic wave 
propagation in lossless media and does not discuss experimental techniques. 
But it is a timely addition to the literature on elasticity and brings together a 
large number of results from many different courses. 

A useful introduction to the wave equation precedes a detailed account of 
the reflection of elastic waves at an.interface. These interfaces are the boun- 
daries of the ‘ mechanical waveguides’ through which the elastic waves are 
propagated. A detailed analysis of the propagation of both continuous waves 
and pulses in fluid and solid waveguides is given. In many cases the sine and 
cosine solutions of the parallel plate guide are followed by the Bessel function 
solutions for the cylindrical guide. The analysis refers only to lossless, isotropic 
media, but includes both the ultrasonic ‘ narrow bandwidth’ pulses and the 
explosive ‘ wide bandwidth ’ pulses studied extensively by Kolsky. 

Unfortunately it is not an easy book to read. The text is broken up un- 
necessarily by the display of the longer equations over many lines and by the 
inclusion of over thirty pages of * additional ’ references, which are not referred 
to in the text. Some chapters of the book are subdivided ad infinitum, so that 
for example, § 9.2.5 (iii) contains only three sentences. But the author is to be 
congratulated on including so many different examples of elastic wave 
propagation in one short book. E. R. D. 


Infra-red Methods: Principles and Applications. By G. K. T. Coun and 
D. G. Avery. (New York and London: Academic Press, 1960.) [Pp. 203.] 
Price $6.80. 

Tuts book provides a clear and concise introduction to current infra-red 

techniques and is intended for the research student or worker in allied fields. 

The first part covers a wide range of topics and has chapters on sources, optical 

materials, detectors, amplifiers and dispersive systems. The emphasis here is 

on the underlying theory which is treated concisely and for the most part 
adequately. The shorter section deals briefly with a few specialized applications 
and considers calibration of detectors, a prism monochromator, gas analysis 
and optical pyrometry. The book is well supplied with graphs, diagrams and 
useful tables, and since the treatment is general, rather than detailed, the 
authors wisely include at the end of each chapter an excellent list of references. 

The only criticism worth mentioning is that the much greater space given to 
the description of prism spectrometers might suggest that they are technically 
more important than grating instruments. Since the diffraction grating is 
rapidly superseding the prism and is, as the authors point out, a cheaper and 

better dispersive element, the grating spectrometer deserves, if anything, a 

fuller treatment than the prism instrument. Criticisms other than this would 

be of minor faults. Infra-red Methods is not only an adequate introduction to 
the subject but also a compact reference book. es 
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A Degree Physics. Part 1: The General Properties of Matter. By C. J. SMirx. 
(London: Edward Arnold, Ltd., 1960.) Second edition. [Pp. 732.] 
£3 3s. Od. 


Molecular Distillation. By G. Burrows. (Monographs on the Physics and 
Chemistry of Materials ; General Editors: W. Jackson, H. FROHLICH, 
N. F. Mort, E. C. Butuarp.) (Oxford : Clarendon Press, 1960.) [Pp. 214.] 
£1 15s. Od. 


Lectures in Theoretical Physics. Vol. U1. Lectures delivered at the Summer 
‘Institute for Theoretical Physics, University of Colorado, Boulder, 1959. 
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tors, Part 2. Proceedings of an International Conference held in Brussels, 
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and New York: Academic Press, 1960.) [Pp. 639.] Price £6 8s. 6d. 
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(London and New York : Academic Press, 1960.) [Pp.557.] Price £6 8s. 6d. 


Growth of Crystals, Vol. 2. Interim reports between the first (1956) and second 
Conference on Crystal Growth, Institute of Crystallography, Academy of 
Sciences, USSR. Edited by A. V. Stupntkov and N. N. SHerrau’. Trans- 
lated from the Russian. (London: Chapman & Hall Ltd.; New York : 
Consultants Bureau, Inc., 1959.) [Pp. 178.] Price £6 8s. Od. 
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K. U. SNOWDEN Phil. Mag. Ser. 8, Vol. 6, Pl. 31. 


Fig. 4 


<10* cycles after repolishing at 1-75 ~x 10° 5 x 104 cycles ( x 75). 
cycles. Fine slip traces produced the 
mottled effect. Arrows indicate boun- 
daries at +45 degrees to the specimen axis 
which is parallel to the long edge of the 
figure. Strain +0-092% (x75). 


Fig. 5 Fig. 6 


5-30 x 10° cycles, repolished 5-93 x 10° cycles, repolished 
7-5 x 10° cycles ( x 75). 5-30 x 106 cycles ( x 75). 


n figs. 4-8 the specimen axis is parallel to the long edge of the figure. In figs. 4-6 
the two main systems of slip are indicated by double-headed arrows, 
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Fig. 7 Fig. 8 


8:0 x 10° cycles ( x 75). 8-50 x 108 cycles, repolished 
at 8-20 x 10® cycles (x 75). 


Fig. 10 


After 2 x 104 cycles (x 75). 
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bi 


| 
a 


Untreated thin crystal of nC,,H,,. Transmission electron micrograph 
showing Bragg fringes. x 1760. 


Fig. 2 


nC3,H 74 crystal with a striations. Photomicrograph. x 630. 
Fig. 3 


nOggH 74 crystal with b and <110> striations. Photomicrograph. x 630. 


A. KELLER Phil. Mag. Ser. 8, Vol. 6, PI. 34. 


nCzgH 74 crystal with a striations. Electron micrograph. 
Replica. Au—Pd shadowed at ~26°. x 3600. 


Fig. 5 


Detail of a NC sgH a4 crystal with a striations. Electron micrograph. 
Replica, Au-Pd shadowed at 9°. x 5500. 


A. KELLER Phil. Mag. Ser. 8, Vol. 6, PI. 35. 


Fig. 6 


Detail of a nC3,H,, crystal with a (NE-SW) and b (NW-SE) striations. 
Electron micrograph. Replica, Au-Pd shadowed at 12°. x 4300. 


Fig. 7 


Detail of a nCy,H,, crystal with a (vertical) and b (horizontal) striations. 
Electron micrograph. Replica, Au—Pd shadowed at ~26°. x 8900. 


A. KELLER Phil. Mag. Ser. 8, Vol. 6, Pl. 36. 


Fig. 8 


Central detail of a nC3,H,, crystal with b and <110> striation. 
Electron micrograph. Replica, Au—Pd shadowed at 10°. x 4600. 


Crystal of a nC3)H,. with <130> and <110> striation. Electron micrograph. 
Replica, Au—Pd shadowed at 10°. x 4600. 


Fig. 10 


Crystals of nC;,H,, with a striation. Transmission electron micrograph 
showing Bragg extinction contrast. Unshadowed, x 840. 


A. KELLER Phil. Mag. Ser. 8, Vol. 6, Pl. 37. 
Fig. 11 


Thin crystals of nC3,H,, with cracks along the a striations. Transmission 
electron micrograph taken in a hot stage. x 1630. 


Fig. 13 


Detail of a crystal of nCsgH4, Crystal of nOz,H,, with b striations 
with pleat-like b striations. showing in Bragg extinction 
Transmission electron micro- contrast. Transmission electron 


graph. x 1960. micrograph. x 1850. 


D. GC, BASSETT and A. KELLER Phil. Mag. Ser. 8, Vol. 6, Pl. 38. 


Polyethylene lozenge with Bragg fringes dividing the crystal into four quadrants. 
Notice the different contrast in opposite pairs of quadrants. Electron- 
micrograph. (x 1800). 


Polyethylene lozenge with rectangular bands shown up by diffraction 
contrast. Electronmicrograph. (x 2500.) 


D. C. BASSETT and A. KELLER Phil. Mag. Ser. 8, Vol. 6, PI. 39. 


A central detail of a polyethylene crystal as in fig. 2, showing ridge and pleat 
striations. Electronmicrograph. Au—Pd shadowed. (x 9000.) 


Fig. 4 


Detail of a crystal as in figs. 2 and 3 showing a regular ridge in a polyethylene 
crystal. Electronmicrograph. Au—Pd shadowed. (x19 000.) 
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Polyethylene crystal with striations on the upper layers only. 
Electronmicrograph. Au—Pd shadowed. 


( x 8500.) 


D. C. BASSETT and A. KELLER Phil. Mag. Ser. 8, Vol. 6, Pl. 41. 


Fig. 6 


Truncated polyethylene lozenge with b ridges in the {100} sectors. Notice the 
{100} sector boundaries. Electronmicrograph. Au—Pd_ shadowed. 
( x 8000.) 


Detail of a truncated polyethylene lozenge with regular ridges along b in the 
{100} sectors. The overgrowth consists of crystals formed at a lower 
temperature. Electronmicrograph. Au—Pd shadowed. (x 23 000.) 


—— 


D. C, BASSETT and A. KELLER Phil. Mag. Ser. 8, Vol. 6, PI. 42. 


Truncated polyethylene lozenge with two {110} sectors in diffraction constrast. 
Notice the central triangular fold which is a characteristic feature of this 
type of lozenge. Electronmicrograph. (x 2800.) 


Fig. 9 


A non-planar polyethylene lozenge prepared so as to avoid a collapse on drying 
down (see text). Electronmicrograph. Au—Pd shadowed, ( x 6500.) 


D. C. BASSETT and A. KELLER Phil. Mag. Ser. 8, Vol. 6, Pl. 43. 


Well-defined fringes in a polyoxymethylene crystal. The contrast is provided 
mainly by diffraction. Electronmicrograph. (x 4500.) 


Fig. 11 


Detail of a crystal as in fig. 10, revealing the underlying morphology 
causing the fringes. Electronmicrograph. (x 4000.) 
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Sodium L,, freshly evaporated on to copper. Broken line is for /(#)a #1. 
Inset is theoretical N(#) vs EH curve according to Jones and Mott (1937). 
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Boron K spectrum from amorphous boron. Plotted points assume J(#)oo#1/2, 
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Silicon L,, from powdered AR grade Si at red heat. Points plotted for 
L(E)a Ht?, 
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(a) Decrease of intensity with time for peak of Li K band. 
(b) Change of intensity with time at 253 A (peak of Li Line). 
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Fig. 4 (continued) 
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(d) 


(c) Li K band and Li Line, the latter developed to maximum intensity. 
(d) Li K band with no Li line following heavy evaporation of Li. 
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Microstructure after fatiguing showing slip in the ferrite. 
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Crack in specimen fatigued in vacuum below the fatigue limit. 
The crack was opened by subsequent extension of the specimen. 


Crack in specimen fatigued in air below the fatigue limit. 
The crack was opened by subsequent extension of the specimen. 
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Triangular etch pits in positive orientation on an octahedral face after etching 
in air at 850°C. = x175. 


Fig. 2 


Triangular etch pits in negative orientation on an octahedral face after etching 
in air at 1250°o, x 470. 
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Fig. 3 


Round pits on octahedral face after heating in air at 950°C. x 875. 


Fig. 4 


Square pits on cubic face after heating in air at 775°C. x 1500. 
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Boat-shaped. pits on dodecahedral face after heating in air at 800°c. x 470. 
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Boat-shaped pit on dodecahedral face after heating in air at 1400°c. x 1500. 


Triangular etch pits on octahedral face in positive orientation after heating in 
nitrogen, oxygen and water vapour mixture at 850°C. x 175. 
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Matt type surface on octahedral face after heating in dry oxygen at 
850°o. 175. 


Fig. 11 


Triangular etch pits on octahedral face in negative orientation after heating in 
nitrogen, oxygen and water vapour mixture at 1250°c. x 470. 
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Etch pits formed on the (010) face of a crystal of gypsum after treatment for 
3 min at 20°o with 1% solution of sodium hexametaphosphate. x 330. 
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J. E. BRIGHT and M. J. RIDGE Phil. Mag. Ser. 8, Vol. 6, Pl. 56. 


The same field as shown in fig. 1, after several hours in an air oven at 98°c, x 330. 
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Fig. 3 
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The field shown in fig. 2 after 60 » had been removed by grinding and the new 
surface washed with distilled water. x 330. 
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The field as shown in fig. 3 after 4 hours at 98°c. x 330. 
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A row of etch pits on the (010) face showing a Geen p ally s in the region of a 
cleavage step. x 300. 
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Fig. 2 (continued) 
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‘Dislocation loops produced. by low temperature pyramidal glide. The suc- 
cessive stages in the process of splitting into circular loops were taken 
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at 20sec intervals at an estimated temperature of —40° co. 
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99:9999°, Aluminium foil, strained in the electron microscope. The axis of 


tension is vertical in the micrograph. Dislocation movement is from 
right to left. x 10000. 


M. J. MAKIN et al. Phil. Mag. Ser. 8, Vol. 6, Pl. 63. 


Dislocation loops in copper irradiated at —195°c (3-1 x 10" fast n cm) 
and examined at —50°c. (x 100 000.) 
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Irradiated copper (2 x 1018 fast n em~?) annealed at 306°c for 63 min showing. 
small defects and large loops. Foil thinned after annealing. (x 100 000.) 


Fig. 4 


As in fig. 3, but showing a zone (3500 A wide) free of large loops along a grait 
boundary. The small defects show no denuded zone. e 100 000.) o j 


